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Abstract. There has been substantial progress in recent years in the quantitative understanding of 
the nonequiHbrium time evolution of quantum fields. Important topical applications, in particular in 
high energy particle physics and cosmology, involve dynamics of quantum fields far away from the 
ground state or thermal equilibrium. In these cases, standard approaches based on small deviations 
from equilibrium, or on a sufficient homogeneity in time underlying kinetic descriptions, are not 
applicable. A particular challenge is to connect the far-from-equilibrium dynamics at early times 
with the approach to thermal equilibrium at late times. Understanding the "link" between the early- 
and the late-time behavior of quantum fields is crucial for a wide range of phenomena. For the 
first time questions such as the explosive particle production at the end of the inflationary universe, 
including the subsequent process of thermalization, can be addressed in quantum field theory 
from first principles. The progress in this field is based on efficient functional integral techniques, 
so-called «-particle irreducible effective actions, for which powerful nonperturbative approximation 
schemes are available. Here we give an introduction to these techniques and show how they can 
be applied in practice. Though we focus on particle physics and cosmology applications, we 
emphasize that these techniques can be equally applied to other nonequilibrium phenomena in 
complex many body systems. 

Based on summer school lectures presented at HADRON-RANP 2004, March 28 - April 3, 2004, 
Rio de Janeiro, Brazil and at QUANTUM FIELDS IN AND OUT OF EQUILIBRIUM, September 
23 - 27, 2003, Bielefeld, Germany. 
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1. MOTIVATIONS AND OVERVIEW 

Cosmology is time evolution. Solid theoretical descriptions exist for the temporal and 
spatial local thermal equilibrium related to the time evolution of a homogeneous and 
isotropic metric and for small perturbations of this. Furthermore, there are powerful 
numerical techniques for A''-particle simulations. In contrast, research concerning the 
interplay between fluctuations and the time evolution of fields is still scarce. Their 
interaction is of crucial importance for the generation of density fluctuations during 
the inflationary phase of the early universe and of entropy at the end of inflation. The 
corresponding temperature fluctuations in the cosmic microwave background radiation 
have led to spectacular high-precision measurements of cosmological parameters. It is 
the dynamics of fluctuations which decides the question whether the baryon asymmetry 
was created during a cosmological phase transition. Back reactions of large density 
fluctuations on the evolution of the cosmic scale factor are also possible. 

A frequently employed strategy is to concentrate on classical statistical field theory, 
which can be simulated numerically. It gives important insights when the number of 
field quanta per mode is sufficiently large such that quantum fluctuations are suppressed 
compared to statistical fluctuations. However, classical Ray leigh- Jeans divergences and 
the lack of genuine quantum effects — such as the approach to quantum thermal equi- 
librium characterized by Bose-Einstein or Fermi-Dirac statistics — limit their use. A 
coherent understanding of the time evolution in quantum fleld theory is required — a 
program which has made substantial progress in recent years with the development of 
powerful theoretical techniques. For the first time questions such as the explosive parti- 
cle production at the end of the inflationary universe, including the subsequent process 
of thermalization, can be addressed in quantum field theory from first principles. Ther- 
malization leads to the loss of a substantial part of the information about the conditions 
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in the early universe. The precise understanding of phenomena out of equilibrium play 
therefore a crucial role for our knowledge about the primordial universe. Important ex- 
amples are the density fluctuations, nucleosynthesis or baryogenesis — the latter being 
responsible for our own existence. 

The abundance of experimental data on matter in extreme conditions from relativistic 
heavy-ion collision experiments, as well as applications in astrophysics and cosmol- 
ogy urge a quantitative understanding of far-from-equilibrium quantum field theory. 
The initial stages of a heavy-ion collision require considering extreme nonequilibrium 
dynamics. Connecting this far-from-equilibrium dynamics at early times with the ap- 
proach to thermal equilibrium at late times is a challenge for theory. The experiments 
seem to indicate early thermalization whereas the present theoretical understanding of 
QCD suggests a much longer thermal equilibration time. To resolve these questions, it 
is important to understand to what "degree" thermalization is required to explain the 
observations. Different quantities effectively thermalize on different time scales and a 
complete thermalization of all quantities may not be necessary. For instance, an approx- 
imately time-independent equation of state p = p{e), characterized by an almost fixed 
relation between pressure p and energy density £, may form very early — even though 
the system is still far from equilibrium! Such prethermalized quantities approximately 
take on their final thermal values already at a time when the occupation numbers of indi- 
vidual momentum modes still show strong deviations from the late-time Bose-Einstein 
or Fermi-Dirac distribution. Prethermalization is a universal far-from-equilibrium phe- 
nomenon which occurs on time scales dramatically shorter than the thermal equilibration 
time. In order to establish such a behavior it is crucial to be able to compare between the 
time scales of prethermalization and thermal equilibration. Approaches based on small 
deviations from equilibrium, or on a sufficient homogeneity in time underlying kinetic 
descriptions, are not applicable in this case to describe the required "link" between the 
early and the late-time behavior. 

A successful description of the dynamics of quantum fields away from equilibrium 
is tightly related to the basic problem of how macroscopic irreversible behavior arises 
from time-reversal invariant dynamics of quantum fields. This is a fundamental question 
with most diverse applications. The basic field theoretical techniques, which are required 
to understand the physics of heavy-ion collision experiments or dynamics in the early 
universe, are equally relevant for instance for the description of the dynamics of Bose- 
Einstein condensates in the laboratory. 

1.1. How to describe nonequilibrium quantum fields from first 

principles? 

There are very few ingredients. Nonequilibrium dynamics requires the specification 
of an initial state at some given time tQ. This may include a density matrix poito) in 
a mixed (Trp^(?o) < 1) or pure state (Trp^(?o) = 1). Nonequilibrium means that the 
initial density matrix does not correspond to a thermal equilibrium density matrix: 

poito) 7^ p^^°^^ with for instance p^^""^ ~ e~^^ for the case of a canonical thermal 
ensemble with inverse temperature /3. In contrast to close-to-equilibrium field theory. 
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the initial density matrix p£)(?o) may deviate substantially from thermal equilibrium. 
This preempts the use of (non-)linear response theory, which is based on sufficiently 
small deviations from equilibrium, or assumptions about the validity of a fluctuation- 
dissipation relation. Since time-translation invariance is explicitly broken at initial times, 
there will be no assumption about a sufficient homogeneity in time underlying effective 
kinetic descriptions. In their range of applicability these properties should come out of 
the calculation. 

Completely equivalent to the specification of the initial density matrix poito) 
is the knowledge of all initial correlation functions: the initial one-point function 
Tr{pz)(?o)^(?o,x)}, two-point function Tr{p£)(ro)^'(?o,x)^'(?o,y)}, three-point function 
etc., where ^(?,x) denotes a Heisenberg field operator. Typically, the "experimental 
setup" requires only knowledge about a few lowest correlation functions at time 
whereas complicated higher correlation functions often build up at later times. The 
question that nonequilibrium quantum field theory addresses concerns the behavior of 
these correlation functions for t > Iq from which one can extract the time evolution of 
other quantities such as occupation numbers. This is depicted schematically in the figure 
below. Of course, one could equally evolve the density matrix in time and compute 
observable quantities such as correlations from it at some later time. However, this 
is in general much less efficient than directly expressing the dynamics in terms of 
correlations, which we will do here (cf. Sec. 6). 

Once the nonequilibrium initial state is specified, the time-evolution is completely 
determined by the Hamiltonian or, equivalently, the dynamics can be described in terms 
of a functional path integral with the classical action S. From the latter one obtains the 
effective action T, which is the generating functional for all correlation functions of 
the quantum theory (cf. Sec. 2). There are no further ingredients involved concerning 
the dynamics than what is known from standard vacuum quantum field theory. Here we 
consider closed systems without coupling to an external heat bath or external fields. 
There is no course graining or averaging involved in the dynamics. In this respect, the 
analysis is very different from irreversible phenomenological approaches. The fact that 
the dynamics is obtained from an effective action automatically guarantees a number 
of crucial properties for the time evolution. The most important one is conservation of 
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energy. The analogue in classical mechanics is well known: if the equations of motion 
can be derived from a given action then they will not admit any friction term without 
further approximations. 

It should be stressed that during the nonequilibrium time evolution there is no loss of 
information in any strict sense. The important process of thermalization is a nontrivial 
question in a calculation from first principles. Thermal equilibrium keeps no memory 
about the time history except for the values of a few conserved charges. Equilibrium 
is time-translation invariant and cannot be reached from a nonequilibrium evolution on 
a fundamental level. It is striking that we will observe below that the evolution can 
go very closely towards thermal equilibrium results without ever deviating from them 
for accessible times. The observed effective loss of details about the initial conditions 
can mimic very accurately the irreversible dynamics obtained from phenomenological 
descriptions in their range of applicability. 



1.1.1. Standard approximation methods fail out of equilibrium 

For out-of-equilibrium calculations there are additional complications which do not 
appear in vacuum or thermal equilibrium.' The first new aspect concerns secularity: 
The perturbative time evolution suffers from the presence of spurious, so-called secu- 
lar terms, which grow with time and invalidate the expansion even in the presence of a 
weak coupling. Here it is important to note that the very same problem appears as well 
for nonperturbative approximation schemes such as standard l/N expansions, where 
denotes the number of field components.^ Typically, secularity is a not a very difficult 
problem and for a given approximation there can be various ways to resolve it. There 
is a requirement, however, which poses very strong restrictions on the possible approxi- 
mations: Universality, i.e. the insensitivity of the late-time behavior to the details of the 
initial conditions. If thermal equilibrium is approached then the late-time result is uni- 
versal in the sense that it becomes uniquely determined by the energy density and further 
conserved charges. To implement the necessary nonlinear dynamics is demanding. Both 
requirements of a non-secular and universal behavior can indeed be fulfilled using ef- 
ficient functional integral techniques: so-called n-particle irreducible effective actions, 
for which powerful nonperturbative approximation schemes are available. They provide 
a practical means to describe far-from-equilibrium dynamics as well as thermalization 
from first principles. 

An illustrative example from classical mechanics. It is instructive to consider for a 
moment the simple example of the time evolution of a classical anharmonic oscillator. 



This does not concern restrictions to sufficiently small deviations from thermal equilibrium, such as 
described in terms of (non-)linear response theory, which only involve thermal equilibrium correlators in 
real time. 

^ Note that restrictions to mean-field type approximations such as leading-order large-A^ are insufficient. 
They suffer from the presence of an infinite number of spurious conserved quantities, and are known 
to fail to describe thermalization. Secularity enters the required next-to-leading order corrections and 
beyond. (Cf. Sec. 4.) 
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This exercise will explain some general problems, which one encounters using perturba- 
tive techniques for initial-value problems, and indicates how to resolve them. Below we 
will consider an illustrative "translation" of the outcome to the case of nonequilibrium 
quantum fields. One of the benefits will be that important concepts such as n-particle 
irreducible effective actions appear here in a very intuitive way, before they will be thor- 
oughly discussed in later sections. 

Our damped oscillator with time-dependent amplitude y{t) is characterized by an 
infinite number of anharmonic terms: 

y+y = -£y-{£y)^~{£yy-{£yy , e<l (1) 

Here each dot denotes a derivative with respect to time t. The presence of anharmonic 
terms to arbitrarily high order is reminiscent of the situation in quantum field theory 
(QFT), where the presence of quantum fluctuations can induce self-interactions to high 
powers in the field. Consider the ideal case for perturbative estimates, i.e. the presence 
of a small parameter £ which suppresses the contributions from higher powers in y. 
The example is chosen such that it can be easily solved without further approximations 
numerically by summing the geometric series with L«=3(£)')" = " For 

the above second-order differential equation the intial- value problem is defined by giving 
the amplitude and its first derivative at initial time t = 0. We employ e.g. y{0) = 1, 
y{0) = — e/2 and consider the time evolution for t >0. One finds the expected damped 
oscillator behavior displayed below for e = 0. 1: 




The question is whether an accurate approximate description of the full);(?) is possible 
if contributions from higher powers of £ are neglected? In view of the presence of the 
small expansion parameter this is certainly possible, however, standard perturbation 
theory fails dramatically to provide a good description. A perturbative expansion of y 
in £, 

Jpert(0 = yo{t) + eyi (t) + yj2(0 + ^(C^) , (2) 

yields the standard hierarchy of equations: 

yo{t)+yoit) = o , yi{t)+yi{t)= -yo , ... (3) 
These can be iteratively solved as 

yo{t) = y' + c.c. , yi{t) = -y't + c.c. , ... (4) 
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Doing this to second order in e one finds: 



>'pert(0 = 




Beyond the lowest order, the solution contains secular terms which grow with powers 
of the time t and one arrives at the important conclusion: 

The perturbative expansion is only valid for et 



We now consider an alternative ("nPI type") expansion: For this we classify the terms 
of the equation of motion (1) according to powers of the small parameter e. The lowest 
order takes into account all terms of the equation to order £*^, which to this order gives 
the same as in perturbation theory: 

>SU^S = ^ >S(') = i^"-Hc.c. (6) 
The next order takes into account all terms of (1) to order e: 

^•(2) , ^(2) _ _ .(2) (2) , . _ }_ it,/T^-et/2 , (J. 

Since there are no £^ terms appearing in (1), this also agrees with the second order of 
the expansion. The next order would take into account all terms of the equation to order 
and so on. This provides a systematic approximation procedure in terms of powers 
of the small parameter £. However, in contrast to the perturbative expansion the "nPI 
type" expansion turns out to be non-secular in time. One explicitly observes that both 

the lowest order approximation and y^pi(0 do not exceed ^{e^) for all times. Therefore, 
the result (7) provides an approximation to the equation (1) up to quantitative corrections 
of order at all times. The already very good agreement with the exact result can be 
checked explicitly. 

It is instructive to compare the "nPI type" solution to the perturbative one by expand- 
ing (7) in powers of e: 

yS(0 = ^e'^v^-"/^+c.c. 

- ^^'■^(l-f^ + y[r-.Y])-f^(£3)-Fc.c. 

= ypertW + ^(e^)+C.C. 

We conclude that 

• the "wPI" result corresponds to the perturbative one up to the order of approxima- 
tion (here up to i?{e^) corrections). 

• the "nPI" result resums all secular terms to infinite perturbative order in £. 
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Stated differently: infinite summation of perturbative orders is required to obtain a 
uniform approximation to the exact solution, i.e. that the error stays of a given order 
at all times. 

What was the reason for the "success" of this alternative expansion? It is based on a 
standard procedure for differential equations in order to enlarge the convergence radius 
of an expansion. In the differential equation one neglects contributions from higher 
powers of the expansion parameter. However, in contrast to perturbation theory one 
does not expand in addition the variable y. This procedure is sometimes called "self- 
consistent" since at each order in the approximation the dynamics is solely expressed 
in terms of the dynamical degree of freedom itself. For instance, in Eq. (7) both on the 
left and on the right hand side of the differential equation appear derivatives of the same 

(2) 

variable y^pj. In contrast, in the perturbative hierarchy (3) the presence of "external" 
oscillating "source terms" such as jq driving y\ leads to secular behavior. The procedure 
of the "nPI type" expansion is simple and turns out to be very efficient to achieve 
non-secular behavior: For given dynamical degrees of freedom truncate the dynamics 
according to powers of a small expansion parameter. Of course, the choice of the degrees 
of freedom is always based on physics. The price to be paid for this expansion scheme 
is that at some order one necessarily has to solve nonlinear equations without further 
approximations. However, as we will see next, it is precisely the nonlinearity which is 
required to be able to obtain universality in the sense mentioned above. 

In order to illustrate properly the aspect of universality the single oscillator exam- 
ple (1) is too simple. In order to include more degrees of freedom we add a "three- 
momentum" label to the variables and consider: 

h^l [(l+3'p)(l+)'q)3'k3'p-q-k - Jpyq(l+yk)(l+yp-q-k)] ■ (8) 
Jqk 

Here the integrals on the r.h.s. involve momenta q and k above some suitable value, and 
one observes that the time derivative of y^ is proportional to a "gain" and a "loss" term. 
This is like in a Boltzmann equation, which describes the rate of scattering of particles 
into the state with momentum p minus the rate of scattering out of that momentum 
(cf. Sec. 4.1.3). If the nonlinear equation (8) is solved without further approximations 
then the late-time solution is given by the well-known result 

1 

^ eiSdpl-M) _ 1 ' 

with real parameters /3 and /i, which for the Boltzmann equation are identified as tem- 
perature and chemical potential respectively. The form of the solution (9) is universal, 
i.e. completely independent of the details of the initial condition for the solution of (8). 
We emphasize that the Boltzmann type equation (8) is "self-consistent" in the sense 
mentioned above. In contrast, a non-"self-consistent" approximation will not show the 
desired universality in general. For instance, consider (8) in a linearized approximation 

jp = (H-jp)cJp"-jpO^, (10) 
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with time-independent ~ /qk(l + Jq(0))3'k(0)3'p-q-k(0) and equivalently for o^. 
Eq. (10) has the solution 



+ 



ypio)-^ 



e 'p . 



(11) 



with = Op — <7p. For late times ^ 1 this always depends on the chosen initial 
^p/lp of course, only useful if they are chosen to be already the equilibrium 

values. 



7.7.2. n-Particle irreducible expansions: universality and non-secularity 

One should not misunderstand the following illustrative "translation" of the above 
mechanics examples to QFT. It is intended to give intuitive insight. The questions of 
secularity and universality are of course settled by actual nonequilibrium calculations in 
QFT, which is the topic of the main body of this text. For a moment, let us transcribe the 
above results to the language of QFT. The amplitude y of the above oscillator example 
plays the role of the one-point function, i.e. the field expectation value or macroscopic 
field. The amplitude squared plays the role of the two-point function or propagator, the 
cubic amplitude that of the three-point function or proper three-vertex, etc.: 

y{t) 0(x) = (<J>(x)) 

(macroscopic field) 

y\t) " ^ " G2{x,y) = {T^x)^y)) - H^)^{y) ^ G{x,y) 

(propagator) 

y^{t) " ^ " G3 (jc, J, z) or, with G3 ~ GGGV3 : 

V3 {x, y, z) (proper three- vertex) 

y'^{t) V4{x,y,z,w) (proper 4- vertex) 



where the symbol T denotes time- ordering. The knowledge of all n-point functions 
provides a full description of the quantum theory. In contrast to the classical example, 
the information contained e.g. in the two-point function cannot be recovered from the 
one-point function: G 7^ etc. In QFT all n-point functions (j), G, V3, V4 . . . constitute 
the set of dynamical "degrees of freedom". 

The equations of motions for all n-point functions are conveniently encoded in a 
generating functional or effective action. There are different functional representa- 
tions of the effective action. The so-called n-particle irreducible (nPI) ejfective action 
r[^, G, V3, . . . , Vn] is expressed in terms of 0, G, V3, . . . , Vn and is particularly efficient for 
a description of suitable approximation schemes in nonequilibrium quantum field the- 
ory. It determines the equations of motion for G, V3, . . . , V„ by first-order functional 
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derivatives or so-called stationarity conditions: 



5rUG^v^_ 5rUG^v^_ 5r[fG^i^_ 

5<^ " ' 5G " ' 573 " 

The above "nPI" approximation scheme for the oscillator example reads now: 

• Classify the contributions to r[^, G, Vs, . . .] according to powers of a suitable ex- 
pansion parameter (e.g. coupling/loops, or 1 /N). In terms of the "variables" ^, G, 
V3 . . . there is a definite answer for the approximate effective action to a given order 
of the expansion. 

• The equations of motion for ^, G, V3 ... are then determined from the approximate 
effective action according to (12) without any further assumptions. 

In general this leads to nonlinear integro-differential equations, which typically require 
numerical solution techniques if one wants to provide the "link" between the early and 
the late-time behavior. We emphasize that nonlinearity is a crucial ingredient for late- 
time universality. The reward is that one can compute far-from-equilibrium dynamics as 
well as the subsequent approach to thermal equilibrium in quantum field theory from first 
principles (cf. Sec. 4). In turn this can be used to derive powerful effective descriptions 
and to determine their range of applicability. 

In order to be useful in practice it is crucial to observe that it is not necessary to 
calculate the most general r[^, G, Vs, . . . , V^] for arbitrarily large n. There exists an 
equivalence hierarchy between nPI effective actions. For instance, in the context of a 
loop expansion one has: 

r(iioop)[^] = r("°°p)[<^),D] = ... , 

p(2ioop)[^l ^ r(2i°°p)[(^,D] = r(2i°°p)[<^,D,y3] = 

r(31oop) ^ p(31oop) [(^ ^ £)] ^ r(31°°P) [(^) , D, Vg] = r(31°°P) , D, , V4] = . . . , 



where r*^"^'°°P) denotes the approximation of the respective effective action to n-th 
loop order in the absence of external sources. E.g. for a two-loop approximation all 
nPI descriptions with n > 2 are equivalent and the 2PI effective action captures already 
the complete answer for the self-consistent description up to this order. In contrast, a 
self-consistently complete result to three-loop order requires at least the 3PI effective 
action in general, etc. There are, however, other simplifications which can decrease 
the hierarchy even further such that lower nPI effective actions are often sufficient in 

practice. For instance, for a vanishing macroscopic field one has rg°°P)[(^) = 0,G] = 
p(^3ioop)|-^ = 0, G, V3 = 0,V4] = ... in the absence of sources. Typically the 2PI, 3PI 
or maybe the 4PI effective action captures already the complete answer for the self- 
consistent description to the desired/computationally feasible order of approximation. 
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2. A^-PARTICLE IRREDUCIBLE GENERATING FUNCTIONALS I 



In this section we discuss the construction of the one-particle irreducible (IPI) and the 
two-particle irreducible (2PI) effective action. The latter provides a powerful starting 
point for systematic approximations in nonequilibrium quantum field theory. To be 
specific, we consider first a quantum field theory for a real, A^-component scalar field cpa 
(a = I,.. .,N) with classical action S[(p]:^ 

S[(P] = i^-d^(Pa{x)d^^a{x) - '^%{X)%{X) - ^ {%{x)(Pa{x)f^ • (13) 

Here summation over repeated indices is implied and we use the shorthand notation 
= Jc^ diX^ J d^x with x = {x^,x). The time path ^ will be specified later. 
The generating functional W [J, R] for connected Green's functions in the presence of 

two source terms ~ Ja{x) and ~ Rab{x,y) is given by 



Z[J,R] = exp(?W[y,i?]) 



j exp S[(p\ + j Ja{x)(pa{x) + Rab{x,y)(pa{x)(pb{y) 



.(14) 



We define the macroscopic field and the connected two-point function Gab by varia- 
tion of W in the presence of the source terms: 

Before constructing the 2PI effective action, we consider first the IPI effective action. It 
is obtained by a Legendre transform with respect to the source term which is linear in 
the field, 

T^[<P] = W[J,R] - [ ^^jll^Ux) = W[J,R] - [ Ux)Ja{x) . (16) 

Jx OJa{X) Jx 

We note that: 

1. r^=*^[^] corresponds to the standard IPI effective action. 

2. Fori? 7^ it can be formally viewed as the IPI effective action for a theory governed 
by the modified classical action, 5^ [9], in the presence of a non-constant "mass 
term" ~ Rab{x,y) quadratic in the fields: 

S'^[(p]=S[(p] + l- f Rab{x,y)%{x)(pb{y) ■ (IV) 

^ Jxy 

This can be directly observed from (14) and the fact that the standard IPI effective action 
is obtained from the same defining functional integral in the presence of a linear source 



^ We always work in units where h = c~ \ such that the mass of a particle is equal to its rest energy (mc^) 
and also to its inverse Compton wavelength (mc/H). 
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term only. As a consequence, it is straightforward to recover for r^[0] all "textbook" 
relations for the IPI effective action taking into account R. For instance, r^[^] to one- 
loop order is given by 

pi?(iloop)[^| = 5^[,^] + lTrln[Goi((^))-/i?], (18) 

which is the familiar result for the IPX effective action with 5[(^] — > S^[(j)] and Gq ^(0) — > 
Gq ^ (^) — iR. Similarly, one finds that the exact inverse propagator is obtained by second 
functional field differentiation. 



(19) 



= ^ [^oiib i^^y)- iRab y) - ^ab 3^) 

Here the classical inverse propagator /Gq ^) = 5^S[^]/ 5^a{x)5^b{y) reads 

^Go_ab i^^r^<l>) = -(^^^ + '^^ + -^^^- (■^) (^) ^ Sab Six-y) 

X 

-jj;^(l>a{x)(l)b{x)5{x-y), (20) 

and E^^(;c, j) is the proper self-energy, to which only one-particle irreducible Feynman 
diagrams contribute, i.e. diagrams which cannot be separated by opening one line. These 
relations will be used below. 

We now perform a further Legendre transform of r^[^] with respect to the source R 
in order to arrive at the 2PI effective action: 



mJ,R] - ^i^^^^)^^^y) + Gab{x,y) 



5Rab{x,y) 2 

= r^[(l>]~^-I^^Rab{x,y)Ux)(l>b{y)-lT:rRG. (21) 

Here we have used that the relation between and J is 7?-dependent, i.e. inverting 
= dW[J,R]/5J yields J = 7^((^)). From the above definition of r^[^] one therefore 
finds 

5r''[(^] 5W[IR] ^ f 5W[J,R] 5Uz) f. , . 57,(z) 



5Rab{x,y) 5Rab{x,y) Jz SJdz) 5Rab{x,y) Jz 5Rab{x,y) 

= iWL. (22) 

dRab{x,y) 
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Of course, the two subsequent Legendre transforms which have been used to arrive at 
(21) agree with a simultaneous Legendre transform of W [/, R] with respect to both source 
terms: 

Jx OJa{x) Jxy ORab{x,y) 

= W[J,R]- ju^)Ja{x)-^-j^Rab{x.y)Ux)Uy)-\^^GR (23) 
From this one directly observes the stationarity conditions: 



-Ja{x) - j Rat{x,y)(l)b{y) , (24) 

(25) 



d^a{x) 



SGahix.y) 2 

which give the equations of motion for and G in the absence of the sources, i.e. 7 = 
and R = 0. 

To get familiar with Eq. (21) or (23), we may directly calculate r[^, G] to one-loop 
order using the above results for r^[<|>]. Plugging (18) into (21) one finds to this order 



r[(^,G] -S[^] + ^Tr\n[GQ\^)-iR]-lTTRG. (26) 
If we set = Gq ^ (0) — iR then we can write 



r[0,G] ~ 5[0] + -TrlnG" ^ + -Tr Go ^G + const, (27) 

with Tr G ^ G ~ Tr 1 = const. To verify this one can check from the stationarity condition 
(25) that indeed to this order 

To go beyond this order it is convenient to write the exact r[0, G] as the one-loop type 
expression (27) and a 'rest': 



r[0, G] = S[0] + ^TrlnG^i + ^TrGg \(^) G-fFjf^, G] + const 



(28) 



Here we have added an irrelevant constant which can be adjusted for normalization. To 
get an understanding of the 'rest' term r2[0, G] we vary this expression with respect to 
G, which yields 

G~b\x,y) = Goli,{x,y;^)-iRab{x,y)-Lab{x,y;^,G), (29) 
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where we have written 



gr2[0,G] 



Comparing with the exact expression (19) one observes that 

Lat{x,y;^,G) =Ef^(x,y;(^)) , 



(30) 



(31) 



which relates the functional G-derivative of r2[0, G] to the proper self-energy. As men- 
tioned above, to the proper self-energy only IPI diagrams contribute with propagator 
lines associated to the effective classical propagator (Gq ^ — iR)^^ . The mapping be- 
tween Ea^(x, ^, G) and E^^(jc,);; ^) is provided by (29), which can be used to express 
the full propagator G as an infinite series in terms of the classical propagator Go and E: 



(Go 1 -/7?)-i + (Go 1 -/7?)-iE(Go 1 -/i?)-i 
+(Go-^-/i?)"^E(Go^-/i?)"^E(Go^-/^)"^ + ... 



(32) 



where we employ an obvious matrix notation. As a consequence there is a direct corre- 
spondence between a IPI diagram with propagator lines associated to G and an infinite 
set of IPI diagrams with propagator lines associated to the classical propagator Gq. Be- 
low we will consider some explicit examples. Most importantly, from the fact that to 
E((^, G) only IPI diagrams contribute one can conclude that 

• r2[0, G] only contains contributions from two-particle irreducible (2PI) diagrams. 

A diagram is said to be two-particle irreducible if it does not become disconnected 
by opening two lines. Suppose r2[^,G] had a two-particle reducible contribution. 
The latter could be written as f GGf', where GG denotes in a matrix notation two 
propagator lines connecting two parts f and f ' of a diagram. Then E((^, G) would contain 
a contribution of the form f Gf ' since it is given by a derivative of r2 with respect to G. 
Such a structure is one-particle reducible and cannot occur for the proper self-energy. 
Therefore two-particle reducible contributions to r2[^, G] are absent. 

Diagrammatically, the graphs contributing to E((^,G) are obtained by opening one 
propagator line in graphs contributing to r2[0, G]. This is exemplified for a two- and a 
three-loop graph for r2 and the corresponding self-energy graphs: 



E - 5r2/5G: 



Because of (32) each 2PI diagram corresponds to infinite series of IPI diagrams, and the 
above 2PI two- and three-loop diagrams contain e.g. the full series of so-called "daisies" 
and "ladder" resummations: 





Introduction to Nonequilibrium Quantum Field Theory 



15 



2.1. Loop or coupling expansion of the 2PI effective action 



Loop or coupling expansions of the 2PI effective action proceed along the same lines 
as the corresponding expansions for the standard IPI effective action, with the only 
difference that 

• the full propagator G is associated to propagator lines of a diagram 

• and only 2PI graphs are kept. 

To be specific, for the considered A'^-component scalar field theory the diagrams are 
constructed in the standard way from the effective interaction 




which is obtained from the classical action (13) by shifting (Pa{x) ^a{^) + <Pa(^) 
and collecting all terms cubic and quartic in the fluctuating field (Pa{x). As for the 
IPI effective action, in addition to the quartic interaction there is an ejfective cubic 
interaction for non-vanishing field expectation value, i.e. 0a (x) ^ 0. 

Since r[0, G] is a functional, which associates a number to the fields and G, only 
closed loop diagrams can appear. We consider for a moment the real scalar theory for 
A^^ = 1. To lowest order one has r2[0, G] =0 and we recover the one-loop result given in 
Eq. (27). At two-loop order there are two contributions 

rf°°P)[fG] = -i3[-i^)jG\x,x) 

-i6 ^ £ ( - i^Hx)) ( - i^m) G\x,y) , (34) 

where we have made explicit the different factors coming from the overall — z in the 
defining functional integral for r[0, G] (cf. Eq. (14)), the combinatorics and the vertices. 
The 2PI loop expansion exhibits of course much less topologically distinct diagrams 
than the respective IPI expansion. For instance, in the symmetric phase (0 = 0) one 
finds that up to four loops only a single diagram contributes at each order. At fifth order 
there are two distinct diagrams which are shown along with the lower-loop graphs in 
Fig. 1. For later discussion, we give here also the results up to five loops for general A'^. 
For = by virtue of 0{N) rotations the propagator can be taken to be diagonal: 

Gab{x,y) = G{x,y)5ah- (35) 
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FIGURE 1. Topologically distinct diagrams in the 2PI loop expansion up to five-loop order for (p =0. 
The suppressed prefactors are given in Eq. (37). 



One finds to five-loop order: 



(0 



p(5l00p) r „n 

1=2 



<2) 



8 3 



48 3A^ 



xy 



.(4) _ XHN + 2){N + ^) f ^2, .^2, .^2 



48 21N^ 



I G^{x,y)G^{x,z)G^{z,y) 

Jxyz 



.(5) 



iX^{N + 2){N^ + 6N + 20) f , . 



(36) 



128 



81A^3 



J xyzw 



{N + 2){5N + 22) f . 



32 81A^3 



G {x,y)G{x,z)G{x,w)G {z,w)G{y,z)G{y,w) . 



xyzw 



2.2. Renormalization 



The 2PI resummed effective action r[0, G(0)] is defined in a standard way by suitable 
regularization, as e.g. lattice regularization or dimensional regularization, and renormal- 
ization conditions which specify the field theory to be considered. We employ renormal- 
ization conditions for the two-point function, F^^), and four-point function, r^'^\ given 
by 



r('\x,y) ^ 
r^'\x,y,z,w) ^ 



52r[0,G(0)] 



5(l>{x)5(j>{y) 



5(l){x)5(l){y)5(l){z)5(l){w) 



(37) 
(38) 



where we consider a one-component field for notational simplicity. Without loss of 
generality we use renormalization conditions for = which in Fourier space read: 



-m 



i?5 



■A 

dp' 



r(2)( 



P )\p=0 



z2r(4)( 



Pl,P2,P3)\pi=p2=P3=0 



-1, 



(39) 
(40) 
(41) 
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with the wave function renormalization Here the renormalized mass parameter 
corresponds to the inverse correlation length. The physical four-vertex at zero momen- 
tum is given by Xr. 



2.2.1. 2PI renormalization scheme to order 

In order to impose the renormalization conditions (39)-(41) one first has to calculate 
the solution for the two-point field G(^) for 0=0, which encodes the resummation and 
which is obtained from the stationarity condition for the 2PI effective action (25). 

The renormalized field is 

^R = Z-'I^^. (42) 

It is convenient to introduce the counterterms relating the bare and renormalized vari- 
ables in a standard way with 

Zm^^m^ + W, + , dZ^Z-l, (43) 

and we write 

G{(I))=ZGr{(I)r). (44) 
In terms of the renormalized quantities the classical action (13) reads 

S = l^(^^d^(l>Rd^(l>R-^-mlf^-^(l>^^+^-^^ . (45) 

Similarly, one can write for the one-loop part TrlnG^ = TrlnG^^ up to an irrelevant 
constant and with G = G(0): 

1 f 

^TrGoi((^))G((^)) = --J^{n, + ml + 5Z,n, + 5ml)GRix,y;^R)l^y 

-^^^^ jfR{x)GR{x,x- <^r) . (46) 

Here dZ\ , dm^ and dX\ denote the same counterterms as introduced in (43), however, 
approximated to the given order. To express T2 in terms of renormalized quantities it is 
useful to note the identity 

T2[^,G{mx='^2[^R^GR{^R)]\x,+8X, (47) 

which simply follows from the standard relation between the number of vertices, lines 
and fields by counting factors of Z. Therefore, one can replace in T2 the bare field and 



For the Fourier transform of the «-th derivative one has 
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propagator by the renormalized ones if one replaces bare by renormalized vertices as 
well. We emphasize that mass and wavefunction renormalization counterterms, 5Z and 
5m^, do not appear explicitly in Yi- The counterterms in the classical action (45), in 
the one-loop term (46) and beyond one-loop contained in have to be calculated for a 
given approximation of Yi- For an explicit example, we consider here the 2PI effective 
action to order with 

r2[0i?,G/?(0i?)] = -^Xr j G\{x,x\^R)+i^Xl j ^R{x)G\{x,y;^R)^R{y) 

where the last term contains the respective coupling counterterm at two-loop. There are 
no three-loop counterterms since the divergences arising from the three-loop contribu- 
tion in (48) are taken into account by the lower counterterms. 

One first has to calculate the solution Gr{^r) obtained from the stationarity condition 
(25) for the 2PI effective action. For this one has to impose the same renormalization 
condition as for the propagator (39) in Fourier space: 

iG^\p^'^ ^r)\p=Q,^r=o = -ml , (49) 
for given finite renormalized "four-point" field 

5^iG]^\x,y;^R) 



VR{x,y;z,w) 



d(^R{z)(^R{w 

For the above approximation we note the identity 

5^Y[^r,Gr{^r)] 



(50) 



5^R{x)^R{y) 

for 



= zG^ {x,y\^R)\^i^=Q (51) 



8Z = dZ\ , dm = difii , 8X1 = 8X2, (52) 

such that (49) for Gr is trivially fulfilled because of (39). In contrast to the exact theory, 
for the 2PI effective action to order a similar identity does not connect the proper 
four- vertex with Vr.^ Here the respective condition for the four-point field Vr in Fourier 
space reads 

VR{puP2,P3)\pi=p2=pi=0 = -^R- (53) 

Note that this has to be the same than for the four-vertex (41). For the universality 
class of the 0^ theory there are only two independent input parameters, which we 



We emphasize that for more general approximations the equation (51) may only be valid up to higher 
order corrections as well. This is a typical property of self-consistent resummations, and it does not affect 
the renormalizability of the theory. In this case the proper renormalization procedure still involves, in 
particular, the conditions (49) and (53). 



Introduction to Nonequilibrium Quantum Field Theory 



19 



take to be ruR and Xr, and for the exact theory Vr and the four-vertex agree. The 
renormalization conditions (39)-(41) for the propagator and four- vertex, together with 
the scheme (49)-(53) provides an efficient fixing of all the above counter terms. In 
particular, it can be very conveniently implemented numerically, which turns out to be 
crucial for calculations beyond order Xr. 

We emphasize that the approximation (48) for the 2PI effective action can only be ex- 
pected to be valid for sufficiently small (^r <^ iriR/gR. If the latter is not fulfilled there are 
additional ^(A|) contributions at three-loop ~ and ~ A^^^- This approximation 
should therefore not be used to study the theory in the spontaneously broken phase or 
near the critical temperature of the second-order phase transition. Quantitative studies of 
the latter can be performed in practice using 1/A^ expansion of the 2PI effective action 



2.2.2. Renormalized equations for the two- and four-point functions 

From the 2PI effective action to order we find with 5Xi = 8X2 from (52) for the 
two-point function: 



According to (51) this expression coincides with the one for the propagator 
5^T[^R, Gr{^r)\/ 5^R{x)5^R{y) at ^r = 0. It is straightforward to verify this using 



which is valid since the variation of Gr{^r) with ^r does not contribute due to the 
stationarity condition (25). The four-point field (50) in this approximation is given by 



to NLO. 



- il+5Zi)n^ + ml + 5ml 

+^{XR + 5Xi){GR{x,x-(^R) + (^^{x))jd{x- y) 

+^XiGl{x,y;^R)M^)My) + ^XlGl{x,y;^R) . (54) 



dr[^R,GRi^R)] _ dr[^R,GR] 
5^r{x) ~ 5^r{x) 



(55) 



VR{x,y;z,w) 



— {Xr + dXi)d{x — y)d{x — z)d{x — w) 




(56) 
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Inserting the chain rule formula 



8^R{z)5^Riw) 



GR{x,u;^R)VR{u,v;z,w)GR{v,y;^R) 



(57) 



one observes that Eqs. (56) and (54) form a closed set of equations for the determination 
of the counterterms dZi, 8m\ and 8Xi. Together with (52) one notes that 8X would be 
undetermined from these equations alone. This counterterm is determined by taking into 
account the equation for the physical four-vertex, which is obtained from the above 2PI 
effective action as 



5^T[^r,Gr{<^r)] 



5^R{x)5^R{y)5^R{z)5^R{w) 



<l>R=0 



-{?iR + d?i)d{x-y)d{x-z)d{x-w) 



+VR{x,y-z,w) + VR{x,z;y,w) + VR{x,w;y,z) 



2 5iGj^ {x,y;^R) ^ diGj^ {x,z;^r) ^ 5iGj^ {x,w;^r) 



5Gr{z,w) 



where from (54) one uses the relation 



5GR{y,w) 



SGR{y,z) 



(t>R=0 



(58) 



diGj^\x,y;^R) 



5Gr{z,w) 



--{XR + dX,)5ix-y)5ix-z)dix~w) 



+ '-XlGl {x, y;^R = Q)5{x-z)5{y-w). (59) 



We emphasize that the counterterm 5X plays a crucial role in the broken phase, since 
it is always multiplied by the field 0o and hence it is essential for the determination of 
the effective potential. It is also required in the symmetric phase, in particular, when one 
calculates the momentum-dependent four-vertex using Eq. (58). 

It is instructive to consider for a moment the 2PI effective action to two-loop order 
for which much can be discussed analytically. In this case Z = 1 and the renormalized 
vacuum mass mr of Eq. (39) or, equivalently, of Eq. (49) is given to this order by 



1 



m 

T 

2 



1 niR \\ m 



(60) 



where = m| + dnti and X = Xr-\- dX\ and we have used (52). Here we have employed 
dimensional regularization and evaluated the integral in d = A — £ for Euclidean mo- 
menta k. The bare coupling in the action (13) has been rescaled accordingly, X fi^X, 
and is dimensionless; p.^ = 4;re^^^/i^ and Je denotes Euler's constant. Below we will 
consider lattice regularizations for comparison and to go beyond two-loop order. 
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Similarly, the zero-temperature four-point function resulting from Eq. (53) for the 2PI 
effective action to order Xr is given by 



with (52). We emphasize that the same zero-temperature equation is obtained starting 
from the renormalization condition for the proper four- vertex (41) with 5 A = 35Ai. 
One observes that all counterterms are uniquely fixed by the renormalization procedure 
discussed in the previous section. 

Though dimensional regularization is elegant for analytical computations, it turns 
out that high momentum cutoff regularizations are often more efficient for numerical 
implementations. We will discuss below cutoff regularizations that are obtained by 
formulating the theory on a discrete space-time lattice. In particular, it is often convenient 
to carry out the numerical calculations using unrenormalized equations, or equations 
where only the dominant (quadratically) divergent contributions in the presence of 
scalars are subtracted. In order to obtain results that are insensitive to variations of the 
cutoff, it is then typically sufficient to consider momentum cutoffs that are sufficiently 
large compared to the characteristic energy-momentum scale of the process of interest. 
It should be also stressed that a number of important renormalized quantities such as 
renormalized masses or damping rates can be directly inferred from the oscillation 
frequency or the damping of the time-evolution for the unrenormalized propagator G. 



2.3. 2PI effective action for fermions 

The construction of the 2PI effective action for fermionic fields proceeds along very 
similar lines than for bosons. However, one has to take into account the anti-commuting 
(Grassmann) behavior of the fermion fields. The main differences compared to the 
bosonic case can be already observed from the one-loop part (r2 = 0) of the 2PI effective 
action. For vanishing field expectation values it involves the integrals: 

/(/) 
^xj/^Xj/e'^o = -/ln(detAQ ^) = -/TrlnAg^, 

Bosons: -iln J ^(pe'^^ ^ -iln{detGQ^y^ = ^TrlnGo^ (62) 

where ^'q'^^ = / d'^xd'^y \j/{x)iAQ^{x,y) \]/{y) denotes a fermion action that is bilinear in the 
Grassmann fields. For Dirac fermions with mass m^-^^ the free inverse propagator reads 

iA^'ix^y) = [i^,-m^f^] 5ix-y) , (63) 

where ^ = y^dpi with Dirac matrices 7^ (jU = 0, . . . , 3), and 1// = i//^'/'. For the bosons 5o 
is given by the quadratic part of (13). Comparing the two integrals one observes that the 
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factor 1/2 for the bosonic case is replaced by — 1 for the fermion fields because of their 
anti-commuting property. With this difference, following along the lines of Sec. 2 one 
finds that the 2PI effective action for fermions can be written in complete analogy to (28). 
Accordingly, for the case of vanishing fermion field expectation values, (^) = (^) = 0, 
one has: 

' ' ' ' (64) 



r[A] = -?TrlnA"^-/TrAQ ^A + r2[A] +const 



-1, 



Here r2[A] contains all 2PI diagrams with lines associated to the time ordered propa- 
gator A(x,y) = (r^(x)^(y)). As for the IPI effective action diagrams get an additional 
minus sign from each closed fermion loop. The trace "Tr" includes integration over time 
and spatial coordinates, as well as summation over field indices. 

As for the bosonic case of Eq. (25), the equation of motion for A in absence of external 
sources is obtained by extremizing the effective action: 



dA{x,y) 



0. 



Using (64) this stationarity condition can be written as 
with the proper fermion self -energy: 



E(-«(x,y;A) 



. dr2[A] 
5A{y,x) 



(65) 



(66) 



(67) 



2.3.1. Chiral quark-meson model 

As an example we consider a quantum field theory involving two fermion flavors 
("quarks") coupled in a chirally invariant way to a scalar a-field and a triplet of pseu- 
doscalar "pions" (a = 1 , 2, 3). The classical action reads 

S = J d'^x^^\[ri^\ir+^[d^od^o + d^7r"d^'n"] 

+ ^¥[(y + irs^"^"] V-V{a^ + 7t^)}, (68) 
Nf J 

where = 7r"n". Here denote the standard Pauli matrices and h/Nf is the Yukawa 
coupling. The number of fermion flavors is Nf = 2 and Nj is the number of scalar 
components. The above action is invariant under chiral SUi{2) x SUr{2) ~ 0(4) trans- 
formations. For a quartic scalar self-interaction. 
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this model corresponds to the well known linear a-model, incorporating the chiral 
symmetry of massless two-flavor QCD. 

Here we do not consider the possibility of a spontaneously broken symmetry, thus 
field expectation values vanish. For this model the 2PI effective action is then a func- 
tional of fermion as well as scalar propagators. The scalar fields form an 0(4) vec- 
tor ^a{x) = (cr(x), 7r(x) ) and we denote the full scalar propagator by Gabi^.y) with 
a,Z? = 0, . . . ,3. In addition to its Dirac structure the fermion propagator carries flavor 
indices i,j for Nf = 2 flavors. The 2PI effective action for this Yukawa model is given 
by 

r[G,A] = ^TrlnG"^ + ^TrGo -/TrlnA^^ -/TrAo ^A + TzfG, A] +const, (70) 
with the free scalar and fermion inverse propagators 

iGoi,b y) = -{D, + m^)d{x-y) dab , i^Q^j (x, y) = {i^, - rn^-f^ )S{x-y) dij . (71) 

The equation of motions for the scalar and fermion propagators are obtained from the 
stationary conditions 

Sr[G,A] _^ _ |IWA1 =0. (72) 



5Gab{x,y) ' dAij{x,y) 

The first non-zero order in a loop expansion of r2[G, A] consists of a purely scalar 
contribution corresponding to the two-loop diagram of Fig. (1), as well as a fermion- 
scalar contribution depicted here graphically: 

©solid line: fermion propagator A 
dashed Une: scalar propagator G 

Without loss of generality in the absence of spontaneous chiral symmetry breaking, the 
effective action and its functional derivatives can be evaluated for Gab taken to be the 
unit matrix in C>(4)-space. Similarly, the most general fermion two-point function can 
be taken to be proportional to unity in flavor space and we can write: 

Gab (x, y) = G(jc, y) 5ab , A,-j (jc, y) = A{x, y) 5ij . (73) 
The two-loop approximation then reads 



r 

^ 2 



[0,AJlG,6=G4i,A,7=A5,i - -g- 3 J^G {X,X) 

-ih^^ [ tr[A{x,y)A{y,x)]G{x,y), (74) 

^ Jxy 

where the trace "tr" acts only in Dirac space. According to (30) and (67) the self-energies 
can then be obtained by 

y . . ^. f ^^2\g8„i„A5,j dG{x',y') dab^^2\GS„hMj 

^abix^y) = 2i / ■ — - = 22—2^-7 — (75) 

Jyy oG(y,y) oG"*(x,y) Nj oG{x,y) 
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where we have used G = GatSab/Nj with 5aa = Nj, and equivalently for the fermion 
propagator with da = Nf. 



2.4. Two-particle irreducible expansion 

In this section we discuss a systematic nonperturbative approximation scheme for the 
2PI effective action. It classifies the contributions to the 2PI effective action according 
to their scaling with powers of l/N, where A'^ denotes the number of field components: 

r2[(^),G] = rLO[,^,G]+rNLO[<^,G]+r^NLO[(^),G] + ... 

Each subsequent contribution r2^, F^'^^, r^'^^o down by an additional factor 

of l/N. The importance of an expansion in powers of l/N stems from the fact that 
it provides a controlled expansion parameter that is not based on weak couplings. It 
can be applied to describe physics characterized by nonperturbatively large fluctuations, 
such as encountered near second-order phase transitions in thermal equilibrium, or for 
extreme nonequilibrium phenomena such as parametric resonance. For the latter cases a 
2PI coupling or loop expansion is not applicable. The method can be applied to bosonic 
or fermionic theories alike if a suitable field number parameter is available, and we 
exemplify it here for the case of the scalar 0(A'^) -symmetric theory with classical action 
(13). We comment on its application to the chiral quark- meson model below. 



2.4.1. Classification of diagrams 

We present a simple classification scheme based on O(A^) -invariants which 
parametrize the 2PI diagrams contributing to r[^,G]. The interaction term of the 
classical action in Eq. (13) is written such that S[^] scales proportional to A'^. From the 
fields 0a alone one can construct only one independent invariant under 0{N) rotations, 
which can be taken as tr = = ^a^a ^ N . The minimum 0o of the classical effective 
potential for this theory is given by (j)^ = N{—6m^/X) for negative mass-squared irp- 
and scales proportional to N. Similarly, the trace with respect to the field indices of the 
classical propagator Gq is of order A'^. 

The 2PI effective action is a singlet under 0{N) rotations and parametrized by the two 
fields and Gab- To write down the possible 0{N) invariants, which can be constructed 
from these fields, we note that the number of ^-fields has to be even in order to construct 
an 0(A'^)-singlet. For a compact notation we use {^^)ab = ^a^t- AH functions of and 
G, which are singlets under 0{N), can be built from the irreducible (i.e. nonfactorizable 
in field-index space) invariants 

0^, tr(G") and tr(00G"). (77) 
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FIGURE 2. Graphical representation of the (/)-dependent contributions for T2 = 0. The crosses denote 
field insertions ~ 0„0„ for the left figure, which contributes at leading order, and ^ ^a^b for the right 
figure contributing at next-to-leading order 



We note that for given A'^ only the invariants with n < N axe irreducible — there cannot 
be more independent invariants than fields. We will see below that for lower orders in 
the l/N expansion and for sufficiently large A^^ one has n < A''. In particular, for the next- 
to-leading order approximation one finds that only invariants with n <2 appear, which 
makes the expansion scheme appealing from a practical point of view. 

Since each single graph contributing to r[^,G] is an C>(A'^) -singlet, we can express 
them with the help of the set of invariants in Eq. (77). The factors of A'^ in a given graph 
have two origins: 

• each irreducible invariant is taken to scale proportional to A'^ since it contains exactly 
one trace over the field indices, 

• while each vertex provides a factor of l/N. 

The expression (28) for the 2PI effective action contains, besides the classical action, 
the one-loop contribution proportional to TrlnG^^ -l-TrGg ^(^)G and a nonvanishing 
r2[0, G] if higher loops are taken into account. The one-loop term contains both leading 
order (LO) and next-to-leading order (NLO) contributions. The logarithmic term cor- 
responds, in absence of other terms, simply to the free field effective action and scales 
proportional to the number of field components A^^. To separate the LO and NLO con- 
tributions at the one-loop level consider the second term TrGp ^ (^)G. From the form of 
the classical propagator (20) one observes that it can be decomposed into a term pro- 
portional to tr(G) ~ A^and terms ~ {X/6N) [tr(^^)tr(G) +2tr(^0G)]. This can be seen 
as the sum of two "2PI one-loop graphs" with field insertion ~ ^a^a and ~ ^a^b^ re- 
spectively, as shown in Fig. 2. Counting the factors of A'^ coming from the traces and the 
prefactor, one sees that only the first contributes at LO, while the second one is NLO. 

According to the above rules one draws all topologically distinct 2PI diagrams and 
counts the number of closed lines as well as the number of lines connecting two field in- 
sertions in a diagram following the indices. For instance, the left diagram of Fig. 2 admits 
one line connecting two field insertions and one closed line. In contrast, the right figure 
admits only one line connecting two field insertions by following the indices. Therefore, 
the right graph exhibits one factor of A^^ less and becomes subleading. Similarly, for the 
two-loop graph below one finds: 



a a a b a b 




LO : NLO 
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(trG)VA^~A^ trG^/Nr^N° 

The same can be applied to higher orders. We consider first the contributions to 
r2[0 = 0, G] = r2[G], i.e. for a vanishing field expectation value and discuss ^ 7^ be- 
low. The LO contribution to r2[G] consists of only one two-loop graph, whereas to NLO 
there is an infinite series of contributions which can be analytically summed: 



rLO[G] 



Gaa{x,x)Gbb{x,x) , 



-NLO 



[G] 



-Trln[5(G)]^ 



8{x-y)+i— Gab {x, y) Gab {x, y) . 



B{x,y;G) 

In order to see that (79) with (80) sums the following infinite series of diagrams 



(78) 

(79) 
(80) 




one can expand: 

Trln[5(G)] = J (^^Gab{x,x)Gab{x,x) 



\j ; {^^^"b{x,y)Gab{x,y)^ (^i-^Ga'b'{y,x)Ga'b'{y,x) 



+ 



(81) 



The first term on the r.h.s. of (81) corresponds to the two-loop graph with the index 
structure exhibiting one trace such that the contribution scales as trG^/N ~ A^*^. One 
observes that each additional contribution scales as well proportional to {\xG^ /N)" ~ 
for all n > 2. Thus all terms contribute at the same order. 

The terms appearing in the presence of a nonvanishing field expectation value are 
obtained from the effectively cubic interaction term in (33) for 7^ 0. One first notes 
that there is no ^-dependent graph contributing at LO. To NLO there is again an infinite 
series of diagrams ~ which can be summed: 



T-LO 

^ 2 



[0,G] ^ rLO[G] 



(82) 



T-NLO 
^ 2 



I{x,y;G) 



T-NLO 
^ 2 



6A^ 



[^ = Q^G]+'—j^ I{x,y; G) U^) Gab{x,y) hiy) , (83) 
A f 

Gab{x,y)Gab{x,y)-i — J I{x,z;G)Gab{z,y)Gab{z,y) . (84) 
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The series of terms contained in (83) with (84) corresponds to the diagrams: 




The functions I{x,y; G) and the inverse of B{x,y; G) are closely related by 

B-\x,y- G) = 5{x-y)- il{x,y- G) , (85) 

which follows from convoluting (80) with and using (84). We note that B and / do 
not depend on 0, and r2[0, G] is only quadratic in at NLO. 

The 2PI can be straightforwardly applied to other theories as well. For instance, 
the chiral quark-meson model of Sec. 2.3.1 admits a l/Nf expansion. For the latter 
the LO contribution scales ~ A''^ and is given by a purely scalar two-loop term as in 
Eq. (74). Note that this two-loop contribution contains a LO part ~ A^^ as well as a 

NNLO part ^ N^- At NLO the first fermion contributions to T2 appear, which scale 
~ Nf (cf. Eq. (74)). If no small coupling is available the expansion parameter Nf^l'm 
this case might not be suitable for a quantitative estimate at low orders. See however the 
precision tests of Sec. 5 for the scalar 0{N) model, which exhibit a good convergence 
already for moderate values ofN>2. 



2.4.2. Symmetries and validity ofGoldstone 's theorem 

We emphasize that by construction each order in the 2PI 1/A/^-expansion respects 
0{N) symmetry. In particular, this is crucial for the validity of Ward identities. For the 
case of spontaneous symmetry breaking note that Goldstone's theorem is fulfilled at any 
order in the 2PI l/N expansion. According to (77) the 2PI effective action can be written 
as a function of the 0{N) invariants: 

T[^M='^[fMG")MHGP)\ . (86) 

In the case of spontaneous symmetry breaking one has a constant 7^ and the propa- 
gator can be parametrized as 

Gab{^)=GL{^^)P^h + GT{^^)Pl, (87) 

where = ^a^b/^^ and Pj^^ = 5ab — P^b the longitudinal and transverse projectors 
with respect to the field direction. The 2PI resummed effective action r[^,G(^)] is 
obtained by evaluation at the stationary value (25) for G. The mass matrix can 
then be obtained from 



(88) 



=const 



If r[(^, G(^)] is calculated from (86) and (87) one observes that indeed it depends only 
on one invariant, ^^: r[0, G(^)] = r[^^]. The form of the mass matrix ^ab can now be 
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inferred straightforwardly. To obtain the effective potential U (^^/2), we write 



n0']Ueo„st = ^^+i^(<^>V2), (89) 

where ^d+\ is the d+l dimensional Euclidean volume. The expectation value of the 
field is given by the solution of the stationarity equation (24) which becomes 

M^^*„</'^0. (90) 

d<^a 

where U' = dU /d{^^/2) and similarly for higher derivatives. The mass matrix reads 

- ^Tif/^^ = + = + ^'^"^^'^^ + • (91) 

In the symmetric phase {^a = 0) one finds that all modes have equal mass squared 
•^ab ~ ^'^ab- In the broken phase, with ^ 0, Eq. (90) implies that the mass of the 
transverse modes ~ U' vanishes identically in agreement with Goldstone's theorem. 
Truncations of the 2PI effective action may not show manifestly the presence of massless 
transverse modes if one considers the solution of the stationarity equation (24) for the 
two-point field G encoding the resummation. It is important to realize that the 2PI 
resummed effective action r[^, G(^)] complies fully with the symmetries. 



3. NONEQUILIBRIUM QUANTUM FIELD THEORY 

Out of equilibrium dynamics requires the specification of an initial state. This may 
include a density matrix at a given time poitQ = 0) in a mixed (Trp^(O) < 1) or pure state 
(Trp^(O) = 1). Nonequilibrium means that the initial density matrix does not correspond 

to a thermal equilibrium density matrix: p£)(0) ^ p^^°^^ with for instance p^^"^^ ~ e~^^ 
for the case of a canonical thermal ensemble. Once the initial state is specified, the 
dynamics can be described in terms of a functional path integral with the classical 
action S as employed in the previous sections above. The corresponding nonequilibrium 
ejfective action is the generating functional for all correlation functions with the initial 
correlations determined by po(0). 



3.1. Nonequilibrium generating functional 



All information about nonequilibrium quantum field theory is contained in the 
nonequilibrium generating functional for correlation functions: 

Z[J,R-Pd] = Tr{p^(0)r^e'(^^-^W*W+^^''^(-'^)*W*W)|. (92) 

Here denotes time-ordering along the time path 'io appearing in the source term inte- 
grals with = J'^dx^ J d'^x as specified below. In complete analogy to Eq. (14) we have 
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introduced the generating functional with two source terms, 7 and R, in order to construct 
the corresponding nonequilibrium 2PI effective action below. We suppress field indices 
in the notation, which can be directly recovered from (14). The nonequilibrium cor- 
relation functions, i.e. expectation values of time-ordered products of Heisenberg field 
operators 4>(x), are obtained by functional differentiation. For instance the two-point 
function reads (Trpo = 1): 



Tr{po(0)r^4>W4>();)} = {T^^{xmy)) 



5^Z[J,R;Pd\ 



i5J{x) i5J{y) 



■ (93) 

J=R=0 



There is a simple functional integral representation of Z[J,R;pd] similar to (14). One 
writes the trace as 

Z[J,R;Pd] = |d(p(i)(x)d(p(2)(x)(,p(i)|p^(0)|(p(2)) 

((p(2) |re<i:.^W*W+i{.v^(-^,y)*W*W) 1^(1)) ^ (94) 

where the matrix elements are taken with respect to eigenstates of the Heisenberg field 
operators at intial time, 4>(? = 0,x)|(p('^) = q)^'\x)\q)^'^), i = 1,2. The source-dependent 
matrix element can be expressed in terms of a path integral using standard techniques if 
one considers a finite, closed real-time contour 

11 



Contour time ordering along this path corresponds to usual time ordering along the 
forward piece and antitemporal ordering on the backward piece . Note that any 
time on is considered later than any time on One can then use 

((p(2) |re'(i:v-/W*W+2/./(^')')*W*W) |(p(i)) 

(p(2)(x)=(p(0-,x) 

= y ^'(pe'('^[<Pl+/A-/W'PW+3Lv/^(^')')<PW'PW) ^ (95) 

(p(i)(x)=(p(0+,x) 

which is the same relation as employed to obtain standard path integral expressions for 
vacuum or equilibrium matrix elements. Note that the closed time contour starting from 
the initial time ?o = and ending at t^ is required because the density matrix element 
(?'^^-*|P£)(0)|?>*-^'') is taken on both sides with respect to states at time t^. Using (95) 
in (94) one finds that the expression for Z[/,7?;p/)] directly displays the ingredients 
entering nonequilibrium quantum field theory — the quantum fluctuations described 
by the functional integral with action 5, and the statistical fluctuations encoded in the 



Introduction to Nonequilibrium Quantum Field Theory 



30 



weighted average with the initial-time elements: 

(pi'-) 

Z[J,R;Pd] =yd(p(i)d(p(2)((p(i)|p^(0)|<p(2)) J ^',pe<5['P]+X.^W<PW+3Xv,^(^>3')<PW<PW) 

<p(i) 



initial conditions quantum dynamics 

Of course, this generating functional can be equally applied to standard vacuum physics. 
In this case the closed time contour just ensures the normalization Z|j=^=o = 1' which 
can simplify calculations. In particular, the density matrix is time independent for the 
vacuum as well as for the thermal equilibrium case. In contrast to thermal equilibrium, 
the nonequilibrium density matrix cannot be interpreted as an evolution operator in 
imaginary time, as is possible for instance for the canonical p^'^i) ~ e^^^ . 

In the literature formulations of closed time path generating functionals often exhibit 
an infinite time interval ] — oo, oo[. In this case the number of field labels has to be doubled 
to distinguish the fields on the underlying closed contour. Causality implies that for 
any n-point function with finite time arguments contributions of an infinite time path 
cancel for times exceeding the largest time argument of the n-point function. To avoid 
unnecessary cancellations of infinite time path contributions we always consider finite 
time paths. The largest time of the path is kept as a parameter and is evolved in the time 
evolution equations as described below. We stress that the initial time of the path has to 
be finite — a system that can thermalize will be already in equilibrium at any finite time 
if the initial time is sent to the infinite past. 



3.2. Initial conditions 

To understand in more detail how the initial density matrix enters calculations, we 
consider first the example of a Gaussian density matrix whose most general form can be 
written as 



{(p(l)|Pz,(0)|<p(2)) 

1 



exp<{/(^((p(i)-(p(2); 



(7^ + 1 



;(p(i)-0)2+((p(2)-0y 



(97) 



where we neglect the spatial dependencies for a moment. Note that for homogeneous 
field expectation values taking into account spatial dependencies simply amounts to 
adding a momentum label in Fourier space (cf. also Sec. 3.4.3). In order to see that 
this is the most general Gaussian density matrix, we first note that (97) is equivalent to 
the following set of initial conditions for one- and two-point functions: 



(^) = Tr{po(0)4.(0}|,=o . =Tr{pz)(0)<9,4.(0}|,=o 



(98) 
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^2 = Tr{pz)(0)4>W4>(0} 1,^,^0 

^7] = ^-Tr{poiO){dMtm')+mdMn)}\r=,=o-<l'<l'^ (99) 

^' + ^ = Tr{pz)(0)<9,4.(0<9,'4.(O}|,=,,=0-<^<^- 

It is straightforward to check explicitly the equivalence between the initial density matrix 
and the initial conditions for the correlators: 



Trpz)(0) 



d(p {(p\Pd{0)\(P) 

1 f°° r 1 ^ 

==y_d(pexp|-2^((p-0)2| = l 



Tr{pz,(0)4>(0)} 



^27t^2 J- 
1 



27t^2 J- 



1 ^ j.\2l <P^<P- 



2^2 



(100) 
0, (101) 



etc. Similarly, one finds since only Gaussian integrations appear that all initial n-point 
functions with « > 2 are given in terms of the one- and two-point functions. Note also 
that the anti- symmetrized initial correlator involving the commutator of ^> and is 
not independent because of the field commutation relation. 

The crucial observation is that higher initial time derivatives are not independent as 
can be observed from the exact field equation of motion, which reads for the real scalar 
4>'^-theory: 



(<924.) = -m2(4>)-A(4.3). 



(102) 



Since (4>-^) is given in terms of one- and two-point functions for Gaussian P£)(0) also 
second and higher time derivatives are not independent. We conclude that the most 
general Gaussian density matrix is indeed described by the five parameters appearing 
in (97). In particular, all observable information contained in the density matrix can be 
conveniently expressed in terms of correlation functions and their derivatives (98) and 
(99). 

For further interpretation of the initial conditions we note that 



Trp2(0) 



1 



d<?) / dcp\(p\pD{Q)W){(?'\PD{Om = -- 



(103) 



The latter shows that for C7 > 1 the density matrix describes a mixed state. A pure state 
requires a = 1 or, equivalently, using (99) this condition can be expressed in terms of 
initial-time correlators as 



Tr{pz,(0)4>(0*(O}|,=,=o-<^><^' 
1 



Tr{pz)(0)^,4.(04'*(0}|,=,=o-<^<^ 
^Tr {po(0) {dMtmt')'rmdMt'))]\^=,^o ' 



(104) 
1 

4 ■ 
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In field theory such a Gaussian initial condition is associated with a vanishing initial 
particle number, as will be discussed in Sec. 4. For C7 = 1 the "mixing term" in (97) is 
absent and one obtains a pure-state density matrix of the product form: 

PD{o) = \a){a\ (105) 

with Schrodinger wave function 

In order to go beyond Gaussian initial density matrices one can generalize the above 
example and parametrize the most general density matrix as 

{(P^^^Pd{0)\(P^^^) = J^e'^'^^ , (107) 
with normalization factor o/K and f'^'[(p] expanded in powers of the fields: 

M(p] = ao+ f ai{x)(p{x) + ]- [ a2ix,y)(p{x)(p{y) + ^ [ a3ix,y,z)(p{x)(p{y)(p{z) 

J X ^ J xy J ■ J xyz 

/ aA{x,y,z,w)(p{x)(p{y)(p{z)(p{w) + ... (108) 

Here (p(0+,x) = (p(')(x) and (p(0^,x) = (p^'^\x). We emphasize that (108) employs a 
compact notation: since the density matrix P£)(0) is specified at initial time tQ = 0, all 
time integrals contribute only at time to of the closed time contour. As a consequence, the 
coefficients ai{x), a2{x,y), aT,{x,y^z), ■ ■ .vanish identically for times different than ^o- 
For instance, up to quadratic order one has 



ai(jc)(p(x) = ^|aj^)(x)(p(i)(x) + af)(x)(p(2)(x) 

a2{x,y)cp{xMy) ^ f {a^i'i)(x,y)(p(i)(x)(p(i)(y) + (x,y)(p(i)(x)(p(2)(y) 

xy Jxy 

+a(2' ^) (x, y) (x) (y) + ap) ^(2) ^(2) 1 _ 



Note that ocq is an irrelevant normalization constant. For a physical density matrix the 
other coefficients are of course not arbitrary. Hermiticity implies a[^^ = — ap"**, aj^'^"* = 
— a^'^^* and Ctj^'^"* = —a^'^^*, which can be directly compared to the discussion above. 



3.3. Nonequilibrium 2PI effective action 

Using the parametrization (107) and (108) for the most general initial density matrix, 
one observes that the generating functional (96) introduced above can be written as 

Z[J,R;pd] = j ^(pe't^['Pl+-^t-'W'PW+iM(-'^'>0<PW'PW+Tr/xv^ _ 

(109) 
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Here we have neglected an irrelevant normalization constant and rescaled the sources in 
(96) as J{x) J{x) — a\{x) and R{x,y) R{x,y) — a2{x,y). The sources can therefore 
be conveniently used to absorb the lower linear and quadratic contributions coming from 
the density matrix specifying the initial state. This absorbtion of tti in J and a2 in R 
exploits the fact that the initial density matrix is completely encoded in terms of initial- 
time sources (X\, a2, OCj,, ... for the functional integral. 

The generating functional (109) can be used to describe situations involving arbitrarily 
complex initial density matrices. However, often the initial conditions of an experiment 
may be described by only a few lowest n-point functions. For many practical purposes 
the initial density matrix is well described by a Gaussian one. For instance, the initial 
conditions for the reheating dynamics in the early universe at the end of inflation are de- 
scribed by a Gaussian density matrix to high accuracy. Clearly, the subsequent evolution 
builds up higher correlations which are crucial for the process of thermalization. These 
have to be taken into account by the quantum dynamics, which is not approximated by 
the specification of an initial density matrix! 

From (109) one observes that for Gaussian initial density matrices, for which = 
a4 = . . . = 0, one has 

z[7,i?;p^^'"''^j = Z[J,R]. (110) 

As a consequence, in this case the nonequilibrium generating functional corresponds 
to the 2PI generating functional introduced in Eq. (14) for a closed time path. We can 
therefore directly take over all steps from Sec. 2 in order to construct the nonequilibrium 
2PI effective action and obtain the important result: 



Nonequilibrium 2PI effective action = r[(^, G] with closed time path 



All previous discussions of Sec. 2 remain unchanged except that for nonequilibrium 
the time integrals involve a closed time path! We emphasize again that the use of the 
nonequilibrium 2PI effective action represents no approximation for the dynamics — 
higher irreducible correlations can build up corresponding to a non-Gaussian density 
matrix for times t > to. It only restricts the "experimental" setup described by the 
initial conditions for correlation functions. Non-Gaussian initial density matrices pose 
no problems in principle but require taking into account additional initial-time sources. 
This is most efficiently described in terms of nPI effective actions for n> 2. 



3.4. Exact evolution equations 

Without approximation the 2PI effective action r[^,G] contains the complete infor- 
mation about the quantum theory. The functional representation of the nonequilibrium 
r[0, G] employs a one-point (^) and a two-point field (G), whose physical values have to 
be computed for all times of interest. The equations of motion for these fields are given 
by the stationarity conditions (24) and (25). 

We will consider first the scalar field theory in the symmetric regime, where 
= 0, G] = r[G] is sufficient. We come back to the case of a nonvanishing field 
expectation value as well as to fermionic and gauge fields below. The equation of 
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motion (29) for G reads 

G-\x,y) = Go\x,y)-L{x,y;G)-iR{x,y), (111) 

where the self-energy Z{x,y;G) is given by (30) The form of the equation (111) is 
suitable for boundary value problems as e.g. appear for thermal equilibrium. However, 
nonequilibrium time evolution is an initial value problem. The equation of motion can 
be rewritten as a partial differential equation suitable for initial-value problems by 
convolution with G, using /,G^^(jc,z)G(z,j) = 5{x — y): 

jGQ\x,z)G{z,y)- j\L{x,z) + iR{x,y)]G{z,y) = 5{x-y). (112) 

Here we employ the notation 5{x — y) = d<^{x^ —y^)d{x — y). For the scalar theory the 
classical propagator (cf. Sec. 2) is GQ^{x — y) = i [Dx + m^] 5{x — y) and one obtains the 
evolution equation for the time-ordered propagator. 



(a + m2) G{x,y) + ij^[i:{x,z;G)+iR{x,y)] G{z,y) = -i5{x-y) 



(113) 



3.4.1. Spectral and statistical components 

In the following we introduce a decomposition of the two-point function G into spec- 
tral and statistical components. The corresponding evolution equations for the spectral 
function and statistical propagator are fully equivalent to the evolution equation for G, 
but have a simple physical interpretation. While the spectral function encodes the spec- 
trum of the theory, the statistical propagator gives information about occupation num- 
bers. Loosely speaking, the decomposition makes explicit what states are available and 
how often they are occupied. 

For the real scalar field theory there are two independent real-valued two-point 
functions, which can be associated to the expectation value of the commutator and the 
anti-commutator of two fields, 

commutator: p (jc, j) = z ( [4> (jc) , 4> (j) ] ) , (114) 

anti-commutator: F(x,y) = ^({4>(x),4>(y)}) . (115) 

Here p{x,y) denotes the spectral function and F{x,y) the statistical two-point function. 
The decomposition identity for spectral and statistical components of the propagator 
reads: 



G{x,y)=F{x,y)--p{x,y)sign<^{x -y 



(116) 



The identity is easily understood by making the time-ordering for the propagator ex- 
plicit: 



G{x,y) = ($(x)4.(3;))0^(x°-/) + (4.(j)*(x))0^(/-x«) 
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- ^ i{Mx) My)]) i&A^' -y')-&Ay'-^')) 

sign^(/-/). 

For real scalar fields the real functions obey F{x,y) = F{y,x) and p{x,y) = —p{y,x). 
We note from (114) that the spectral function p encodes the equal-time commutation 
relations 

p(jc,j)|,o^yO = , d^op{x,y)\^o^yO = d{x-y). (117) 

The spectral function is also directly related to the retarded propagator p (jc,);)0(x^ — 
or the advanced one —p{x^y)&{y^ —x^). However, it is important to realize that out 
of equilibrium there are only two independent two-point functions — no more. These 
can be associated to F and p, which has the advantage compared to e.g. an advanced 
propagator that the non-analyticity entering through time-ordering is always explicit. 

To obtain a similar decomposition for the self-energy, we separate E in a "local" and 
"nonlocal" part according to 

E(x, y; G) = -iL^^^ {x; G)d{x-y)+ E(jc, y;G). (118) 

Since E^^) just corresponds to a space-time dependent mass-shift it is convenient for the 
following to introduce the notation 

M^{x;G) =m^ + L^^\x;G). (119) 

To make the time-ordering for the non-local part of the self-energy, 'L{x,y;G), explicit 
we can use the same identity as for the propagator (116) to decompose: 

L{x,y)=LF{x,y)-'-Lp{x,y)sign^{x'^-y'^). (120) 

Though the discussion is given here for a vanishing field expectation value, it should 
be emphasized that the same decompositions (116) and (120) apply also for ^ 7^ 0. The 
r.h.s. of (1 15) would then receive an additional contribution subtracting the disconnected 
part ~ 0^. Equivalently, one can always view (116) as defining F and p from the 
connected propagator G irrespective of the value of ^. 

We emphasize that the equivalent decomposition can be done for fermionic degrees of 
freedom as well. The fermion propagator in terms of spectral and statistical components 
reads in the same way as for bosons: 

A{x,y) = F^f\x,y) - ^p^-^Hx^y) sign^(xO -/) . (121) 

However, in contrast to bosons for fermions the field anti-commutator corresponds to 
the spectral function, 

anti-commutator: p^-^\x,y) = i{{^{x),^'{y)}) , (122) 

commutator: F'^-f\x,y) = ^{['¥{x),^'{y)]) . (123) 
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This can be directly observed from the time-ordered fermion propagator A{x,y) = 
{'¥{x)^'{y))@^{x^ -y^) - {^'{y)'¥{x))@^ {y^ - x^). Here the minus sign is a conse- 
quence of the anti-commutation property for fermionic fields (cf. Sec. 2.3). In analogy 
to the bosonic case, the equal-time anti-commutation relations for the fields are again 
encoded in p(-^) {x,y). For instance, for Dirac fermions one has 

r°p(x,3;)Lo=,o=/5(x-y) (124) 
with the Dirac matrix 'Z', and the two-point functions have the hermiticity properties 

(p(3;,x))^ = -fp{x,y)f, {F{y,x))^ = fF{x,y)f. (125) 
The corresponding decomposition for the fermion self-energy reads :^ 

E(/)(x,y) = i:^/\x,y)-^z'/\x,y)sign^{x'-y^). (126) 



3.4.2. Detour: Thermal equilibrium 

To see the above decomposition in a probably more familiar context, we consider 
for a moment thermal equilibrium. This is done for illustrational purposes only, and 
we emphasize that the notion of an equilibrium temperature is nowhere implemented in 
nonequilibrium quantum field theory. If a nonequilibrium evolution approaches thermal 
equilibrium at late times then this a prediction of the theory and not put in by hand. 

The 2PI effective action in thermal equilibrium is given by the same expression (28) 
if the closed time path is replaced by an imaginary path ^ = [0, — //3]. Here /3 denotes 
the inverse temperature. Since thermal equilibrium is translation invariant, the two-point 
functions depend only on relative coordinates and it is convenient to consider the Fourier 
transforms F('^i\(0,P) and p('^i)(£0,p) with 

P^''\-^y) = J |^e-'--(^°-^V^-p(x-y)^(eq)(^^p) (127) 

and equivalently for the thermal spectral function. 

The periodicity ("KMS") condition characterizing thermal equilibrium for the prop- 
agator in imaginary time is given by G{x,y)\^o^Q = G{x,y)\^^_ip. Employing the de- 
composition identity (116) for the propagator G, one can write the periodicity condition 
as 

FM(ft),p) + ^p^'^\co,p) = e-P'' (f^'^\co,p) - ^p^^'i)(a),p)') . (128) 



" If there is a local contribution to the proper self-energy, one separates in complete analogy to the scalar 
equation (118). The decomposition (126) is taken for the non-local part of the self-energy, while the local 
contribution gives rise to an effective space-time dependent fermion mass term. 
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To see this note that = comes first on the imaginary path ^ = [0, — z/3], while 
— —iji comes latest such that sign<^(jc^ — y^) contributes opposite signs on the left 
and on the right of the equation. Eq. (128) can be rewritten in a more standard form as a 
fluctuation-dissipation relation for bosons:^ 



(129) 



with nBE(tt)) = {e^^ — 1) ^ denoting the Bose-Einstein distribution function. Eq. (129) 
relates the spectral function to the statistical propagator. While p'^^°^ encodes the in- 
formation about the spectrum of the theory, one observes from (129) that the function 
F^'^i) encodes the statistical aspects in terms of the particle distribution function hbe- 
In the same way one obtains for the Fourier transforms of the spectral and statistical 
components of the self-energy the thermal equilibrium relation 

E5f'i)((«,p) = -/Q+nBE(«)) 4''^(«,P). (130) 

We note that the ratio Ep^''''(£0,p)/2(U plays in the limit of a vanishing co-dependence 
the role of the decay rate for one-particle excited states with momentum p. 

For fermions the anti-periodicity condition of the fermionic propagator in ther- 
mal equilibrium, A(jc,y)|^.o=o ~ ~^{^^y)\}0=^ifi^ implies a corresponding fluctuation- 
dissipation relation. The difference is that the Bose-Einstein distribution in (129) is 
replaced by the Fermi-Dirac distribution /ifdIw) = (e^® + 1)^^ according to 1/2 + 
n'SE{(0) 1/2 — nFD(ft)) in the respective relation. 

It is important to realize that out of equilibrium F and p are not related by the 
fluctuation-dissipation relation! In contrast to the nonequilibrium theory, the relation 
(129) is a manifestation of the tremendous simplification that happens if the system is 
in thermal equilibrium. An even more stringent reduction occurs for the vacuum where 
nB^{o}) = 0. In this respect, nonequilibrium quantum field theory is more complicated 
since it admits the description of more general situations. Of course, the nonequilibrium 
theory encompasses the thermal equilibrium or vacuum theory as special cases. We leave 
now this equilibrium detour and return to the nonequilibrium case. 



3.4.3. Nonequilibrium evolution equations 

Out of equilibrium we have to follow the time-evolution both for the statistical propa- 
gator, F, as well as for the spectral function, p. The evolution equations are obtained 
from (113) with the help of the identities (116) and (120). Most importantly, once 
expressed in terms of F and p the time-ordering is explicit and the respective sign- 



In our conventions the Fourier transform of the real-valued antisymmetric function p{x,y) is purely 
imaginary. 
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functions appearing in the time-ordered propagator can be conveniently evaluated along 
the time contour . 

With the notation (118) the time evolution equation for the time-ordered propagator 
(113) reads 

[a + m2 {x\ G) ] G{x, y) + / ^E(x, z\G)G{z,y) = -ib{x-y), (131) 

where we have set 7? = 0. The influence of the initial-time sources encoded in R is 
discussed below. For the evaluation along the time contour involved in the integration 
with = /^dz° /d^z we employ (116) and (120): 



E(x,z;G)G(z,3;) = i ^{yF {x,z)F{z.y) 



-^E^(x,z)p(z,3;)sign^(zO-/)-^EP(^,z)F(z,3;)sign^(xO-z«) (132) 

--YP (x, z)p (z, y) sign^ (x^ - z°) sign^ (z^ - 

The first term on the r.h.s. vanishes because of integration along the closed time contour 
^ (cf. Sec. 3.1). To proceed for the second term one splits the contour integral such that 
the sign-functions have a definite value, for instance 

dzOsign^(zO -/) = r dzO(-l) -I- / V = -2 r dzO (133) 

Jo Jy h 

for the closed contour with initial time t^ = 0. To evaluate the last term on the r.h.s. of 
Eq. (132) it is convenient to distinguish the cases 

(a) 0^(.yQ-/) = l : 

/ dzOsign^(x«-z°)sign^^(z«-/) = f dz«(-l)+ f dz«+ /V(-l), (134) 

(b) 0^(3;O-.yO) = l : 

/ dzOsign^(/-zO)sign^(zO-/) = f dz\-l)+ F dz°+ /V(-l). (135) 

One observes that (a) and (b) differ only by an overall sign factor ~ sign<^(;c^ —y^)- 
Combining the integrals therefore gives: 



I ji:{x,z;G)G{z,y) = j d'^z^f^ dz''Zp{x,z)F{z,y) - dz«E^(x,z)p(z, j) 

--sign^,(/-/) j^^ dz°Ep(x,z)p(z,3;) \ . (136) 
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One finally employs 



U,G{x,y) = U,F{x,y) - Uign.^{x'' -y'')U,p{x,y) - i5{x-y) (137) 

such that the 5-temi cancels with the respective one on the r.h.s. of the evolution equation 
(131). Here we have used 

-^-d^^[p{x,y)^ign^{x^-y^)\ = -^-^ign.^{x^ -y^)d},p{x,y) 

-i&^{x°-y^)d^p{x,y) 



-id{x-y) , 

where ( 1 17) is employed for the last line and to observe that a term ~ p (x, y) 5^ {x^ — y^) 
vanishes identically. Comparing coefficients, which here corresponds to separating real 
and imaginary parts, one finds from (136) and (137) the equations for F{x,y) and p{x,y). 
Using the abbreviated notation Jl^ dz = Jl^ dz^ J°°^ d^z we arrive at the coupled evolution 
equations for the spectral function and the statistical propagator: 



[D, + M\x)]p{x,y) = - dzi:p{x,z)piz,y), 

[n, + M^{x)\F{x,y) = - dzLp{x,z)F{z,y) + dzi:Fix,z)p{z,y). 



(138) 

These are causal equations with characteristic "memory " integrals, which integrate over 
the time history of the evolution. We emphasize that the presence of memory integrals 
is a property of the exact theory and in accordance with all symmetries, in particular 
time reflection symmetry. The equations themselves do not single out a direction of 
time and they should be clearly distinguished from phenomenological nonequilibrium 
equations, where irreversibility is typically put in by hand. Since these equations are 
exact they are fully equivalent to any kind of identity for the two-point functions such 
as Schwinger-Dyson/Kadanoff-Baym equations without further approximations. For 
0=0 the functional dependence of the self-energy corrections in (138) is given by 
= M^{F), Lf = 1:f{P,F) and Zp ='Lp{p,F). The case 7^ is discussed below. 
Note that the initial-time properties of the spectral function have to comply with 
the equal-time commutation relations (117). In contrast, for F{x,y) as well as its first 
derivatives the full initial conditions at ?o = need to be supplied in order to solve 
these equations. To make contact with the discussion of initial conditions in Sec. 3.2 
and Eq. (99), we consider for a moment the spatially homogeneous case for which 
F{x,y) = F{x^,y^;x — y) = /[d'^p/(2;r)'^] exp[«p(x — y)]F(jc*^, j°;p) and equivalently 
for p(x,y). In terms of the Fourier components F{t,t';'p) the solution of the integro- 
differential equations (138) requires the following initial conditions: 

^(^^';p)M'=o = ^p , I {dtF{t,t'\v>)+dj,F{t,t'\v>))\t=,,=Q = -^p^??, 
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dA,F{t,t'-v>)\t=t'=Q = 7)1 + ^. (139) 

Here we have used that the required correlators at initial time are identical to those given 
in Eq. (99) for the considered case 0=0, where we had suppressed the momentum 
labels in the notation. Accordingly, these are the very same parameters that have to be 
specified for the corresponding most general Gaussian initial density matrix (97). We 
emphasize that the initial conditions for the spectral function equation are completely 
fixed by the properties of the theory itself: the equal-time commutation relations (117) 

specify p(?,?';p)|f=r'=o = 0, <9/p(r,r';p)|,=,'=o = 1 and dtdtip{t,t'\i^)\t=t'=o = for the 
anti-symmetric spectral function. 

We are now in the position to discuss the role of the initial-time sources, which contain 
the information about the initial-time density matrix. According to the discussion of 
Sees. 3.2 and 3.3 the initial-time sources are fully described by the Fourier components 
of the bilinear source term R{t, t'; att = t' = tQ = for the case 0=0. Mathematically, 
the role of the initial-time sources for the evolution equations is rather simple: Since 
these sources have support only for the initial time tQ and vanish identically for times 
t 7^ tQ, they only fix the initial values for the correlators and their first derivatives.^ For 
simpler differential equations this property is well documented in the literature on the 
theory of Green's functions. In order to see that this indeed holds also for the case 
considered here, recall that the time evolution equations (138) are derived from (113) 
forR = 0. Eq. (139) shows that there is a one-to-one correspondence between the initial- 
time sources parametrizing the density matrix (97) and the required initial conditions for 
the solution of the time evolution equations. To check that no further dependencies on 
R remain in the evolution equations for times t ^ to, one notes that the 7?-dependence 
appears in (113) as a term ~ /.7?(x,z)G(z, j), which identically vanishes forx^ tQ. As 
a consequence, the equation governing G{x,y) (or F and p) cannot explicitly depend on 
R for its time arguments different than to, i.e. for all times of interest. 

The clear separation of the dynamical role of spectral and statistical components 
is a generic property of nonequilibrium field theory. As discussed in Sec. 3.4.2, the 
nonequilibrium theory encompasses standard vacuum theory as a special case where this 
separation is absent. This dichotomy for nonequilibrium time evolution equations is not 
specific to scalar field degrees of freedom. In terms of spectral and statistical components 
the equations for fermionic fields or gauge fields have a very similar structure as well. For 
instance, the respective form of the evolution equations for the fermion spectral function 
p^f\x,y) and statistical propagator F^-f\x,y) defined in (121) can be directly obtained 
from (138) by the l.h.s. replacement: 

[D, + M^]p{x,y) -[i^^-m^f^]p^f\x,y) (140) 

and equivalently for F^-f^ ix,y). On the r.h.s. of (138) then appear the respective fermion 

If) (f) 

propagators and self-energies Ep and Lp. as defined in (126). This can be directly 
verified from the equation of motion for the time-ordered fermion propagator (66) by 



We consider additional external sources to be absent. 
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convoluting it with A. For a free inverse propagator as in (71) for Dirac fermions this 
yields 



— rn^f) 



^{x,y) - i / T.^f\x,z)^{z.y) = i^{x-y) . 



(141) 



Following along the lines of the above discussion for scalars one finds for the fermion 
case the coupled evolution equations: 



■m 



if) 



dzY!-P{x,z)p^f\z.y)., 



F^f\x,y) 



dz4^\x,z)F(^Hz,y)- I dzL^P{x,z)p^f\z,y).(U2) 



if). 



Similarly, the nonequilibrium evolution equations for gauge fields can be obtained as 
well. For instance, denoting the full gauge field propagator by 



D^'ix^y) = F^\x,y) - '-p^J {x^y) sign^{x^ -/) 

for a theory with free inverse gauge propagator given by 

iDQl^{x,y) = [g^yn-{l-^-')d^dy]^5{x~y) 



(143) 



(144) 



for covariant gauges with gauge-fixing parameter ^ and vanishing "background" fields, 
one finds the respective equations from (138) by 

p, + M^]p{x,y) -[g\D-{l-^-')d^dy]^pJ^ix,y) (145) 

and equivalently for F^{x,y). Of course, the respective indices have to be attached 
to the corresponding self-energies on the r.h.s. of the equations. The derivation of the 
nonequilibrium gauge field evolution equations will be discussed in more detail in Sec. 6 
in the context of higher nPl effective actions with n> 2. 



3.4.4. Non-zero field expectation value 

In the presence of a non-zero field expectation value, 7^ 0, the form of the scalar 
evolution equations for the spectral and statistical function (138) remain the same. 
However, the functional dependence now includes = M^{^,F), Lf = Lf{^,p,F) 
and Zp = Zp{^,p,F). Note that the local self-energy correction (119) described by 
does not depend on the spectral function because the latter vanishes for equal-time 
arguments. For the A^-component scalar field theory (13) one has with (j)^ = ^a^a- 

Ml,{x-(^,F) = (^m^ + A_[F,,{x,x)+f{x)]^dab 
X 

+ —[Fab{x,x) + (l)a{x)(l)b{x)], (146) 
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with the respective field indices attached to (138). In this case the evolution equations 
for the spectral function and statistical two-point function (138) are supplemented by a 
differential equation for given by the stationarity condition (24), which yields the field 
evolution equation: 



8Tn 



5(f)a{x) 



(147) 



For comparison the corresponding equation for the classical field theory, 5S[^]/8^a{x), 
is obtained from (147) by the replacement Mab{x;^ =0,F) — > m^6ab and r2 = 0. For 
the above field evolution equation we have used that the one-loop type contribution to 
the 2PI effective action reads 



^TrGoH<^))G: 



1 



Sab+:^^a{x)^h{x) ]Fba{x,x) (148) 



for the classical inverse propagator (20), such that one can write 

m2,(x;<^^0,F) 



2 d[iTrG,'{^)G/2] 



m — 



5(j)a{x) 



(149) 



The solution of (147) requires specifying the field and its first derivative at initial time. 
In the context of the above discussion for homogeneous fields this just corresponds to 
specifying (98), which together with the required initial conditions (139) for the two- 
point function completes the correspondence with the initial-time density matrix (97). 



3.4.5. Lorentz decomposition for fermion dynamics 

The evolution equations for the fermion spectral function and statistical propagator 
(142) are rather complicated in general. For Dirac fermions each two-point function con- 
tains 16 complex components, which is often supplemented by additional field attributes 
such as "flavor" or "color". However, it is often not necessary in practice to consider all 
components due to the presence of symmetries, which require certain components to 
vanish identically without loss of generality. Depending on the symmetry properties of 
the initial state and interaction these terms remain zero under the nonequilibrium time 
evolution, which can dramatically simplify the analysis. It is typically very useful to 
decompose the fields p'^f\x,y) and F'^f\x,y) of Eq. (142) into terms that have definite 
transformation properties under Lorentz transformation. In order to ease the notation, in 

this section we will write p^^^ ^ p, F^-^^ i-^ F, "L^P h-> A and E^"* i— ^ C. Using a standard 
basis we write 

P=Ps + ijspp + JiiPv + JiiYsPa + 2^^ivpr ' (150) 

where o^y = §[7^1 7 7v] and = i'f'y^'f-'f' . The 16 (pseudo-)scalar, (pseudo-)vector and 
tensor components 

ps = trp, pp = -iixy^p, p^=ixfp, pl = ixy^fp, p^" = t~ra^V , (151) 
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are complex two-point functions. Here we have defined tr = ^tr where the trace acts in 
Dirac space. Equivalently, there are 16 complex components for F, A and C for given 
other field attributes. Using the hermiticity properties (125) they obey 

p[^\x.y) = -{pf{y.x)y , Fp{x,y)=(Ff\y,x))\ (152) 

where Y = {S,P,V,A,T}. Inserting the above decomposition into the evolution equa- 
tions (142) one obtains the respective equations for the various components displayed 
in (151). 

For a more detailed discussion, we first consider the l.h.s. of the evolution equations 
(142). In fact, the approximation of a vanishing r.h.s. corresponds to the standard mean- 
field or Hartree-type approaches frequently discussed in the literature. However, to 
discuss thermalization we have to go beyond such a "Gaussian" approximation: it is 
crucial to include direct scattering which is described by the nonvanishing contributions 
from the r.h.s. of the evolution equations. This is discussed in detail in Sec. 4. Starting 
with the l.h.s. of (142) for the evolution equation for the spectral function one finds: 



tr[(/^ 


' - 'nf)p] 


= {id^p!;^)-mfps , 




' - mf)p] 


= -i {id^p^) - mf pp , 


tr 


' - mf)p] 


= {id^ps)+i {idvPr^) -nifpl^ , 




'-mf)p] 


= i{id^pp) + \£^^y' {idyPTjs) -mfP^ , 


tr[a^^(/^ 


' - mf)p] 


= -i{id^p^-id'pl^)+e^'y' [idyp^^s] -nifPr 



(153) 



The corresponding expressions for the l.h.s. of the evolution equation for the statistical 
propagator F follow from (153) with the replacement p ^ F . We turn now to the inte- 
grand on the r.h.s. of (142) for the evolution equation of p. For the various components 
(151) one finds: 



fr[Ap] = AsPs-ApPp+A^Py^^-A^^Pa,^ 
1 
2 



+^Afpr,Mv, (154) 



-ztr [ysAp] = Aspp+Apps-iA'^pA.^ + iA^APv-ii 

+^e^^^^Ar,^vPr,y5, (155) 
fi-[/Ap] = Aspl^+A^ps-iApp^^ + iA^^pp + iAy^yPj^ 

+/A^>Kv + ^e^''^^ (AA,vPr,y5+Ar,vypA,5) , (156) 
tr[75/Ap] = Asp^+A^^ps-iApp!^+iA^pp + iAA,vp7 

+iA>'/ pAy + U^"^^ {Av,vPT,rS+AT,vrPv,5) , (157) 
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-i {A^yp'y -Alpli)+e^'y' {Ay,,p^^s -AajPv,s) 

+i (Aj pl-Alp^) + i {A^/pt/ -A^pr/) • (158) 

With the above expressions one obtains the evolution equations for the various Lorentz 
components in a straightforward way using (142). We note that the convolutions ap- 
pearing on the r.h.s. of the evolution equation (142) for F are of the same form than 
those computed above for p. The respective r.h.s. can be read off Eqs. (154)-(158) by 
replacing p ^ F for the first term and A ^ C for the second term under the integrals 
of Eq. (142) for F. We have now all the relevant building blocks to discuss the most 
general case of nonequilibrium fermionic fields. However, this is often not necessary in 
practice due to the presence of symmetries, which require certain components to vanish 
identically. 

In the following, we will exploit symmetries of the action (68) for the chiral quark- 
meson model described in Sec. 2.3.1 and of initial conditions in order to simplify the 
fermionic evolution equations: 

Spatial translation invariance and isotropy: We will consider spatially homogeneous 
and isotropic initial conditions. In this case it is convenient to work in Fourier space and 
we write 

p{x,y)^p{x^y^;x-y) = J ^ e'> p (/,/; p) , (159) 

and similarly for the other two-point functions. Moreover, isotropy implies a reduction 
of the number of independent two-point functions: e.g. the vector components of the 
spectral function can be written as 

pO(xO,/;p) =p«(x«,/;p), pv(/,/;p) = vpv(xO,/;p), (160) 

where p = |p| and v = p/p. 

Parity: The vector components Py (jc°, j^; p) and pvix^^y^; p) are unchanged under a 
parity transformation, whereas the corresponding axial- vector components get a minus 
sign. Therefore, parity together with rotational invariance imply that 

pO(.rO,/;p)=PA(x«,/;p)=0. (161) 

The same is true for the axial- vector components of F, A and C. Parity also implies the 
pseudo-scalar components of the various two-point functions to vanish. 

CP-invariance: For instance, under combined charge conjugation and parity transfor- 
mation the vector component of p transforms as 

p'y{x\y^p) p'y{y\x^p), 
Pv(xO,/;p) — . -pv(/,/;p), 

and similarly for and Ay. The F-components transform as 

F^{x\y';p) — -F«(/,/;;.), 
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and similarly for Cy and Cy. Combining this with the hermiticity relations (152), one 
obtains for these components that 



Rcp^{x^y-p) 

ReFO(xO,/;p) 
ReAOy{x^,y^;p) 
RcC^y{x^y-p) 



3pv(xO,/;;?)=0, 
3Fv/(xO,/;p)=0, 
3Av(xO,/;p)=0, 
3Cv(xO,/;;7)=0, 



(162) 



for all times and and all individual momentum modes. 

Chiral symmetry: The only components of the decomposition (150) allowed by chiral 
symmetry are those which anticommute with 75. In particular, chiral symmetry forbids 
a mass term for fermions {mf = 0) and we have 



P5(^°,/;p 



0, 



(163) 



and similarly for the corresponding components of F , A and C. We emphasize, however, 
that in the presence of spontaneous chiral symmetry breaking there is no such simplifi- 
cation. 

In conclusion, for the above symmetry properties we are left with only four indepen- 
dent propagators: the two spectral functions py{y}^ ,y^;-g) and py{x^ ,y^;-g) and the two 



corresponding statistical functions Fy{x ,y ■,t^) and ,y ,p). They are either purely 
real or imaginary and have definite symmetry properties under the exchange of their time 
arguments x^ ^ y^. These properties as well as the corresponding ones for the various 
components of the self-energy are summarized below: 

Py^ Ay'. 



py,Ay 

V ' 
Fy, Cy 



imagmary, symmetnc; 
real, antisymmetric; 
imaginary, antisymmetric; 
real, symmetric. 

The exact evolution equations for the spectral functions read (cf. Eq. (142)):^ 
■ :p^{x'y-p)=ppy{x',y^-p) 



dx^ 
+ 

. d 

+ 



dz'[AUx''.z'';p)p^{z',y';p)-Ayix',z';p)pviz^y'';p) 
py{x',y';p)=pp^{x',y';p) 



(164) 



dz^ 



A''yix'^,z';p)pyiz''y;p)-Ayix',z'';p)p\}iz''y;p) 



(165) 



We note that the following equations do not rely on the restrictions (162) imposed by CP-invariance: 
they have the very same form for the case that all two-point functions are complex. 
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Similarly, for the statistical two-point functions we obtain (cf. Eq. (142)): 



+ 



aO(/,zO;p)fO(zO,/;p)-Av.(/,zO;p)Fv.(z«,/.^^ 
'(^y{x^z'^-p)p^y{z^y^;p)-Cv{x\z'';p)pv{z\y'';p) 



(166) 



ij^Fy{x\y^-p)=pF^{x''J-p) 



+ 



dz^ 



dz^ 



aO(xO,z«;p)Fv/(z^/;p)-Ak(x°,zO;p)F«(z«,/;p) 
'c^y{x^z'^;p)pv{z\y';p)-Cy{x\z'';p)p^y{z^y^p) 



(167) 



The above equations are employed in Sec. 4 to calculate the nonequilibrium fermion 
dynamics in a chiral quark-meson model. 

The time evolution for the fermions is described by first-order (integro-)differential 
equations for F and p: Eqs. (164)-(167). As pointed out above, the initial fermion 
spectral function is completely determined by the equal-time anticommutation relation 
of fermionic field operators (cf. Eq. (124)). To uniquely specify the time evolution for F 
we have to set the initial conditions. The most general (Gaussian) initial conditions for F 
respecting spatial homogeneity, isotropy, parity, charge conjugation and chiral symmetry 
can be written as 



Fv{t,t' ,p)\,^t'=o 
^v(^^^p)l^=r'=0 



1 
2 
0. 



^-ni{p) 



(168) 
(169) 



Here nQ{p) denotes the initial particle number distribution, whose values can range 
between and 1 . 



3.5. Nonequilibrium dynamics from the 2PI loop expansion 

In Sec. 3.4.3 we have derived coupled evolution equations (138) for the statistical 
propagator F and the spectral function p. A systematic approximation to the exact 
equations can be obtained from the loop or coupling expansion of the 2PI effective 
action, as discussed in Sec. 2.1. This determines all the required self-energies using (25) 
and the decomposition identities (116) and (120) for scalars, and (121) for fermions. 
(Approximations for gauge field theories will be discussed in Sec. 6 in the context of 
nPI effective actions with n > 2.) We emphasize that all classifications of contributions 
are done for the effective action. Once an approximation order is specified on the level 
of the effective action, there are no further classifications on the level of the evolution 
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equations. This is a crucial aspect, which ensures the "conserving" properties of 2PI 
expansions (cf. the discussion in Sec. 1.1). The purpose of this section is to present the 
relevant formulae, whose physics will be explained in detail in the later sections. 

We consider first the case of the scalar 0{N) symmetric field theory with classical 
action (13) and a vanishing field expectation value. The case ^ 7^ is treated in Sec. 3.6. 1 
below. From the 2PI effective action to three-loop order (36) one finds with Fab{x,y) = 
F{x,y)5ab and pab{^-,y) = p{^-,y)^ab the two-loop self-energies: 



eW(x) 
^F{x,y) 



^~6Ar^^^'^)' 



F^{x,y)-^p^{x,y) 
F\x,y)-\p\x,y) 



(170) 
(171) 

(172) 



which enter (1 19) and (138). The two-loop contribution to the effective action adds only 
to E(°) and corresponds to a space-time dependent mass shift in the evolution equations. 
As is discussed in detail below in Sec. 4. 1 , the time evolution obtained from the two-loop 
2PI effective action (Hartree, or similarly leading-order large-N approximations) suffers 
from the presence of an infinite number of spurious conserved quantities, which are not 
present in the fully interacting theory. As an important consequence no thermalization 
can be observed to that order. A crucial ingredient for the description of nonequilibrium 
evolution comes from the three-loop contribution to the 2PI effective action as described 
by (171) and (172). 

For the chiral Yukawa model for Nf = 2 fermion flavors and A'^^ scalar fields with 
classical action (68) as described in Sec. 2.3.1, one obtains from the two-loop 2PI 
effective action the fermion self-energies entering Eqs. (164)-(167): 



(2^ 



<(/,/;q)p(x^/;p-q 



,0 ,,0. 



+p'f(xO,3;«;q)F(/,/;p-q) 



(173) 



C?(xO,/;p) = -h^ 



1 



d^^ 
(2^ 



-p^(xO,/;q)p(z",/;p-q) 



<(xO,/;q)F(xO,/;p-q) 



.0 „o., 



(174) 



Here py {x^ , y*^ ; q) and Fy {x^ , 3^*^ ; q) are the vector components of the fermion two-point 
functions as introduced in Sec. 3.4.5. The dynamics of the scalar two-point functions 
p{x^,y^;'p) and F{x^,y^;'p) is described by the evolution equations (138) with the scalar 
self-energies to this loop order: 



N} + 2 r d^g 
6Nj .1 {2k)' 



■F{x^,x^;q) 



(175) 
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Sp(^«,/;p) = -^/^p(f(x«,/;q)F,,,(/,/;p-q), (176) 

,2 . j3. 



EH^°,/;p) = -^/^[<(/,/;q)Fv,,(A/;p-q) 

(177) 



-ip(f(xO,/;q)pv,^(xO,/;p-q) 



3.6. Nonequilibrium dynamics from the 2PI l/N expansion 

We consider first the case of the scalar 0{N) symmetric field theory with classical 
action (13) and a vanishing field expectation value such that Fab{x,y) = F{x,y)5ab and 
Pab{x-,y) = p{x,y)5ab- The case 7^ is treated in Sec. 3.6.1. In the 1/A^-expansion of 
the 2PI effective action to next-to-leading order, as discussed in Sec. 2.4, the effective 
mass term M^{x; G) appearing in the evolution equations (138) is given by 

M^{x\F) = m^ + X^^^F{x,x). (178) 
6A' 

One observes that this local self-energy part receives LO and NLO contributions. In 
contrast, the non-local part of the self-energy (118) is nonvanishing only at NLO: 
E(x,y;G) = —X/{3N)G{x,y)I{x,y) and using the decomposition identities (116) and 
(120) one finds 

^F{x,y) = -—{F{x,y)lF{x,y)--p{x,y)Ip{x,y)y (179) 

M^^y) = -—[F{x,y)Ip{x,y)+p{x,y)lF{x,y)). (180) 
Here the summation function (84) reads in terms of its statistical and spectral compo- 



nents:'° 



~6 1 / ^^hi.^^z){F\z,y)--^p\z.y)) -2 j dzlF{x,z)F{z,y)p{z.y)V{\%l) 
Ip{x,y) = -F{x,y)p{x,y)-- j dzlp{x,z)F{z,y)p{z.y). (182) 



^'^ This follows from using the decomposition identity for the propagator (116) and I{x,y) = lFix,y) 
j/p(jc,y)signc^(x'' — y"). Cf. also the detailed discussion in Sec. 3.4.3. 
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using the abbreviated notation Jifdz = Jl~ dz^ J^^d'^z. Note that F{x,y), p{x,y), 
T,p{xjy), Ep(x,y), Ifix^y) and Ip{x,y) are real functions. 

For the Yukawa model of Sec. 2.3.1 the results from the i/Nf expansion of the 
2PI effective action to NLO can be directly inferred from (173)-(177). Recall that all 
classifications are done on the level of the effective action, as explained for fermions at 
the end of Sec. 2.4.1. The only difference between the NLO evolution equations and the 
result (173)-(177) from the two-loop effective action is that for the former the NNLO 
part ~ NV^ in (175) is dropped. 



3.6.1. Nonvanishing field expectation value 

We consider the scalar N -component field theory with classical action (13). In the 
presence of a nonzero field expectation value (^^ the most general propagator can no 
longer be evaluated for the diagonal configuration (35). Restoring all field indices the 
evolution equations (138) read 

[Uxdac + Ml^{x)\pcb{x,y) = - dzi:P,{x,z)pcb{z,y), 

[n,5ac+Mi{x)]F,b{x,y) = -J^ dzLP,{x,z)F,i,{z,y) 

+ dzLf,,{x,z)Pcb{z.y). (183) 
Jo 

The statistical propagator and spectral function components have the properties 
Kbi^^y) = Pab{x,y) = Fba{y,x) and p*b{x,y) = pab{x,y) = -pba{y,x). At NLO in the 
2PI l/N expansion M^^{x) = M^^(x;0,F) is given by (146) and the self-energies 
Kb 3^) = Kb {^^y'AiP, F) and (x, y) = {x,y;^,p, F) are obtained from (83) as 

A- 1 
Kbix.y) = ~^{^^(^'^) [<l>aix)(l)biy)+l^ab{x,y)] - -Ip{x,y)pabix,y) 

+PF{x,y)Fab{x,y) - ^Pp{x,y)Pabix,y)y (184) 

Kbi^^y) = -^{hi^^y)[^a{x)(l)b{y)+Fab{x,y)]+If{x,y)p^b{x,y) 

+Pp{x,y)Fab{x,y) +PF{x,y)Pab{x,y)y (185) 

The functions If (x, y) = 7f (x, p , F) and Ip (x, y) = Ip{x,y;p,F) satisfy the same equa- 
tions as for the case of a vanishing macroscopic field given by (181) and (182). The 
respective ^-dependent summation functions Pf{x,y) = Pf{x,y;^,p,F) and Pp{x,y) = 
Pp{x,y;^,p,F) are given by 

A f r° 

PF{x,y) = -—^HF{x,y)- dz [Hp{x,z)lF{z,y) +Ip{x,z)HF{z,y)] 
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+ / dz [Hf{x,z)Ip{z,y)+lF{x,z)Hp{z,y 

J 



dz / dvIp{x,z)Hp{z,v)Ip{v,y) 
Jo 

+ dz dvIp{x,z)Hp{z,v)lF{v,y) 
Jo Jo 



+ dz J^^ dvlF{x,z)Hp{z,v)Ip{v,y) 



(186) 



^p(^,j) = -3^<^p(-^,3')-y„ dz[//p(x,z)/p(z,3;)+/p(x,z)//p(z,3;)] 



with 



+ jl dz dv/p(x,z)//p(z,v)/p(v,3;) 



(187) 



Hf{x,y) = -^a{x)Fab{x,y)^b{y) , ^p(^,3') = -^a{x)pab{x,y)h{y) . (188) 

The time evolution equation for the field (147) for the 2PI effective action to NLO (83) 
is given by 



5ab+Mlj^{x-^ = Q,F)j^b{x) 

^ 3N Jo \-^P^^'y'>^^bix,y)+lFix,y)pab{x,y)](j)biy) 
= -J^ dyZP,{x,y;(l) = 0,F,p)(l)b{y). (189) 
The second equality follows from (185). 



3.7. Numerical implementation 

Beyond the leading-order approximation, the time evolution equations (138) and 
(147) are nonlinear integro-differential equations. The approximate self-energies ob- 
tained from a loop expansion or \/N expansion of the 2PI effective action are described 
in Sees. 3.5 and 3.6. Though these equations are in general too complicated to be solved 
analytically^^ without additional approximations, they can be very efficiently imple- 
mented and solved on a computer. Here it is important to note that all equations are 



Cf. Sec. 4.3 for an approximate analytical solution in the context of parametric resonance. 
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explicit in time, i.e. all quantities at some later time tf can be obtained by integration 
over the explicitly known functions for times t <tf for given initial conditions. This is 
a direct consequence of causality. In this respect, solving the nonequilibrium evolution 
equations can be technically simpler than solving the corresponding theory in thermal 
equilibrium. The reason is that for the study of thermal equilibrium the equation of the 
form (111) is typically employed, which has to be solved self-consistently since on its 
l.h.s. and r.h.s. the same variables appear. It involves the propagator for the full range of 
its arguments. In contrast, the time-evolution equations (138) for p(jc,y)|^o=fj ^0=^^ ^^^d 
F(x,y)|^o=,j '^ot depend on the r.h.s. on p{x,y) and F{x,y) for times > ti and 

y'^ >t2. To see this note that the integrands vanish identically for the upper time-limits 
of the memory integrals because of the anti-symmetry of the spectral components, with 
p (x, y) \-^=yO = and Ep (x, y) \ ^.o^y) = 0. As a consequence, only explicitly known quan- 
tities at earlier times determine the time evolution of the unknowns at later times. The 
numerical implementation therefore only involves sums over known functions. 

To be more explicit we consider first a scalar field theory. For spatially homogeneous 
fields it is sufficient to implement the equations for the Fourier components F{tj';'p) 
and p (f, t'; p) and we consider ^ = 0. The numerical implementation with ^ 7^ follows 
along the very same lines for the equation (147). Spatially inhomogeneous fields pose 
no complication in principle but are computationally more expensive. As an example 
we consider a time discretization t — nat, t' — mat with stepsize a, such that F{t,t') ^ 
F{n,m), and 



where we have suppressed the momentum labels in the notation. The second derivative is 
replaced by a finite-difference expression, which is symmetric in at <-> —at. It is obtained 
from employing subsequently the lattice "forward derivative" [F (n -|- 1 , m) — F (n, m)] / at 
and "backward derivative" [F{n,m) —F{n — l,m)]/af. The integral is approximated 
using the trapezoidal rule with the average function value [F{n,m) +F{n + l,m)]/2 in 
an interval of length at. The above simple discretization leads already to stable numerics 
for small enough stepsize at, but the convergence properties may be easily improved 
with more sophisticated standard estimators if required. 

As for the continuum the propagators obey the symmetry properties F(n,m) = 
F{m,n) and p{n,m) = —p(m,n). Consequently, only "half" of the (n,m)-matrices have 
to be computed and p{n,n) = 0. Similarly, since the self-energy Ep and the summation 
function Ip appearing in the l/N expansion to NLO (cf. Eq. (182)) are antisymmet- 
ric in time one can exploit that Zp{n,n) and Ip{n,n) vanish identically. For the case of 
the NLO time evolution the equations (138) with (179)-(182) are used to advance the 
matrices F(n,m) and p(n,m) stepwise in the "n-direction" for each given m. As ini- 
tial conditions one has to specify F(0,0;p), F(l,0;p) and F(l, l;p), while p(0,0;p), 
p(l,0;p) and p(l, l;p) are fixed by the equal-time commutation relations (117). The 
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(190) 








(191) 



time discretized versions of (138) read: 



+ l,m;p) = 2F(n,m;p) —F{n— l,m;p) 

2 f 2 2-1 N + 2 f 1 
-Uf <p +m + A J^F{n,n;k) |>F(n,m;p) 

-a3|Ep(n,0;p)F(0,m;p)/2-Ef(n,0;p)p(0,m;p)/2 (192) 

m— 1 , 

+ (Ep(n,/;p)F(/,m;p) -Ef (n,/;p)p(Z,m;p) 
1=1 ^ 

n-l 

(n > m) and similarly for p. These equations are explicit in time: Starting with n= I, 
for the time step n + l one computes successively all entries with m = 0, . . . , n + 1 from 
known functions at earlier times. At first sight this property is less obvious for the non- 
derivative expressions (181) for If and (182) for Ip whose form is reminiscent of a gap 
equation. However, the discretized equation for Ip, 

/p(n,m;q) = ^ J^I^F{n,m;q-k)p{n,m;k) 

n-i 1 
-at Ip{nJ;q)F{l,m;q-k)p{l,m;k)y (193) 

I=m+l J 

shows that all expressions for Ip{n,m) are explicit as well: Starting with m = n where 
Ip vanishes one should lower m = n, . . . ,0 successively. For m = n—\ one obtains an 
explicit expression in terms of F{n,m) and p{n,m) known from the previous time step 
in n. For m = n — 2 the r.h.s. then depends on the known function /p (n, n — 1 ) and so on. 
Similarly, for given Ip{n,m) it is easy to convince oneself that the discretized Eq. (182) 
specifies If{n,m) completely in terms of 7f (n,0), . . . ,7^ (n,m — 1), which constitutes an 
explicit set of equations by increasing m successively from zero to n. 

It is crucial for an efficient numerical implementation that each step forward in time 
does not involve the solution of a self consistent or gap equation. This is manifest in 
the above discretization. The main numerical limitation of the approach is set by the 
time integrals ("memory integrals") which grow with time and therefore slow down the 
numerical evaluation. Typically, the influence of early times on the late time behavior 
is suppressed and can be neglected numerically in a controlled way. In this case, it is 
often sufficient to only take into account the contributions from the memory integrals 



For the discretization of the time integrals it is useful to distinguish the cases n > m and n < m. We 
compute the entries F(« + 1 , m) from the discretized equations for n>m except for n + l = m where we 
have to use the equations for « < m. 
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for times much larger then the characteristic inverse damping rate (cf. Sec. 4.1.4 below). 
An error estimate then involves a series of runs with increasing memory time. 

For scalars one can use a standard lattice discretization for a spatial volume with 
periodic boundary conditions. For a spatial volume V = (Nsa)'^ with lattice spacing a 
one finds for the momenta: 




(194) 



where = 0, . . . , A^^ — 1 . This can be easily understood from acting with the correspond- 
ing finite-difference expression (190) for space-components: d^e~^P^ ^ ^-ipx^^ipa _|_ 
g-!p« _2]/a^ = — e^'P-'' 4 sin^ {pa/ 2) /a^ . On the lattice there is only a subgroup of the ro- 
tation symmetry generated by the permutations of Px,Py,Pz the reflections pj^ ^ —p^ 
etc. for d = 3. Exploiting these lattice symmetries reduces the number of independent 
lattice sites to {Ns + 1 ) {Ns + 3) (A^^ + 5)/48. We emphasize that the self-energies are cal- 
culated in coordinate space, where they are given by products of coordinate-space corre- 
lation functions, and then transformed back to momentum space. The coordinate-space 
correlation functions are available hy fast Fourier transformation routines. 

The lattice introduces a momentum cutoff n/a, however, the renormalized quantities 
are insensitive to cutoff variations for sufficiently large n/a. We emphasize here again 
that it is often convenient to carry out the numerical calculations using unrenormalized 
equations, or equations where only the dominant (quadratically) divergent contributions 
in the presence of scalars are subtracted. Cf. the discussion at the end of Sec. 2.2. In 
order to study the infinite volume limit one has to remove finite size effects. Here this is 
done by increasing the volume until convergence of the results is observed. We finally 
note that the nonequilibrium equations are very suitable for parallel computing on PC 
clusters using the MPI standard. 

For fermionic field theories the numerical implementation is more involved than for 
purely scalar theories. Consider the evolution equations for the Yukawa model (164)- 
(167), together with the self-energies (173)-(174). The structure of the fermionic equa- 
tions is reminiscent of the form of classical canonical equations. In this analogy, Fy {tj') 
plays the role of the canonical coordinate and Fy{t,t') is analogous to the canonical mo- 
mentum. This suggests to discretize Fy{t,t') and py{t,t') at t — t' = Inat (even) and 



" For time evolution problems the volume which is necessary to reach the infinite volume limit to a 
given accuracy can depend on the time scale. This is, in particular, due to the fact that finite systems 
can show characteristic recurrence times after which an initial effective damping of oscillations can be 
reversed. The observed damping can be viewed as the result of a superposition of oscillatory functions 
with differing phases or frequencies. The recurrence time is given by the time after which the phase 
information contained in the initial oscillations is recovered. Then the damping starts again until twice 
the recurrence time is reached and so on. In the LO approximation one can explicitly verify that the 
observed recurrence times e.g. for the correlation Fix^x) scales with the volume or the number of lattice 
sites to "infinity". We emphasize that the phenomenon of complete recurrences, repeating the full initial 
oscillation pattern after some characteristic time, is not observed once scattering is taken into account. 
Periodic recurrences can occur with smaller amplitudes as time proceeds and are effectively suppressed 
in the large-time limit. 



Introduction to Nonequilibrium Quantum Field Theory 



54 



Fy{t,t') and py{t,t') at t — t' = {In + (odd) time-like lattice sites with spacing at. 
This is a generalization of the "leap-frog" prescription for temporally inhomogeneous 
two-point functions. This implies in particular that the discretization in the time direc- 
tion is coarser for the fermionic two-point functions than for the bosonic ones. This 
"leap-frog" prescription may be easily extended to the memory integrals on the r.h.s. of 
Eqs. (164)-(167) as well. 

We emphasize that the discretization does not suffer from the problem of so-called 
fermion doublers. The spatial doublers do not appear since (164)-(167) are effectively 
second order in x-space. Writing the equations for Fy{t^t' ,1) and pv{tj',x) starting 
from (164)-(167) one realizes that instead of first order spatial derivatives there is a 
Laplacian appearing. Hence we have the same Brillouin zone for the fermions and 
scalars. Moreover, time-like doublers are easily avoided by using a sufficiently small 
stepsize in time at /as- 

The fact that Eqs. (164)-(167) with (173)-(174) and the respective scalar ones contain 
memory integrals makes numerical implementations expensive. Within a given numeri- 
cal precision it is typically not necessary to keep all the past of the two-point functions 
in the memory. A single PHI desktop workstation with 2GB memory allows us to use a 
memory array with 470 time-steps (with 2 temporal dimensions: t and t'). For instance, 
we have checked for the presented runs in Sec. 4.2 that a 30% change in the memory 
interval length did not alter the results. For a typical run 1-2 CPU-days were neces- 
sary. Much shorter times can be achieved with parallel computing on PC clusters using 
the MPI standard. The shown plots in that section are calculated on a 470 x 470 x 32-^ 
lattice. (The dimensions refer to the t and t' memory arrays and the momentum-space 
discretization, respectively.) By exploiting the spatial symmetries described in Sec. 3.4.5 
the memory need could be reduced by a factor of 30. We have checked that the infrared 
cutoff is well below any other mass scales and that the UV cutoff is greater than the 
mass scales at least by a factor of three. In the evolution equations we analytically sub- 
tract only the respective quadratically divergent terms. To extract physical quantities we 
follow the time evolution of the system for a given lattice cutoff and measure e.g. the 
renormalized scalar mass from the oscillation frequency of the correlator zero-modes to 
set the scale. 



4. NONEQUILIBRIUM PHENOMENA 

Our aim in this section is to study quantum field theories which capture important as- 
pects of nonequilibrium physics and which are simple enough that one can perform a 
precise quantitative treatment. We employ first the two-particle irreducible \/N expan- 
sion for the scalar A^-component field theory (13). Considering subsequent orders in the 
expansion corresponds to include more and more aspects of the dynamics. This allows 
one to study the influence of characteristic ingredients such as scattering, off-shell and 
memory effects on the time evolution. We note that time-reflection invariance for the 
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equations and conservation of energy is preserved at each order in the 1 /A^-expansion.^"^ 
It should be stressed that during the nonequilibrium time evolution there is no loss of 
information in any strict sense. Here we consider closed systems without coupling to an 
external heat bath or external fields. There is no course graining or averaging involved in 
the dynamics. The important process of thermalization is a nontrivial question in a cal- 
culation from first principles. Thermal equilibrium itself is time-translation invariant and 
cannot be reached from a nonequilibrium evolution on a fundamental level. It is strik- 
ing that we will observe below that scattering drives the evolution very closely towards 
thermal equilibrium results without ever deviating from them for accessible times. 

Here we study the time evolution for various initial condition scenarios away from 
equilibrium. One corresponds to a "quench" often employed to mimic the situation of 
a rapidly expanding hot initial state which cools on time scales much smaller than the 
relaxation time of the fields. Initially at high temperature we consider the relaxation 
processes following an instant "cooling" described by a sudden drop in the effective 
mass. Another scenario is characterized by initially densely populated modes in a narrow 
momentum range around ±pts. This so-called "tsunami" initial condition is reminiscent 
of colliding wave packets moving with opposite and equal momentum. A similar non- 
thermal and radially symmetric distribution of highly populated modes may also be 
encountered in the so-called "color glass" state at saturated gluon density. Of course, 
a sudden change in the two-point function of a previously equilibrated system or a 
peaked initial particle number distribution are general enough to exhibit characteristic 
properties of nonequilibrium dynamics for a large variety of physical situations. We 
will first consider results for the 1 + 1 dimensional quantum field theory, since technical 
aspects such as the infinite volume limit, renormalization and large times can be all 
implemented with great rigour. The study also has the advantage that the physics of off- 
shell effects can be very clearly discussed since they play a more important role for one 
spatial dimension as compared to three dimensions. The conclusions we will draw are 
not specific to scalar theories in low dimensions, and we will compare below with the 
corresponding results in 3 -|- 1 dimensions for scalars interacting with fermions. 

For the 2PI l/N expansion we will first consider the time evolution at leading order 
(LO) and compare it with the next-to-leading order dynamics (NLO). The "mean-field 
type" LO approximation has a long history in the literature. However, a drastic con- 
sequence of this approximation is the appearance of an infinite number of conserved 
quantities, which are not present in the interacting finite -A'^ theory. These additional con- 
stants of motion can have a substantial impact on the time evolution, since they strongly 
constrain the allowed dynamics. As an important consequence of the LO approximation, 
the late-time behavior depends explicitly on the details of the initial conditions and the 
approach to thermal equilibrium cannot be observed in this case as is shown below. The 
extensive use of this approximation was based on the hope that deviations from the LO 
behavior are not too sizeable for not too late times. By taking into account the NLO 
contributions it turns out, however, that for many important questions there are substan- 
tial corrections long before thermalization sets in. In contrast, we will find that various 
approximations of the 2PI effective action which go beyond "mean-field type" (LO, 



These properties are also taken care of by the employed numerical techniques. 
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Hartree or any Gaussian) dynamics show thermalization and give comparable answers. 
The 2PI approximations beyond mean-field share the property that the spurious con- 
stants of motion are absent, which emphasizes the important role of conservation laws 
for the nonequilibrium dynamics: fake conserved quantities keep the information about 
the initial conditions for all times and can spoil any effective loss of details about the 
early-time behavior. 



4.1. Scattering, off-shell and memory effects 



4.1.1. LO fixed points 

For simplicity we consider spatially homogeneous field expectation values ^a{t) — 
(^»a(?,x)), such that we can use the Fourier modes Fab{t,t';'p) and pa^(?,r';p).'^ The so- 
lution of (138) in the limit of a free field theory (A = 0) describes modes which oscillate 
with frequency ^/p^ + rrP-/27t for unequal-time functions F(r,0;p) and p(?,0;p), un- 
less they are not identically zero. Equal-time correlation modes F{tj;'p) oscillate either 
with twice that frequency or they are constant in time (p(f,f;p) = according to Eq. 
(117)). The latter corresponds to solutions which are translationally invariant in time, 
le.F{t,t';p)=F{t-t';-p) and p{t,t';p) =p{t-t';-p). 

The LO contribution to the 2PI effective action (78) adds a time-dependent mass shift 
to the free field evolution equation. The resulting effective mass term, given by (146) 
for N ^ oo/is the same for all Fourier modes and consequently each mode propagates 
"coUisionlessly". There are no further corrections since according to (79) and (83) the 
self-energies E/r and Ep are ff{l/N) and vanish in this limit! The evolution equations 
(138) for this approximation read: 



d^ + -^d^{t)+M\t;(l) = 0,F 



(195) 



with 



6A^' 



M^{t;^,F) 







6A^ 



(196) 



where /p = / d'^p/{2nY. In this case one observes that the evolution of F and is 
completely decoupled from p. Similar to the free field theory limit, to LO the spectral 
function does not influence the time evolution of the statistical propagator. The reason is 
that in this approximation the spectrum consists only of "quasiparticle" modes of energy 
Q}p{t) = -\/p2 +M^{t) with an infinite life-time. The associated mode particle numbers 



^ Note that a spatially homogeneous field expectation value contains all fluctuations arising from an 
inhomogeneous (^(x) as well. From a practical point of view the extraction of physics related to inhomo- 
geneous field fluctuations is, of course, more difficult since it is contained in higher correlation functions. 
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are conserved for each momentum separately. In contrast, in the interacting quantum 
field theory particles decay, get created, annihilated and the notion of a conserved 
particle number is clearly absent for real scalar fields. Note that there are also no memory 
integrals appearing on the r.h.s. of (195), which incorporate in particular direct scattering 
processes. 

As an example we consider a "tsunami" type initial condition in 1 + 1 dimension. The 
initial statistical propagator is given by 

F{0,0;p)drd,>F{tj';p)\,^r'=o = [^^(0) + 1/2]^ , (197) 

(^)(0) = <9,<^) (01^=0 = 0. (198) 

Here Fab{t,t'\ p) = F{tj'; p)5ab, which is valid for all times with these initial conditions. 
The mass term M^(0) is given by the gap equation (196) in the presence of the initial 
non-thermal particle number distribution 

np(0)=^e"^^'^'"^'^^'. (199) 

As the renormalization condition we choose the initial renormalized mass in vacuum, 
niR = M(0)|„(o)=o = 1, as our dimensionful scale. There is no corresponding coupling 
renormalization in 1 + 1 dimension. In these units the particle number is peaked around 
\p\ = Pts = 5mR with a width determined by a = 0.5mR and amplitude £/ = 10. We 
consider the effective coupling A / {6in^) = 1 . 

On the left of Fig. 3 we present the time evolution of the equal-time correlation 
modes F{tj;p) for different momenta: zero momentum, a momentum close to the 
maximally populated momentum pts and about twice pts- One observes that the equal- 
time correlations at LO are strictly constant in time. This behavior can be understood 
from the fact that for the employed "tsunami" initial condition the evolution starts 
at a time-translation invariant non-thermal solution of the LO equations. There is an 
infinite number of so-called fixed point solutions which are constant in time. If the real 
world would be well approximated by the LO dynamics then this would have dramatic 
consequences. The thermal equilibrium solution, which is obtained for an initial Bose- 
Einstein distribution instead of (199), is not at all particular in this case! Therefore, 
everything depends on the chosen initial condition details. This is in sharp contrast to the 
well-founded expectation that the late-time behavior may be well described by thermal 
equilibrium physics. Scattering effects included in the NLO approximation indeed drive 
the evolution away from the LO non-thermal fixed points towards thermal equilibrium, 
which is discussed in detail in Sec. 4.1.2 below. 

A remaining question is what happens at LO if the time evolution does not start 
from a LO fixed point, as was the case for the "tsunami" initial condition above. In 
the left graph of Fig. 4 we plot the evolution of M^{t) in the LO approximation as a 
function of time t, following a "quench" described by an instant drop in the effective 
mass term from 2M^(0) to M^(0). The initial particle number distribution is np{0) — 
l/(exp[v/p2 + 2M2(0)/ro] - 1) with Tq = 2M(0) and (^(0) = 0(0) = 0. The sudden 
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FIGURE 3. LEFT: Comparison of the LO and NLO time dependence of the equal-time correlation 
modes F{t,t\ p) for the "tsunami" initial condition (197)-(199). The importance of scattering included in 
the NLO approximation is apparent: the non-thermal LO fixed points become unstable and the "tsunami" 
decays, approaching thermal equilibrium at late times. RIGHT: Effective particle number distribution for 
a "tsunami" in the presence of a thermal background. The solid line shows the initial distribution which 
for low and for high momenta follows a Bose-Einstein distribution, i.e. Log[l + l/«p(0)] ~ e;,(0)/7b. 
At late times the non-thermal distribution equilibrates and approaches a straight line with inverse slope 

Tea > Tq. 



change in the effective mass term drives the system out of equilibrium and one can 
study its relaxation. We present M^{t) for three different couplings A = Aq = 0.5M^(0) 
(bottom), A = lOAo (middle) and A = 40 Ao (top). All quantities are given in units of 
appropriate powers of the initial-time mass M(0). Therefore, all curves in the left graph 
of Fig. 4 start at one. The time-dependent mass squared M^{t) shoots up in response to 
the "quench" and stays below the value 2M(0)^ of the initial thermal distribution. The 
amplitude of initial oscillations is quickly reduced and, averaged over the oscillation 
time-scale, the evolution is rapidly driven towards the asymptotic values. 

We compare the asymptotic values with the self-consistent solution of the LO mass 
equation (196) for constant mass squared M^^^ and given particle number distribution 

Mg\p ^m^ + ^ (np{0) + ^ . (200) 



gap 



The result from this gap equation is Mg^p = {1.01,1.10, 1 .29}M^(0) for the three values 
of A, respectively. For this wide range of couplings the values are in good numerical 
agreement with the corresponding dynamical large-time results inferred from Fig. 4 as 
{1.01, 1.11, 1.31}m2(0) at t = 200/M(0). We explicitly checked that at t = 400/M(0) 



Here the logarithmic divergence of the one-dimensional integral is absorbed into the bare mass param- 
eter using the same renormalization condition as for the dynamical evolution in the LO approximation, 
i.e. «2 + I / d£ („^,(o) + 1 ) (p2 +m2(0))- ^ m2(0). 
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FIGURE 4. LEFT: Shown is the time-dependent mass term M^{t) in the LO approximation for three 
different coupHngs following a "quench". All quantities are given in units of appropriate powers of the 
initial-time mass M(0). RIGHT: Time dependence of the equal-time zero-mode F{t,t;p = 0) after a 
"quench" (see text for details). The inset shows the mass term M^{t), which includes a sum over all 
modes. The dotted Unes represent the Hartree approximation (LO+), while the solid lines give the NLO 
results. The coupling is X/6N = 0. nM-{0) for = 4. 



these values are still the same. One concludes that the asymptotic behavior at LO is well 
described in terms of the initial particle number distribution «p(0). The latter is not a 
thermal distribution for the late-time mass terms with values smaller than 2M^{0). 

The LO approximation for p) and p{t,t'; p) becomes exact'^ for t,t' 0, since 

the memory integrals on the r.h.s. of the exact Eqs. (138) vanish at initial time. For very 
early times one therefore expects the LO approximation to yield a quantitatively valid 
description. However, from Fig. 4 one observes that the time evolution is dominated 
already at early times by the non-thermal LO fixed points. As discussed above, the 
latter are artefacts of the approximation. Though the precise numerical values of the 
LO fixed points are pure artefacts, we emphasize that the presence of approximate fixed 
points governing the early-time behavior is a qualitative feature that can be observed 
also beyond LO (cf. below). The question of how strongly the LO late-time results 
deviate from thermal equilibrium depends of course crucially on the details of the initial 
conditions. Typically, time- and/or momentum- averaged quantities are better determined 
by the LO approximation than quantities characterizing a specific momentum mode. 
This is exemplified on the right of Fig. 4, which shows the equal-time zero-mode 
F{tj;p = 0) along with M^{t) including the sum over all modes. Here we employ 
a "quench" with a larger drop in the effective mass term from 2.9M^(0) to M^(0). 
The initial particle number distribution is ^^(0) = 1/ {exp[^/ p^ +M^{0) /Tq] — 1) with 
Tq = 8.5M(0). In the figure the dotted curves show the dynamics obtained from an 
"improved" LO (Hartree) approximation, L0+, that takes into account the local part of 
the NLO self-energy contribution. The resulting equations have the very same structure 



This is due to the fact that we choose to consider Gaussian initial conditions. 
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FIGURE 5. LEFT: Shown is the evolution of the unequal-time correlation F{t,Q;p~Q) after a 
"quench". Unequal-time correlation functions approach zero in the NLO approximation and correlations 
with early times are effectively suppressed (X/6N ^ (5/6A^~ 0.083) M^(0) for = 10). In contrast, there 
is no decay of correlations with earlier times for the LO approximation. RIGHT: The logarithmic plot of 
\p{t,0;p = 0)1 and |F(f,0;/7 = 0)| as a function of time t shows an oscillation envelope which quickly 
approaches a straight line. At NLO the correlation modes therefore approach an exponentially damped 
behavior. (All in units of M(0).) 

as the LO ones, however, with the LO and NLO contribution to the mass term M^{t) 
included as given by Eq. (178) below. The large-time limit of the mass term in the L0+ 
approximation is determined by the L0+ fixed point solution in complete analogy to the 
above discussion. We also give in Fig. 4 the NLO results, which are discussed below. 

The effective loss of details about the initial conditions is a prerequisite for thermal- 
ization. At LO this is obstructed by an infinite number of spurious conserved quantities 
(mode particle numbers), which keep initial-time information. An important quantity 
in this context is the unequal-time two-point function F(?,0; p), which characterizes the 
correlations with the initial time. Clearly, if thermal equilibrium is approached then these 
correlations should be damped. On the left of Fig. 5 the dotted line shows the unequal- 
time zero-mode F(f,0; p = 0) following the same "quench" at LO as for Fig. 4 left. One 
observes no decay of correlations with earlier times for the LO approximation. Scatter- 
ing effects entering at NLO are crucial for a sufficient effective loss of memory about 
the initial conditions, which is discussed next. 



4.1.2. NLO thermalization 

In contrast to the LO approximation at NLO the self-energies 'Lf and Ep ~ ^{l/N) 
do not vanish. For the initial conditions (197) all correlators are diagonal in field index 
space and 0=0 for all times. In this case the evolution equations derived from the 
NLO 2PI effective action (79) are given by (138) with the self-energies (178)-(182). 
As discussed in Sec. 3.4.3, the NLO evolution equations are causal equations with 
characteristic "memory " integrals, which integrate over the time history of the evolution 
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taken to start at time to = without loss of generality. 

We consider first the same "tsunami" initial condition (197) as for the LO case 
discussed above. The result is shown on the left of Fig. 3 for A'^ = 10. One observes a 
dramatic effect of the NLO corrections! They quickly lead to a decay of the initially high 
population of modes around pts- On the other hand, low momentum modes get populated 
such that thermal equilibrium is approached at late times. In order to make this apparent 
one can plot the results in a different way. For this we note that according to (197) the 
statistical propagator corresponds to the ratio of the following particle number at initial 
time t = 0: 

2 

and the corresponding mode energy. Here we have defined: 



'^pW + ^T = [F{t,t;p)K{t,t;p)-Q\t,t;p)] ' , (201) 



K{t,t';p) = dtd,,F{t,t';p) , Q{t,t';p) = ^ [dtF{t,t';p) +dt,F{t,t';p)] . (202) 



£p(0 = (;)^) , (203) 



Since we employed 2(0, 0;p) = for the initial conditions (197) the initial mode energy 
is given by 

' K{t,t;p) ^ 

such that F{tj;p) = [np{t) + l/2]/£p{t). For illustration of the results we may use 
(201) and (203) for times t > in order to define an effective mode particle number 
and energy, where we set 2(?,?';p) = in (201) for the moment. The behavior of the 
effective particle number is illustrated on the right of Fig. 3, where we plot Log(l + 
l/^p(0) a function of ep(f). Note that for a Bose-Einstein distributed effective particle 
number this is proportional to the inverse temperature: Log(l + l/n^^{T)) ~ l/T. For 
the corresponding plot of Fig. 3 we have employed an initial "tsunami" at pts/m^ — 
2.5, where we added an initial "thermal background" distribution with temperature 
Tq /mR = A}^ Correspondingly, from the solid line {t = 0) in the right figure one observes 
the initial "thermal background" as a straight line distorted by the non-thermal "tsunami" 
peak. The curves represent snapshots at equidistant time steps AtmR = 10. After rapid 
changes in np{t) at early times the subsequent curves converge to a straight line to high 
accuracy, with inverse slope T^q/To = 1.175. The initial high occupation number in 
a small momentum range decays quickly with time. More and more low momentum 
modes get populated and the particle distribution approaches a thermal shape. 

The crucial importance of the NLO corrections for the nonequilibrium dynamics can 
also be observed for other initial conditions. The right graph of Fig. 4 shows the equal- 
time zero-mode F{tj; p = 0), along with M^{t) including the sum over all modes, fol- 
lowing a "quench" as described in Sec. 4.1.1. While the dynamics for vanishing self- 
energies Zf and Ep is quickly dominated by the spurious LO fixed points, this is no 
longer the case once the NLO self-energy corrections are included. In particular, one 
observes a very efficient damping of oscillations at NLO. This becomes even more pro- 
nounced for unequal-time correlators as shown for F{t,0;p = 0) in the left graph of 



18 



The initial mass term is M(0) /otr = 2.24 and A /6N = 0.5m| for = 4. 
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Fig. 5. We find that the unequal-time two-point functions approach zero for the NLO 
approximation and correlations with early times are effectively suppressed. Of course, 
time-reversal invariance implies that the oscillations can never be damped out to zero 
completely during the nonequilibrium time evolution, however, zero can be approached 
arbitrarily closely. In the NLO approximation we find that all modes F{t,t';p) and 
p t'; p) approach an approximately exponential damping behavior for both equal-time 
and unequal-time correlations. On the right of Fig. 5 the approach to an exponential be- 
havior is demonstrated for \p{t,0\p = 0)\ and \F{t,0\p = 0)\ with the same parameters 
as for the left figure. The logarithmic plot shows that after a non-exponential period at 
early times the envelope of oscillations can be well approximated by a straight line. From 

an asymptotic envelope fit of 0; /> = 0) to an exponential form ~ exp(— 7q'^^™^V) we 

obtain here a damping rate y^^"^"^^ = 0.016A/(0). For comparison, for the parameters 

employed for the right graph of Fig. 4 we find y^^™^^ = 0. 1 1M(0) from an exponential 
fit to the asymptotic behavior of F(f,0; p = 0)- The oscillation frequency of F(?,0; p — 0) 
is found to quickly stabilize around l.lM(0)/2;r ^ 0.18M(0), which is of the same or- 
der than the damping rate. Correspondingly, one observes in Fig. 4 an equal-time zero 
mode which is effectively damped out at NLO after a few oscillations. The oscillation 
frequency is always found to stabilize very quickly and to be quantitatively well de- 
scribed by ep(f)/2;r for the unequal-time modes (ep(?)/;r for the equal-time modes) 
with the effective mode energy ep(?) given by (203). We emphasize that the parameter 
M{t) given by (178) cannot be used in general to characterize well the oscillation fre- 
quencies of the correlation zero-modes. Beyond LO or Hartree the renormalized mass, 
which can be obtained from the zero-mode oscillation frequency, receives corrections 
from the non- vanishing self-energies Hf and Ep. 

The strong qualitative difference between LO and NLO appears because an infinite 
number of spurious conserved quantities is removed once scattering is taken into ac- 
count. It should be emphasized that the step going from LO to NLO is qualitatively very 
different than the one going from NLO to NNLO or further. In order to understand better 
what happens going from LO to NLO we consider again the effective particle number 
(201). It is straightforward by taking the time derivative on both sides of (201) to obtain 
an exact evolution equation for np{t) with the help of the exact relations (138): 



fdt"{ [Lp{t,t";p)F{t",t;p)-ZF{t,t";p)p{t",t;p)] <9,F(r,r';p) 

JtQ l- 

- [Lp {t,t"; p) dfF it", ? ; p) - , t"; p) <9,p r ; p)] F , ? ; p) 



(204) 



Here tQ denotes the initial time which was set to zero in (138) without loss of generality. 
Since E^ ~ ff{l/N) as well as Ep, one directly observes that at LO, i.e. for N oo^ 
the particle number for each momentum mode is strictly conserved: dtnp{t) = at LO. 
Stated differently, eq. (201) just specifies the infinite number of additional constants 
of motion which appear at LO. In contrast, once corrections beyond LO are taken 
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into account then (201) no longer represent conserved quantities. The first non-zero 
contribution to the self-energies occurs at NLO. As a consequence dtnp{t) ^ in general. 

For the LO approximation we have seen that there is an infinite number of time- 
translation invariant solutions (201) to the nonequilibrium evolution equations. Thermal 
equilibrium plays no particular role at LO. However, this is very different at NLO and 
beyond. It is very instructive to consider what is required to find solutions of (204) which 
are homogeneous in time. Of course, any time-translation invariant solution cannot 
be achieved if one respects all symmetries of the quantum field theory such as time 
reflection symmetry. However, we can consider the limit to — oo and ask under what 
conditions time-translation invariant F(?,?';p) =F('^°™)(? — ?';p) etc. represent solutions 
of (204). Using the (anti-) symmetry of the (spectral function) statistical propagator we 
have 



J 27Z 

i7(hom)^^_^/.p^ = [^cos[co{t-t')]F^^°'^\(0,p). (205) 
J 2n 



Note that in our conventions p^^°^\co;'p) = — p^'^°™)(— ft);p) is purely imaginary, while 
F(h°™)(G);p) -F(*^°"^)(-(0;p) is real. Since dtF^^'^'^^t -t';-p)\,=f = 0, time-translation 
invariant solutions of (204) require: 



lim fdt" [Ej?°'"Hr-?";p)<9,F(h°"^)(f"-?;p) 



-Z?°'")(r-r";p)a,p(h°-)(r"-^;p) 



0. (206) 



Inserting the Fourier transforms (205) and the respective ones for the self-energies and 
using 

lim ^^^^MzM = ;r5(a)), (207) 

{t-t())^^ CO 
one finds that the condition (206) can be equivalently written as 



dCQCQ 



Lf^"^^ (CO, p)F(h-) (CO, p) - (CO, p)p (''°-) (CO, ] 



= 0. (208) 



In contrast to the LO approximation, at NLO and beyond we find with ^ that 

time-translation invariant solutions are highly nontrivial. However, one observes from 
the equilibrium fluctuation-dissipation relation (129) with (130) that indeed thermal 
equilibrium singles out a solution which fulfills this condition. We emphasize that in 
equilibrium the fluctuation-dissipation relation is fulfilled at all orders in the 2PI 1 /N 
expansion. In this respect there is no qualitative change going from NLO to NNLO or 
beyond. In particular, very similar results as shown in Fig. 3 are obtained from the loop 
expansion of the 2PI effective action beyond two-loop order, which is demonstrated in 
Sec. 4.1.4. Three-loop is required since the two-loop order corresponds to a Gaussian 
approximation, which suffers from the same infinite number of spurious conserved 
quantities as the LO approximation. 
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We emphasize that the above discussion does not imply that the evolution always 
has to approach thermal equilibrium at late times. There are many initial conditions for 
which the system cannot thermalize. Obvious examples are initial conditions inferred 
from statistical ensembles which do not obey the clustering property. These do not 
exhibit thermal correlation functions in general. 



4.13. Detour: Boltzmann equation 

The nonequilibrium 2PI effective action can be employed to obtain ejfective kinetic 
or Boltzmann-type descriptions for quasiparticle number distributions. The latter are 
widely used in the literature and it is important to understand their merits and limits. 
The "derivation" of such an equation from the 2PI effective action requires a number of 
additional approximations which go beyond a loop or 1/A^ expansion of the 2PI effective 
action. At the end of this procedure we will find an irreversible equation, which therefore 
lost part of the symmetries of the underlying quantum field theory. The comparison 
will allow us, in particular, to clearly point out the role of neglected off-shell effects as 
compared to the quantum field theory description. 

Kinetic descriptions concentrate on the behavior of particle number distributions. 
The latter contain information about the statistical propagator F{x,y) at equal times 

= . Correlations between unequal times are outside the scope of a Boltzmann 
equation. We can therefore start the discussion from the exact equation for the effective 
particle number (204).^'^ For our current purposes it is irrelevant that we will consider 
a spatially homogeneous one-component field theory with vanishing field expectation 
value = 0). The steps leading to the Boltzmann equation are summarized as follows: 

1 . Consider a 2PI loop expansion ofT[G] to three-loop order This leads to self-energy 
corrections up to two loops, diagrammatically given by (^j , -(3 , which lead to 
the time-dependent mass term 

M\t)=m^ + ^ [ F{t,t;p) (209) 
2 Jp 

and 

E^(f,f';p) = -— / F(?,/;p-q 

D Jq,k 

Ep(f,f';p) = -— / p(f,/;p-q 

/ 7q,k 

(210) 



k) 



F{t,t';q)F{t,t';k) ~ -p{t,t';q)p{t,t';k) 
F{tj';q)F{tj';k)-^p{t,t';q)p{t,t';k) 



^ For this it is useful to rewrite the following characteristic term as: 

- / (ipF - p) = (Ef - /Ep/2) {F + ip/2) - (Ef + /Ep/2) {F - ip/2) . 

The difference of the two terms on the rh.s. can be directly interpreted as the difference of a "loss" and a 
"gain" term in a Boltzmann type description. 
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2. Choose a quasiparticle ansatz, i.e. a free-field type form for the two-point functions 
and their derivative. This corresponds to the replacements on the r.h.s. of (204): 



dtF{t",t-i^) 



dtpit",t;p) 



F{t",t;p) 



p{t",t;p 



(zip + 1/2) cos [cOp{t"-t)]/cOp, 
sin [(Op{t"-t)] /(Op, 

(np + l/2)sin[cOp(/'"-0] , 
-cos [(Op{t"-t)] . 



(211) 



(212) 



3. Here the particle number rip and mode energy (Dp may still depend weakly on time: 
One assumes a separation of scales, with a sufficiently slow time variation of np{t) such 
that one can pull all factors of «p(?) out of the time integral on the r.h.s. of (204). The 
particle numbers are then evaluated at the latest time of the memory integral. 

4. Send the initial time to to the remote past: to —o°. Of course, by construction 
the resulting equation is not meant to describe the detailed early-time behavior since 
to — ^ — °°. In the context of kinetic descriptions, one finally specifies the initial condition 
for the effective particle number distribution at some finite time and approximates the 
evolution by the equation with to in the remote past. 

A standard alternative derivation employs first steps 1. and 4., then a first-order 
gradient expansion of the two-point functions in the center coordinate {t + t')/2 and 
then a quasiparticle ansatz. However, for the typically employed first-order gradient 
approximation both approaches are fully equivalent. The current procedure has the 
advantage that one can send the initial time to —°° last, which allows one to discuss 
a few finite-time effects that are typically discarded. Applying the assumptions 1.-3. to 
the exact equation (204) leads after some lengthy shuffling of terms to the expression: 



3 Jsqk 2(t)p2(t)q2(t)k2((^ 




(I) 



+ 3[(1 +np)(l +nq)(l -|-nk)ns-npnqnk(l -l-ns)] 



(H) 



JtQ 

3 [(1 +np)(l +nq)nk«s-«p"q(l +"k)(l +«s)] 




+ 



(III) 




+ 



[{\ + np)nqn]ins - np{\ +nq)(l +nk)(l +«s)] 



(IV) 




(213) 
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The contributions on the r.h.s. of this equation have a typical "gain & loss" structure 
with a simple interpretation: 

(I) describes production and annihilation of four "quasiparticles" (0 ^ 4, 4 ^ 0). 
(II) and (IV) describe 1^3 and 3^1 processes. 

(Ill) describes 2^2 scattering processes, which are the well-known contributions to 
the standard Boltzmann equation. 

We emphasize that only the processes described by (III) lead to a non-zero contribution 
in the limit {t — to) oo. For time-independent mode energies (Op one can perform the 
time integrations in (213) explicitly. This leads using (207) to 5-functions, which enforce 
the strict conservation of the mode energies involved in the scattering. At this point the 
description of the quantum field theory has been fully reduced to the physics of colliding 
"on-shell" quasiparticles. Of course, strict energy and momentum conservation for these 
quasiparticles leads to vanishing contributions (I), (II) and (IV). The latter describe 
processes, which change the total particle number such as 1 — > 3 processes. We will 
see below that in a quantum field theory "off-shell" particle number changing processes 
such as described by (II) can play a very important role! An obvious example is given 
by the fact that the Boltzmann equation including the contribution (III) admits grand 
canonical thermal solutions with nonzero chemical potential /i : 

However, the real scalar quantum field theory is charge neutral and /i = 0! We conclude 
that taking into account only 2^2 scattering processes of "on-shell" quasiparticles ob- 
viously fails to describe the quantum world if for a given initial condition the chemical 
potential turns out to be large. We emphasize that total particle number changing pro- 
cesses are included in the original quantum field theory equation before applying the 
additional steps 2.-4. In particular, they do not only appear as contributions from finite 
initial-time effects as discussed here for illustration. The main "off-shell" effects are 
removed in step 2., when the quasiparticle assumption enters the description. 

The physics of "off-shell" processes can be very nicely observed in the 1 + 1 dimen- 
sional quantum field theory. Note that for one spatial dimension the contribution (III) 
vanishes as well in the limit {t — to) ^ oo\ In this case both momentum and energy con- 
servation, p + q-k-s = Q A a//?2+m2 + ^q^^M'^ - -JWTW- - y/s^^TW = 0, lead 
to an "ineffective" two-to-two scattering where the incoming and outcoming modes are 
unaffected: for instance with p = k and q = s the contribution (III) in Eq. (213) vanishes 
identically. The equation becomes trivial in this case. We conclude that for the 1 + 1 
dimensional quantum field theory there is no dynamics from the Boltzmann equation! 
In contrast, from the same three-loop approximation of the 2PI effective action without 
the steps 2.-4. we will see in the next section that the 1 + 1 dimensional quantum field 
theory thermalizes. 
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FIGURE 6. LEFT: Examples for the time dependence of the equal-time propagator F(f,f;p) with 
Fourier modes p = 0,3,5 from the 2PI three-loop effective action. The evolution is shown for three 
very different nonequilibrium initial conditions with the same energy density (all in initial-mass units). 
RIGHT: Wigner transform p {X^\ (0,p) of the spectral function as a function of O) at X'^ = 35 . 1 for p = Q 
(in units of ihr). Also shown are fits to a Breit-Wigner function (dotted) with (ooiTo) = (1.46,0.37). The 
inset shows a blow-up around 3(Bo. The expected bump from "off-shell" 1 <-> 3 processes is small but 
visible. 



4.1.4. Characteristic time scales 

Nonequilibrium dynamics requires the specification of an initial state. A crucial 
question of thermalization is how quickly the nonequilibrium system effectively looses 
the details about the initial conditions, and what are the characteristic stages of a partial 
loss of information. Thermal equilibrium keeps no memory about the time history except 
for the values of a few conserved charges. As a consequence, for the real scalar field 
theory thermalization requires that the late-time result is uniquely determined by energy 
density. 

In Fig. 6 we show the time dependence of the equal-time propagator F(?, t; p) for three 
Fourier modes p = 0, 3, 5 and three very different initial conditions with the same energy 
density. All quantities are given in appropriate units of the mass at initial time, and we 
consider 1 + 1 dimensions. For the solid line the initial conditions are close to a mean 
field thermal solution, the initial mode distribution for the dashed and the dashed-dotted 
lines deviate more and more substantially from a thermal equilibrium distribution. It is 
striking to observe that propagator modes with very different initial values but with the 
same momentum p approach the same large-time value. The asymptotic behavior of the 
two-point function modes are universal and uniquely determined by the initial energy 
density. 

One observes that after an effective damping of rapid oscillations the modes are 
still far from equilibrium. All correlation functions quickly approach an exponentially 

damped behavior. A characteristic rate y^'^'^"^^^ can be obtained from the zero mode 
of the unequal-time two-point function F{t,0;p = 0) with a corresponding time scale 
proportional to the inverse rate. In this early-time range correlations with the initial 
time are effectively suppressed and asymptotically F{t,Q;p = 0) — 0+. The early- 
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time range is followed by a smooth "drifting" of modes, which is characterized by 
a slow dependence of F{t,t'\p = 0) on {t + t')/2. The presence of such a regime is 
a prerequisite for descriptions based on gradient expansions in the center coordinate 

(t + 1') /2. Clearly, the early-time behavior for times ~ 1/ y^'^™^'' is beyond the scope of 
a gradient expansion. Though the exponential damping at early times is crucial for an 
effective loss of details of the initial conditions, it does not determine the time scale for 
thermalization. One typically finds very different rates for damping and for the late-time 
approach of t;p = 0) to thermal equilibrium. 

One of the "key" properties for the success of the quantum field theoretical description 
is the nontrivial (i.e. not the free-field or "5"-type) dynamical spectral function p . 
To analyze the spectral function (not to solve the dynamics!), we perform a Wigner 
transformation for the modes p {t, t'; p) and write with = (? + 1') /2, = t — t': 

f2X° „ 

ip(X^-(0,p)= ds^e"^' p(X^ + s^/2,X'^-s^/2;p). (215) 

J -2X0 

The / is introduced such that p(X'^;(0,p) is real. Because the spectral function is an- 
tisymmetric, p{X^;co,'p) = — p(X*^; — (U,p), and we will present the positive-frequency 
part only. Since we consider an initial- value problem with t,t' >0, the time integral over 
s^ = t-t' is bounded by ±2X^. 

On the right of Fig. 6 we display the Wigner transform for the zero momentum mode 
as a function of co. One clearly observes that the interacting theory has a continuous 
spectrum described by a peaked spectral function with a nonzero width. The peak is 
located at coo/mR = 1.46 in units of the initial renormalized mass (cf. Sec. 4.1.1), and 
the results are shown for mrX^ = 35.1 with X/mj^ = A?^ The inset shows a blow-up 
around SojQ/mj^ = 4.38. The expected bump in the spectral function is small but visible. 
We stress that this bump in the spectral function is kinematically forbidden for the "on- 
shell" approximation and arises from "off-shell" 1 3 processes. In Fig. 6 we also 
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Here we used a "tsunami" similar to the one discussed in Sec. 4.1.1, however, the shown results are 
insensitive to the initial condition details. 
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present fits to a Breit-Wigner spectral function 



Pbw(X«;CO,p) =[„2_^2(;^0)]2 + „2r2(;^0)' (216) 

with a width T^iX^) = 27^'^™^^ (X^). While the position of the peak can be fitted easily, 
the overall shape and width are only qualitatively captured. In particular, the slope of 
p (0, p) for small CO is quantitatively different. We also see that the Breit-Wigner fits 
give a narrower spectral function (smaller width) and therefore would predict a slower 
exponential relaxation in real time. 

The characteristic time scales observed for the nonequilibrium evolution of modes in 
Fig. 6 can be associated to 

1 . rapid oscillations of correlation functions with period ~ 1 / cOp described by the 
"peak" of the spectral function. 

2. damping of oscillations with inverse rate 1 / jp^'^^^^ described by a nonzero "width " 
Fp = 27p'^™^^ of the spectral function. In equilibrium the "width" is given by 

2(oF('^i^(£0,p) = -Z^p'^\(0,p) ~ ^(X^/N) . (217) 

3. late-time thermalization with inverse rate 1 / y^*'^™^ because of "off-shell" number 
changing processes. For the three-loop or NLO in l/N approximation of the 2PI 
effective action the processes changing the total particle number are perturbatively 
of order ~ A^/A^^ ("slow!"). This can be understood from the fact that the total 
particle number changing processes require a nonzero "width" ~ /N, and this 
width enters in jN) evolution equations for the two-point functions. 

One finds qualitatively the same characteristic ranges for the corresponding nonequilib- 
rium evolution in 3 + 1 dimensions, and we will consider results for the chiral "quark- 
meson" model below. An important quantitative difference between one and three spatial 
dimensions results from the fact that "on-shell" 2^2 processes contribute for the latter. 
However, as for the 1 -I- 1 dimensional case they do not change the total "quasiparticle" 
number. Total number changing "off-shell" or "on-shell" processes are required in gen- 
eral to reach thermal equilibrium. As we have seen above, "off-shell" number changing 
processes appear at NLO in the 2PI XjN expansion. In contrast, they can be achieved 
"on-shell" first at NNLO. E.g. for (/» "^-theory the lowest order "on-shell" (2 4) contri- 
bution appears from the five-loop diagram (cf. Sec. 2.1): 



Perturbatively, this contribution is of the same order ~ /N'^ than the "off-shell" 
particle number changing processes arising at three-loop or NLO in the l/N expansion 
of the 2PI effective action. In this respect it is interesting to observe that precision tests 
for the nonequilibrium dynamics as discussed in Sec. 5 indicate accurate results already 
at NLO in the 2PI expansion. 
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FIGURE 7. Fermion occupation number n'^' {t'-,p) for three different momentum modes as a function of 
time. The evolution is shown for two different initial conditions with same energy density. The long-time 
behavior is shown on a logarithmic scale for t > 30 m^^. 



4.2. Prethermalization 

Prethermalization is a universal far-from-equilibrium phenomenon which describes 
the very rapid establishment of an almost constant ratio of pressure over energy density 
(equation of state), as well as a kinetic temperature based on average kinetic energy. The 
phenomenon occurs on time scales dramatically shorter than the thermal equilibration 
time. As a consequence, prethermalized quantities approximately take on their final 
thermal values already at a time when the occupation numbers of individual momentum 
modes still show strong deviations from the late-time Bose-Einstein or Fermi-Dirac 
distribution. 

Here we consider the nonequilibrium evolution of quantum fields for a low-energy 
quark-meson model, which is described in Sec. 2.3.1. It takes into account two quark 
flavors with a Yukawa coupling ~ to a scalar a-field and a triplet of pseudoscalar pions, 
n. The theory corresponds to the well-known "linear c7-model", which incorporates the 
chiral symmetries of massless two-flavor QCD. The employed couplings in the action 
(68) with (69) are taken to be of order one, and if not stated otherwise h = X = I. We 
emphasize that the main results about prethermalization are independent of the detailed 
values of the couplings. Here we employ the 2PI effective action (70) to two-loop order 
given by (74). All quantities will be given in units of the scalar thermal mass 

Thermalization: In Fig. 7 we show the effective occupation number density of fermion 



For the relation to a nonperturbative expansion of the 2PI effective action to next-to-leading order in 
Nf see Sees. 2.3.1 and 2.4). 

The thermal mass m is evaluated in equilibrium. It is found to prethermalize very rapidly. The employed 
spatial momentum cutoff is A/m = 2.86. 
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FIGURE 8. Fermion and scalar mode temperatures Tl/'''\t) as a function of momentum p for various 
times. 



momentum modes, n^-^\t\ p), as a function of time for three different momenta . The 
plot shows two runs denoted as (A) and (B) with different initial conditions but same 
energy density. Run (A) exhibits a high initial particle number density in a narrow 
momentum range around ±p. This situation is reminiscent of two colliding wave packets 
with opposite and equal momentum. We emphasize, however, that we are considering 
a spatially homogeneous and isotropic ensemble with a vanishing net charge density. 
For run (B) an initial particle number density is employed which is closer to a thermal 
distribution. 

One observes that for a given momentum the mode numbers of run (A) and (B) 
approach each other at early times. The characteristic time scale for this approach is well 
described by the damping time ?damp(/') Irrespective of the initial distributions (A) or 

(B), we find (for pjmr^^^ 1) ^^{^p — 25 m^^ for fermions and ?^g|^p — 2%mr^ for scalars. 
In contrast to the initial rapid changes, one observes a rather slow or "quasistationary" 
subsequent evolution. The equilibration time ?eq — 95 itT^ is substantially larger than 
^damp and is approximately the same for fermions and scalars. Thermal equilibration is 
a collective phenomenon which is, in particular, rather independent of the momentum. 
As we have also observed for the 1 + 1 dimensional scalar theory in Sec. 4.1.4, mode 
quantities such as effective particle number distribution functions show a characteristic 
two-stage loss of initial conditions: after the damping time scale much of the details 
about the initial conditions are effectively lost. However, the system is still far from 
equilibrium and thermalization happens on a much larger time scale. 



" This quantity is directly related to the expectation value of the vector component of the field commu- 
tator {[\if, V/]) in Wigner coordinates and fulfills < n^^\t; p) <l- 
•^'^ The rate l/fdamp(p) is determined by the spectral component of the self-energy. 
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FIGURE 9. The ratio of pressure over energy density w as a function of time. The inset shows the early 
stages for two different coupHngs and demonstrates that the prethermahzation time is independent of the 
interaction details. 



We define mode temperatures (?) by equating tlie mode particle numbers 

(fs) 

fip ' (t) with a time and momentum dependent Bose-Einstein or Fermi-Dirac distribu- 
tion, respectively: 

np{t) = [oxp(cOp{t)/Tpit))±\]-' . (218) 

This definition is a quantum mechanical version of its classical counterpart as defined by 
the squared "generalized velocities". In thermal equilibrium with cOp ^ a/ and 
Tp = Teq equation (218) yields the familiar occupation numbers (jU — 0). Here the mode 

frequency (Op ''' {t) is determined by the peak of the spectral function for given time and 
momentum, as detailed for the scalar theory in Sec. 4.1.4. In Fig. 8 we show the fermion 
and scalar mode temperature as a function of momentum for various times t ^ ?damp- 
One observes that at late times, when thermal equilibrium is approached, all fermion 

and scalar mode temperatures become constant and agree: Tp^\t) — Tp\t) = T^q. In 
contrast, there are sizeable deviations from the thermal result even for times considerably 
larger than the characteristic damping time. 

Kinetic prethermalization: In contrast to the rather long thermalization time, prether- 
mahzation sets in extremely rapidly. In Fig. 9 we show the ratio of pressure over energy 
density, w = p/e, as a function of time. One observes that an almost time-independent 
equation of state builds up very early, even though the system is still far from equilib- 
rium! The prethermalization time fpt is here of the order of the characteristic inverse 
mass scale irT^ . This is a typical consequence of the loss of phase information by sum- 
ming over oscillating functions with a sufficiently dense frequency spectrum. In order 
to see that this phenomenon is not related to scattering or to the strength of the interac- 
tion, we compare with a smaller coupling in the inset and observe good agreement of 
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FIGURE 10. Chemical temperatures for scalars (upper curves) and fermions (lower curves) for different 
values of the coupling h and Teq. We also show the kinetic temperature Ti^„{t) (solid line), which 
prethermalizes on a very short time scale as compared to chemical equilibration. 



both curves. The dephasing phenomenon is unrelated to the scattering-driven process of 
thermalization. 

Given an equation of state, the question arises whether there exists a suitable definition 
of a global kinetic temperature Ti^ii. In contrast to a mode quantity such as Tp{t), a 
temperature measure which averages over all momentum modes may prethermalize. 
Building on the classical association of temperature with the mean kinetic energy per 
degree of freedom, we use here a definition based on the total kinetic energy £'kin(?): 

7idn(0=^ldn(0/Ceq. (219) 

Here the extensive dimensionless proportionality constant Ceq = £'kin,eq/7eq is given 
solely in terms of equilibrium quantities'^. Since total energy is conserved, the time scale 
when "equipartition" is reached (i.e. Ey^i^/E is approximately constant) also corresponds 
to a time-independent kinetic temperature. The latter equals the equilibrium temperature 
Teq if Eyiu/E has reached the thermal value. 

The solid line of Fig. 10 shows 7kin(0 normalized to the equilibrium temperature 
(for Teq/m =1). One observes that an almost time-independent kinetic temperature is 
established after the short-time scale fpt ~ /« ^ The time evolution of bulk quantities 
such as the ratio of pressure over energy density w, or the kinetic temperature T]^n, are 
dominated by a single short-time scale. These quantities approximately converge to the 
thermal equilibrium values already at early times and can be used for an efficient "quasi- 



^ For a relativistic plasma one has E^ia/N ~ e/n = aT. As alternatives, one may consider the weighted 
average T{t) = Y.n{t\p)T{t\p)/Y.n{t;p) where the sum is over aU modes, or a definition analogous to 
Eq. (220) below. 
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thermal" description in a far-from-equilibrium situation! 

Chemical equilibration: In thermal equilibrium the relative particle numbers of dif- 
ferent species are fixed in terms of temperature and particle masses. A system has chem- 
ically equilibrated if these ratios are reached, as observed for the hadron yields in heavy 
ion collisions. Obviously, the chemical equilibration time t^h will depend on details of 
the particle number changing interactions in a given model and tch < ?eq- In our model 
we can study the ratio between the numbers of fermions and scalars. For this purpose 

we introduce the chemical temperatures T^^'''\t) by equating the integrated number 

density of each species, n^f'^\t) = J d?p/ {2n)^ rip '''\t) , with the integrated Bose- 
Einstein/Fermi-Dirac form of distributions: 



The time evolution of the ratios T^^ ' {t) / Teq is shown in Fig. 10 for different values of 
the coupling constant h and the equilibrium temperature T^q. One observes that chemical 

equilibration with T^^ (t) = T^^^ (t) does not happen on the prethermalization time scale, 
in contrast to the behavior of 7kin(?). Being bulk quantities, the scalar and fermion 
chemical temperatures can approach each other rather quickly at first. Subsequently, a 
slow evolution towards equilibrium sets in. For the late-time chemical equilibration we 
find for our model fch — ^eq- However, the deviation from the thermal result can become 
relatively small already for times ? -C ?eq- 

Let us finally consider our findings in view of collisions of heavy nuclei and try to 
estimate the prethermalization time. Actually, tpt is rather independent of the details of 
the model like particle content, values of couplings etc. It mainly reflects a characteristic 
frequency of the initial oscillations. If the "temperature" (i.e. average kinetic energy 
per mode) sets the relevant scale one expects T tpi = const. (For low T the scale will 
be replaced by the mass.) For our model we indeed find r?pt ~ 2 — 2.5. We expect 
such a relation with a similar constant to hold for the quark-gluon state very soon after 
the collision. For T > 400 — 500 MeV we obtain a very short prethermalization time ?pt 
of somewhat less than 1 fm. This is consistent with observed very early hydrodynamic 
behavior in collision experiments, however, one further has to investigate the required 
isotropization of pressure. 



4.3. Far-from-equilibrium field dynamics: Parametric resonance 

In classical mechanics parametric resonance is the phenomenon of resonant amplifi- 
cation of the amplitude of an oscillator having a time-dependent periodic frequency. In 
the context of quantum field theory a similar phenomenon describes the amplification 
of quantum fluctuations, which can be interpreted as particle production. It provides an 



2^ We define fpt by | w(fpt) - Weq|/weq < 0.2 for t > fpt. 
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(220) 



Here g^-f^ = 8 counts the number of fermions and g^'''^ = 4 for the scalars. 



important building block for our understanding of the (pre)heating of the early universe 
at the end of an inflationary period, and may also be operative in various situations in 
the context of relativistic heavy-ion collision experiments. Here we will consider the 
phenomenon as a "paradigm" for far-from-equilibrium dynamics of macroscopic fields 
or one-point functions. Dynamics of correlation functions for vanishing — or similarly 
for constant — macroscopic fields has been described in detail above. There is a wealth 
of phenomena associated to non-constant fields such as parametric resonance, spinodal 
decomposition or in general the dynamics of phase transitions where the field can play 
the role of an order parameter. 

The example of parametric resonance is particularly challenging since it is a nonper- 
turbative phenomenon even in the presence of arbitrarily small couplings. Despite being 
a basic phenomenon that can occur in a large variety of quantum field theories, paramet- 
ric resonance is a rather complex process, which in the past defied most attempts for a 
complete analytic treatment even for simple theories. It is a far-from-equilibrium phe- 
nomenon that involves densities inversely proportional to the coupling. The nonpertur- 
batively large occupation numbers cannot be described by standard kinetic descriptions. 
So far, classical statistical field theory simulations on the lattice have been the only quan- 
titative approach available. These are expected to be valid for not too late times, before 
the approach to quantum thermal equilibrium sets in (cf. also Sec. 5). Studies in quan- 
tum field theory have been mainly limited to linear or mean-field type approximations 
(leading-order in large-A'', Hartree), which present a valid description for sufficiently 
early times. However, they are known to fail to describe thermalization and miss impor- 
tant rescattering effects (cf. also Sec. 4.1). The 2PI 1 /N expansion provides for the first 
time a quantitative nonperturbative approach in quantum field theory taking into account 
rescattering. 

Recall for a moment the classical mechanics example of resonant amplitude amplifi- 
cation for an oscillator with time-dependent periodic frequency. The amplitude y{t) is 
described by the differential equation y-\-ca^{t)y = with periodic (o{t + T) = (o{t) of 
period T. Since the equation is invariant under t ^ t + T there are periodic solutions 
y{t + T) = cy{t). This can be expressed as y{t) = d/^U.{t) with periodic U.{t + T) = 
n(?) . One concludes that for real c > 1 there is an instability with an exponential growth. 
For small elongations a physical realization of this situation is a pendulum with a peri- 
odically changing length as displayed: 

periodic 



In contrast to the mechanics example, in quantum field theory there will be no external 
periodic source. A large coherent field amplitude coupled to its own quantum fluctua- 
tions will trigger the phenomenon of parametric resonance. Mathematically, however, 
important aspects are very similar to the above classical example for sufficiently early 
times: The mechanical oscillator amplitude y plays the role of the statistical two-point 
function F in quantum field theory, and the periodic ft)^(?) plays the role of an effec- 
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tive mass term M^{^{t)) whose time dependence is induced by an oscillating macro- 
scopic field ^{t). Simple linear approximations to the problem, which have been much 
employed in the literature, are even mathematically equivalent to the above mechanics 
example. Accordingly, the well-known Lame-type solutions of the mechanics problem 
will also play a role in the quantum field theory study. Substantial deviations do, how- 
ever, quickly set in with important non-linear effects. 

4.3.1. Parametric resonance in the 0{N) model 

We consider the scalar A'^-component quantum field theory with classical action (13), 
and employ the 2PI 1 /A'^-expansion to next-to-leading order. The relevant equations of 
motion are given by (183)-(189) as described in Sec. 3.6.1. We will describe numerical 
solutions of these equations in 3 + 1 dimensions without further approximations. More- 
over, the approach allows us to identify the relevant contributions to the dynamics at 
various times and to obtain an approximate analytic solution of the nonlinear dynamics 
for the entire amplification range. It should be emphasized that the approach solves the 
problem of an analytic description of the dynamics at nonperturbatively large densities. 
It is necessary to take into account the infinite set of NLO 2PI diagrams as described 
in Sec. 2.4.1. As we will show in the following, each of these eventually contributes to 
the same order in the coupling A such that any finite order in loops or couplings is not 
sufficient. In this sense, the 2PI 1 /A/^-expansion to NLO represents a minimal approach 
for the controlled description of the phenomenon. This justifies the rather involved com- 
plexity of the approximation. 

We have in mind a situation reminiscent of that in the early universe after a period of 
chaotic inflation, driven by a macroscopic (inflaton) field. We consider a weakly coupled 
system that is initially in a pure quantum state, characterized by a large field amplitude 

l6N 

Ut) = o{t)Mo^j^5ai. (221) 

and small quantum fluctuations, corresponding to vanishing particle numbers at initial 
time. Here Mq sets our unit of mass, a(0) = ao is ~ ^(1) and dtO{t) \t=Q — 0. The initial 
statistical propagator contains a "longitudinal" component and (A^ — 1) "transverse" 
components F^: 

= diag{i^|,Fi,...,ia}, (222) 

with (9fF(?,0;p)|f=o = 0, and dtdjiF{t,t'\T^)\t=t'=Q = F^^{Q,Q\v) for ^ pure-state ini- 
tial density matrix (cf. Sec. 3.2). 

To get an overview, in Fig. 1 1 we show the contributions from the macroscopic field 
and the from the fluctuations to the (conserved) total energy as functions of time for 
A^ = 4 and A = 10^6. 27 The total energy Etot is initially dominated by the classical part 



Typical volumes (A^yflj)^ of A^^ = 36^8 with a,^ = 0.4-0.3 lead to results that are rather insensitive to 
finite-size and cutoff effects. 
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FIGURE 11. Total energy (solid line) and classical-field energy E^i (dotted line) as a function of time. 
The dashed line represents the fluctuation part iifiuc, showing a transition from a classical-field to a 
fluctuation dominated regime. 



Eci{t = 0), with reads in terms of the rescaled field (221): 

(223) 

t=t' 

where L3 denotes the spatial volume, and Efluc = ^^tot — Ed- More and more energy 
is converted into fluctuations as the system evolves. A characteristic time — denoted 
as fnonpert in Fig. 11 — is reached when both contributions become of the same size, 
i.e. Efluc — Ed- Before this time, the coherent oscillations of the field lead to a resonant 
enhancement of the statistical propagator Fourier modes 

F^{t,t';po) ^ e^^^'+''\ (224) 

in a narrow range of momenta around a specific value |p| ~ pq: this is parametric reso- 
nance. Apart from the resonant amplification in the linear regime, we identify two char- 
acteristic time scales — denoted as ^source and ^collect below — , which signal strongly en- 
hanced particle production in a broad momentum range. Nonlinear interactions between 
field modes cause the resonant amplification to spread to a broad range of momenta. 
More specifically, the initially amplified modes act as a source for other modes. The rate 
of the source-induced exponential amplification exceeds the characteristic rate 70 for the 
resonant amplification. This is illustrated in Figs. 12 and 13, where the effective parti- 
cle numbers are shown for various momenta as a function of time in the transverse and 
the longitudinal sector. Before the transition to a fluctuation dominated regime around 
^nonpert; One observes a rapid change of the particle numbers due to resonant as well as 
source-induced amplification. The source-induced amplification is crucial for the rapid 
approach to the subsequent, quasistationary regime, where direct scattering drive a very 
slow evolution towards thermal equilibrium. The transition to this slow regime around 
^nonpert Can be very well observed from Figs. 12 and 13. 

The fluctuation dominated regime is characterized by strong nonlinearities. For in- 
stance, from Fig. 1 1 one infers for t ~ ^nonpert that the classical field decay "overshoots" 
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FIGURE 12. Effective particle number density for the transverse modes as a function of time for various 
momenta p < 5po. At early times, modes with p ~ po are exponentially amplified with a rate lyo- Due to 
nonUnearities, one observes subsequently an enhanced growth with rate 670 for a broad momentum range. 
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FIGURE 13. Same as in Fig. 12, for the longitudinal modes. Nonlinear source effects trigger an 
exponential growth with rate 470 for p < 2po. The thick line corresponds to a mode in the parametric 
resonance band, and the long-dashed line for a similar one outside the band. The resonant amplification is 
quickly dominated by source-induced particle production. 



and is temporarily reversed by feed-back from the modes. This can be directly seen in 
the evolution for the rescaled field shown in Fig. 14, which shows a dip around ?nonpert- 
(The particle numbers of Figs. 12 and 13 exhibit correspondingly a reverse behavior.) 
The oscillations in the envelope of o{t) damp out with time, and one observes a slow 
decay of the field and the associated energy at later times. The latter phenomena cannot 
be seen in leading-order or Hartree-type approximations and it is crucial to include the 
next-to-leading order contributions (cf. also Sec. 4.1.1). They are important for a reliable 
description of the system at the end of the resonance stage for finite < 1/A. For real- 
istic inflationary models with typically A ^ 1 this is, in particular, crucial to determine 
whether there are any radiatively restored symmetries. 
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FIGURE 14. The rescaled field (J as a function of time. 



The characteristic properties described above can be understood analytically from the 
evolution equations for the one- and two-point functions. To set the scale we use the 
initial longitudinal mass squared = M^{t = 0) with (cf. Eq. (146)): 



9 A 

m^-\ 

6A^ 



371,(0 + (^-1)71(0 



(225) 
(226) 



where we denote the "tadpole" contributions by 7[| and 7^ from the longitudinal 
and transverse propagator components, respectively, and some A ^ po- Initially, 



7^|(0,0;p) = 1/ [2ft)||(p)] with frequency (0\\ 



+ m2 ( 1 + 3o^)] 1/2, and simi- 



larly for the transverse components F^(0,0;p) where the frequency tt>^(p) contains 
a mass term Oq instead of 3aQ. Parametrically the initial statistical propagator is of 
order one, i.e. ~ ff{N^X'^). Parametric resonance leads to the dominant amplification 
of 7^i(?,?;po) with rate IJq. As a consequence of the exponential amplification of 
the statistical propagator there is a characteristic time at which loop corrections will 
become of order one as well. This is schematically summarized in Fig. 15. The time 

^A~i/2) is denoted by f = ^source- At this time the one-loop 



when F^(f,f;po) - 

diagram with two field insertions indicated by crosses as depicted in Fig. 15 will give a 
contribution of order one to the evolution equation for (?,?;p). For instance, the two 
powers of the coupling coming from the vertices of that diagram are canceled by the 
field amplitudes (221) and by propagator lines associated to the amplified i^^(?,?';po). 
Similarly, at the time t = ^collect the maximally amplified transverse propagator mode has 
grown to F^{t,t';po) ~ ^{N^^^X^^/^). As a consequence, the "setting sun" diagram 
in Fig. 15 becomes of order one and is therefore of the same order as the classical 
contributions. Though the loop corrections become of order one later than the initial 
time, they induce amplification rates that are multiples of the rate 70 which lead to a 
very rapid growth of modes in a wide momentum range. Finally, when the fluctuations 
have grown nonperturbatively large with F±{tj';'po) 



*^A^i) any loop correction 
will no longer be suppressed by powers of the small coupling A . In this case the nonper- 
turbative l/N expansion becomes of crucial importance for a quantitative description of 
the dynamics to later times. In the following, we derive these results from the evolution 
equations (183)-(189) in more detail. 
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(IV) 



tnonpeit ~ 111 ( A ^2 y, 

F^~0(NOA-1) 



O(A-i) 



O(NOAO) 



(III) 

tcollect ' 



2 Unpen /3 +ln(N)/67o 



F^~0(n1^^A~^^^) for N< a 



-1 



~ O(NOaO) 

rate: 6 7^ for F^(p9^Po) 



(11) 

tsource 

F^~0(N"A-i'2) 



tnonpeit 

rO \-l/2x 



O(NOaO) 

rate: 4 7^ for I^,(p<2Pq) 



(I) 




F^(t,t;pQ)~exp(27gt) 


rate: 2 7^ 



t = o , f^~o(noaO) 



FIGURE 15. Schematic overview of the characteristic time scales and the respective relevant diagram- 
matic contributions (see text). 



(I) Early-time (linear) regime: Resonant amplification. At early times the cr-field 
evolution equation receives the dominant, i.e. ^(A*^), contributions from the classical 
action S given in (13). As a consequence, the field dynamics can be described by the 
classical equation of motion. In the following, all quantities are rescaled with appropriate 
powers of Mq to become dimensionless and we set Mq = 1 . The classical field equation 
reads 

d^G{t) + o{t) + a^{t)=0. (227) 

Differential equations of this time have been extensively studied in the litera- 
ture. For the initial condition considered here it has the (anti-)periodic solution 
o{t + n/coo) = —o{t), which can be expressed in term of the Jacobian cosine cn 

as a{t) = aocn[f y''l + cTq, Gq/ ■\J2{1 + cTq)]. Therefore, the solution shows a rapidly 
oscillating behavior with characteristic frequency 

(00 ^ 2^\+o^. (228) 
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The period average of the field amplitude can be also expressed in terms of the initial 
amplitude ctq as C72(r) ~ c7q/2. 

The evolution equations for the two-point functions to order correspond to 

free-field equations with the addition of a time-dependent mass term ~ 3a^(?) for the 
longitudinal and ~ <7^(?) for the transverse modes: 

d^ + p^ + l + a\t)]F^{ty;p) = 
<92 + p2 + i+3a2(0]i^|(r,r';p) = 0, (229) 

and equivalently for the spectral functions p±{tj';'p) and P||(?,?';p). As a consequence 
of this approximation, the two-point functions can be factorized as products of single- 
time functions: 

ia(f,/;p) = \[Mt;p)fl{t';p)+fl{t;p)Mt';p)], 

Pi(r,/;p) = /[/±(?;p)/I(r';p)-/I(?;p)/_L(?';p)] , (230) 

and similarly for the longitudinal components. We emphasize that the simple decom- 
position (230) is no longer correct at higher orders in the coupling. In terms of these 
so-called mode functions the equations of motion read, e.g. for the transverse modes: 

'd^ + p^ + l + a-{t)]Mt;p) =0. (231) 

Up to an overall arbitrary phase, the above initial conditions for the two-point functions 



translate into /j_(0,p) = l/^2ft)^(p), (??/±(?,p)|f=o = — w^(p)/2 and equivalently 
for For this approximation the quantum field theory problem becomes mathemati- 
cally equivalent to the well-known classical mechanics problem described in the above 
introduction. The analytical solution of the Lame-type equation (23 1) is well known and 
can be summarized for our purposes as follows: There is an exponential amplification 
of F^{t,t';'p) for a bounded momentum range < p^ < <7q/2. This corresponds to a 

maximum momentum for amplification p^^^^ + 1 + 0'^(?) ~ 1 + ag = ((yo/2)^. A further 
important result is that there is a separation of scales: (Oq ^ 7o, with the maximum am- 
plification rate 70 ~ 25(0d for p^ = Pq ~ ao/4 with the small number 5 < e^''^ ~ 0.043. 
There is much smaller growth in a narrow momentum range for F\^, which is of no 
importance here since it is quickly dominated by loop corrections as is shown below. 
Time-averaged over ~ 0)^^ , for t,t' ^ Jq^ one finds the result of Eq. (224). Note here 
that also leading-order large-A'^ or Hartree-type approximations show the same bounded 
amplification regime, such that at all times no higher momentum modes would get pop- 
ulated. This is, of course, an artefact of the approximation which is absent once NLO 
corrections are taken into account as is described below. 

The analytic results for the Lame regime agree precisely with the NLO numerical 
results for sufficiently early times. Figs. 12 and 13 show the transverse and longitudinal 
particle numbers for various momenta in the range < p < 5pQ, averaged over the rapid 
oscillation time ~ 1/ coq- For the transverse modes the lowest two momenta shown are 
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inside the resonance band. We define the effective particle numbers as in Eq. (201).^^ As 
shown in Sec. 4.1.2 for the 0{N) model these definitions yield an efficient description 
approaching a Bose-Einstein distributed particle number at sufficiently late times. 

(II) Source-induced (nonlinear) amplification regime: Strongly enhanced particle pro- 
duction for longitudinal modes. The ^(A^) approximation for longitudinal modes 
breaks down at the time 

t^t' = Wee : F^{t,t';po) ~ ^iN^^X-'^^) . (232) 
This can be derived from the ^(A) evolution equations to which one-loop self-energies 



contribute, which diagrammatically are given by V y . The approximate evolu- 

tion equation reads: 

(^d^ + p^+M\t)+3o\t)y\\{t,t';p)c^ 

^^^^^a(o|/Jd/'a(On^(M'';p)ii|(^'',^';p) 

— J^ d/'a(On^(r,r";p)P||(/',r';p)| = RHS. (233) 

Here we have abbreviated Uj_{t,t";-p) = J d^q/ (In)^ F^{t,t";-p - q)A^{tj";q) with 
A = {F, p}, and we used that ^ Dn and ^ pj_-^'^ One observes that indeed for 

~ ^(A^°A~^/2) the r.h.s. of Eq. (233) becomes ~ and cannot be neglected. 
In order to make analytical progress, one has to evaluate the "memory integrals" in 
the above equation. This is dramatically simplified by the fact that the integral is ap- 
proximately local in time, since the exponential growth lets the latest-time contributions 
dominate the integral. (This will be the case for times t < ^nonpert after which exponential 
amplification stops, cf. below.) For the approximate evaluation of the memory integrals 
we consider time-averages over (Oq^ '^Jq^ '■ 



t rt 

II 



At" — ^ / dt" (c~ 1) (234) 

Jt-c/(OQ 

and perform a Taylor expansion around the latest time t (t'): 

P||,j_(?,/';p) ~ d,.p^^^^{t,t"-T^)l=,.{t"-t) = {t-t"), 

Fi±{^,t";p) - Fi^{t,t;p), (235) 



We have employed 2 = as in Sec. 4.1.2 which is, however, irrelevant. 

The latter inequality will be discussed in detail in Sec. 5, where it is shown to indicate the validity of 
classical statistical field theory approximations. 



Introduction to Nonequilibrium Quantum Field Theory 



83 



where we have used the equal-time commutation relations (117). With these approxima- 
tions the r.h.s. of Eq. (233) can be evaluated as: 



(-2 (M — I) 

RHS ~ Xo^{t)^^ — '-T^{t)FM,t';-p) (mass term) (236) 
cOq 3N 

(N —\) 

+ Xo{t)o{t')^^——^U^{t,t';p) (sourceterm) (237) 
cOq oN 

The first term is a contribution from NLO to the mass, whereas the second term repre- 
sents a source. Note that both the LO mass term and the above correction to this mass 
are of the same order in A, however, with opposite sign. To evaluate the momentum 
integrals, we use a saddle point approximation around the dominant p ~ po, valid for 
r , > \ with {t, t', p) ~ (f , t', Po) exp [- 1 1^' |(? + 1') {p - p^f/l] From this one 
obtains for the above mass term: 

The result can be used to obtain an estimate at what time t this loop correction becomes 
an important contribution to the evolution equation. Note that to this order in A it is cor- 
rect to use F^(?,?';po) ~ ) on the r.h.s. of (238). The condition XT^{t = fnonpert) 
~ ^(1) can then be written for A ^ 1 as: 



^nonpert ~ 



(239) 



The same saddle point approximation can be performed to evaluate the source term 
(237): 

Here we only wrote the source term for p = where it has its maximum, although it 
affects all modes with p < 2pQ. Again this can be used to estimate the time t = ^source at 
which A ~ 



1 

2 ^nonpert 



(241) 



One arrives at the important conclusion that the source term (237) becomes earlier of 
order one than the mass term (236): For ^source ^ ^ ^ ^nonpert the source term dominates 
the dynamics! Using these estimates in (233) one finds that the longitudinal modes with 
< p < 2j!7o get amplified with twice the rate IYq: 

FM,t■^^) ~ AF2(r,r;po) ~ Ae^W. (242) 



30 Here Y{p) ~ yo + /o (P - Po) V2 with /q' ~ -645(1 - 68)^1 + al/al 
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Though the non-linear contributions start later, they grow twice as fast! The analytical 
estimates for ^source and rates agree very well with the numerical results shown in Fig. 13. 



(Ill) Collective amplification regime: explosive particle production in a broad mo- 
mentum range for transverse modes. A similar analysis can be made for the transverse 
modes. Beyond the Lame-type ff{X^) description, the evolution equation for re- 
ceives contributions from the feed-back of the longitudinal modes at ^(A) as well as 
from the amplified transverse modes at ^(A^). They represent source terms in the evo- 
lution equation for Fl(?,?';p) which are both parametrically of the form ~ X^F^/N as 
is depicted below: 





A cf.(242) A2 3 

, \ N 



-F 



Following along the lines of the above paragraph and using (239) this leads to the 
characteristic time t = ^collect at which these source terms become of order one: 



^ 2 InA^ 

^collect — T ^ nonpert "I" "7 

3 670 



(243) 



For ? ~ ~ ^collect the dominant transverse mode has grown to 

ia(r,r';po) ~ ^(A^^/^A"2/3)_ (244) 

Correspondingly, for ^collect ^ ^nonpert one finds a large particle production rate ~ 670 
in a momentum range < p < 3po, in agreement with the full NLO results shown in 
Fig. 12. In this time range the longitudinal modes exhibit an enhanced amplification 
as well (cf. Fig. 13). It is important to realize that the phenomenon of source-induced 
amplification repeats itself: the newly amplified modes, together with the primarily 
amplified ones, act as a source for other modes, and so on. In this way, even higher 
growth rates of multiples of 70 can be observed and the "explosive" amplification 
rapidly propagates towards higher momentum modes. We emphasize that the collective 
amplification regime is absent in the LO large-A^^ approximation. Consequently, even for 
the transverse sector the latter does not give an accurate description if ^collect < ^nonpert, 
thatisforA^< A"^ 



Behavior of the field. Around t < ^nonpert is the earliest time when sizeable correc- 
tions to the classical field equation (227) appear, which give (cf. Eq. 189): 

[d^ + 1 + dM'-{t) + o\t)] a{t) = ~J^dt'i:P{t,t';p = 0)\a=oait') , (245) 

where 8M^ ~ '^^^^ '^ T^. Before ^nonpert, where the "memory expansion" discussed above 
is valid, we can discuss the behavior of the field a = a^^^ + 5a perturbatively in terms of 
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a slowly varying small correction 5a to the classical solution o^^^ of (227) and a small 
5M^. If we neglect for a moment the NLO contributions on the r.h.s. the linearized 
equation (245) then yields 

a{t) ^(l ^-^^^^T^it)] . (246) 

One concludes that there is an exponential decrease of the field amplitude at LO since 

T^{t) ^ e^^ot for t 

^ ^nonpert according to (238). The dominant NLO corrections read 

/' dt' Ef , p = 0) I cy=o ait') ~ (^,EJ , t'; p = 0) | ,,=0) Ur' o (0 
Jt-c/(Oo " 2(0^ V II / 



2(D^18A^y (2;r)3(2;r) 



c 



2 A2 



2(o2 ISiV 



Tl{t)o{t). (247) 



Note that the NLO correction to the effective mass term comes with an opposite sign 
than the LO correction. For t ?nonpert one has 5M^ ~ ^{N^X^) and EJ ~ ff{N-^X^), 

i.e. they become of the same order in A. As a consequence, cancellations may lead 
(temporarily) to reverse field decay. This is indeed what one observes from the full 
numerical NLO solution for A/^ = 4 and X = 10^^ shown in Fig. 14. Around ^nonpert 
strong nonlinearities appear, where the field decay "overshoots" and is shortly reversed 
by feedback from modes, overshoots again etc. 

( IV) Fluctuation dominated regime: nonperturbatively large densities ~ 1 / A and 
quasistationary evolution. When the system evolves in time, more and more energy 
is converted into fluctuations. The description in terms of ff{X^) evolution equations 
break down at 

t^t' = ^nonpert : i^±(?,?';Po) ~ ^(A^° A^^) . (248) 

Transverse and longitudinal modes have grown to ff{N^ X^^) in a wide momentum 
range for times t > ^nonpert- This corresponds to nonperturbatively large particle number 
densities n^(p) and n||(p) inversely proportional to the coupling. Because of this para- 
metric dependence there are leading contributions to the dynamics coming from all loop 
orders. As a consequence, shortly after fnonpert a comparably slow, quasistationary evolu- 
tion sets. It is important to note that descriptions based on a finite-order loop expansion 
of the 2PI effective action cannot be applied. To describe this regime it is crucial to 
employ a nonperturbative approximation as provided by the 2PI \/N expansion at NLO. 

In order to discuss this in more detail, we consider the self-energies appearing in 
the evolution equations (138) for Fl(?,?';p) and pj^(?,?';p). From the NLO expressions 
given in (184) and (185) one finds for E 



F . 



li(M';p) = -JfiJjJ^l (249) 
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If {t,t'; q)FL , p - q) - ^ Ip {tj'; q)p^ , t'; p - q) 
+ Pf {t,t'; q)FL , p - q) - ^ Pp [t, t'; q)p± p - q) } • 

The functions If,p and /V,p contain the summation of an infinite number of "chain" 
graphs, where each additional element adds another loop to the graph (cf. also the figures 
in Sec. 2.4. 1). We will argue in the following that each loop order of this infinite number 
of graphs contributes with the same order in the coupling A. 

According to (181) and (182) the sum of "chain" graphs described by If,p can be 
iteratively generated by the relations 



A A 

Ip{ty;q) = -—nF{t/;q) + — dt" Ip{t,t";q)nt 

~Jn II ^^"^^(^^";q)np(?",?';q), (250) 
/p(r,r';q) = -^np(?,?';q) + A^^'j^"/p(^,^";q)np(/V';q), (251) 



and corresponding expressions for /^.p (cf. Eqs. (186) and (187)). Here the "chain 
elements" TIf,p are given by 



1 f d^k f 

nF{t,t';q) = --J j^lF^\{t,t';q-k)F^\{t,t';k) 



+ {N-l)F^{t/-q-k)F^{t/;k)-- Pi|(r,/;q-k 



Up{t,t';q) 



p^\{t/;k) + iN-\)p^{t/';q-k)p^{t,t";k) 



(252) 



3|i^|(?,?';q-k)p||(r,f';k) 
+ {N-l)F^{t,t';q-k)p^{t,t';k)y 



(253) 



(0 



Denoting a given loop order Ihy ly we can write 



lF{t,t';q) = l^ll!\t,t';q) 



(254) 



il.'\t.t';q) 
iPit.t';q) 



-:^^F{t,t';q). 



+ 



J^dt"np{tr;q)nF{t",t';q) 

J^dt''nF{t/'-,q)np{t''j'-q), 
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and similarly for Ip = L/Li^p^- Concentrating on the transverse modes, one observes 
from (252) and (253) that 

n^(r,r';0) ~ ff{NX~^) , Ylp{t,t'\Q) ~ ff{NX-^) . (255) 

Inserting this into (250) or (254) for If, and into the corresponding expression (251) for 
Ip, one finds 

l^p\t,t';0)^ ff{N^X-^), lll\t,t';0)r^ ff{N^X^), (256) 

irrespective of the loop order I. Here it is important to note that averaged over the rapid 
oscillation time the spectral functions are of order one: 

Pi(r,/;p) ~P||(r,r';p) ~ ^(ivOAO). (257) 

As a consequence, with (257) the self-energy (250) receives leading contributions in the 
coupling from all loop orders. Since these contributions are proportional N^^X'^, they 
are subleading in 1/A''. In particular, the expansion in powers of l/N employed here 
remains a valid nonperturbative approximation scheme. It is interesting to note that the 
2PI l/N expansion to NLO can be understood as a "three-loop" approximation with an 
effective vertex. This is denoted schematically below for the graphs appearing at NLO: 




The effect of the nonperturbatively large densities is taken into account by a self- 
consistent vertex correction as indicated by the diagrammatic equation for the four- 
vertex. The question of how self-consistent vertex corrections can be systematically 
described for cases where no 1 /N expansion is applicable will be treated in Sec. 6. We 
finally note that for times t > fnonpert the evolution equations are no longer characterized 
by the effective time-locality in the sense described above such that the "memory 
expansion" cannot be used for an approximate description at late times. 



5. CLASSICAL ASPECTS OF NONEQUILIBRIUM QUANTUM 

FIELDS: PRECISION TESTS 

It is an important question to what extend nonequilibrium quantum field theory can be 
approximated by classical statistical field theory. It is a frequently employed strategy 
in the literature to consider nonequilibrium classical dynamics instead of quantum dy- 
namics since the former can be simulated numerically up to controlled statistical errors. 
Classical statistical field theory indeed gives important insights when the number of field 
quanta per mode is sufficiently large such that quantum fluctuations are suppressed com- 
pared to statistical fluctuations. We will derive below a sufficient condition for the va- 
lidity of classical approximations to nonequilibrium dynamics. The description in terms 
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of spectral and statistical correlation functions as introduced in Sec. 3.4.1 is particu- 
larly suitable for comparisons since these correlation functions possess a well-defined 
classical limit. However, classical Rayleigh- Jeans divergences and the lack of genuine 
quantum effects — such as the approach to quantum thermal equilibrium characterized 
by Bose-Einstein or Fermi-Dirac statistics — limit the use of classical statistical field 
theory. To find out its use and its limitations we perform below direct comparisons of 
nonequilibrium classical and quantum evolution for same initial conditions. One finds 
that classical methods can give an accurate description of quantum dynamics for the case 
of large enough initial occupation numbers and not too late times, before the approach to 
quantum thermal equilibrium sets in. Classical approaches are unsuitable, in particular, 
to determine thermalization rates. 

Classical methods have been extensively used in the past to rule out "candidates" 
for approximation schemes applied to nonequilibrium quantum field theory. Approx- 
imations that fail to describe classical nonequilibrium dynamics should be in general 
discarded also for the quantum case. If the dynamics is formulated in terms of correla- 
tion functions then approximation schemes for the quantum evolution can be straight- 
forwardly implemented as well for the respective classical theory. For instance, the 2PI 
l/A/^-expansion introduced in Sec. 2.4 can be equally well implemented in the classi- 
cal as in the quantum theory. Therefore, in the classical statistical approach one can 
compare NLO results with results including all orders in 1/A^. This gives a rigorous an- 
swer to the question of what happens at NNLO or beyond in this case. In particular, for 
increasing occupation numbers per mode the classical and the quantum evolution can 
be shown to approach each other if the same initial conditions are applied and for not 
too late times. For sufficiently high particle number densities one can therefore strictly 
verify how rapidly the \/N series converges for far-from-equilibrium dynamics! The 
possibility of precision tests are an important aspect of classical statistical field theory 
methods. 



5.1. Exact classical time-evolution equations 

In this section we define the basic classical correlation functions and derive exact 
dynamical equations for them. The evolution equations for the classical spectral function 
and the classical statistical propagator will turn out to be of the same form as Eqs. (138) 
for the corresponding quantum correlators — with the only difference that the quantum 
self-energies are replaced by the respective classical ones. The comparison between 
quantum and classical evolution can be very clearly discussed using the the language 
of correlation functions, since the quantum spectral function p{x,y) and the statistical 
propagator F(x,);) of Eq. (116) both have a well-defined classical equivalent. 

We consider the classical A''-component scalar field (pa with action (13). The classical 
field equation of motion is then given by 



U^ + m^ + —G>i,{x)(pb{x 
bis 



(pa{x) = 0, 



(258) 



whose solution requires specification of the initial conditions ^a(0,x) = (Pa{x) and 
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7ta{0,x) = 7ta{x) with 7ta{x) = d^o(Pa{x). We define the macroscopic or average classical 
field: 

<^a,ciW = (<P«W)ci = J DKD(pW[n, (p](pa{x) . (259) 

Here iy[;r, 9] denotes the normalized probability functional at initial time. The measure 
indicates integration over classical phase-space: 

/ DnDcp = / n nd^«(x)d(Pa(x) . (260) 

Here the theory will be defined on a spatial lattice to regulate the Rayleigh- Jeans diver- 
gence of classical statistical field theory. The connected classical statistical propagator 
Fab,c\{^,y) is defined by 

Fab,ci{^,y) + ^aA^)^b,ci{y) = Mx)(pb{y))c\ = J DTtDcp W [7t, (p](Pa{x)(Pb{y) . (261) 

The classical equivalent of the quantum spectral function is obtained by replacing — z 
times the commutator with the Poisson bracket:^^ 

Pab,c\{^,y) = -( {<PflW, 9i(j)}PoissonBracket )cl • (263) 

As a consequence, one finds the equal-time relations for the classical spectral function: 

PabA^^y) Ix0=y0 = 0, d^OpabA^^y) 1x0=/ = SabS{x - y) . (264) 

Though their origin is different, note that they are in complete correspondence with the 
respective quantum relations (117). 

The evolution equations for the classical statistical correlators can be obtained by 
functional methods in a similar way as for the quantum theory. In order not to be too 
redundant, we employ here a different approach starting from the differential equation 
for the free correlators. We will choose here VK[;r,(p] to be invariant under the 0{N) 
symmetry with 0a = such that 

Fabxi (x, y) = Fci {x, y) dab , (265) 

and equivalently for pabA^^y)-^^ ^^^^ discuss 0^ 7^ below. The unperturbed spec- 
tral function is a solution of the homogeneous equation^^ 

[a + m2]po(x,);)=0, (266) 



^ Recall that the Poisson bracket with respect to the initial fields is 

N . 

{A(x),B(y)}poissonBracket = E / ^''^ 

^ An explicit example for W[n, <p\ would be W[k, <p] = exp 

{/d''x(7r2 + (V9)2 + m2^2j /2} where 

Zq = / DnD^>W\%, We emphasize that none of the following derivations will make use of a specific 
form for W[7r,(p]. 

This is the case both in the quantum and the classical theory. 



5A(jc) 5A(x) 



(262) 
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with initial conditions determined by the equal-time canonical relations (264). Let Fq 
denote the solution of the unperturbed homogeneous problem 



[D, + m^]FQ{x,y)=0, 



(267) 



with initial conditions determined by the initial probability functional. The derivation of 
the time evolution equations of classical correlation functions is conveniently formulated 
by introducing two additional two-point functions: the classical retarded and advanced 
Green functions, which are related to the classical spectral function Pd by 



(268) 



The classical retarded self-energy Hn^c^ix.y) can be defined as the difference between 
the full and free inverse retarded Green functions 



Y.R^,,{x,y) = {G'^y {x^y) - {Gl)-\x,y) . 
where the free retarded Green function 

G^(x,3;) = 0(/->'>o(x,>'), 
solves the inhomogeneous equation 

[U, + m^]Gl{x,y) = 5'^+\x-y), 



(269) 



(270) 



(271) 



with retarded boundary conditions. In the same way we define the advanced self-energy 
Ea,c1- Retarded (advanced) Green functions and self-energies vanish when < y^ 
(x^ > y^). With these definitions we can rewrite (269) and the respective equation for 
the advanced Green function as: 



y^R /^R V 



Cl7 



GA /-'A /-^n V 



where we use a compact notation 

AB = J d''+hA{x,z)B{z,y). 



(272) 
(273) 



(274) 



A combination of Eqs. (272) and (273) gives the Schwinger-Dyson equation for the 
classical spectral function pd = G^j — G^j with 



pel = Pq — Gq- Lr^cI ■ G^i + Gq ■ Ea,c1 ■ G^i . 



(275) 



There is a similar identity for the classical statistical propagator F^, defined in (261). 
Using the definitions (272) and (273) one can write the following identity for the 
statistical function 



"cl 



Fq-Gq- 



.-1 



^cl 



■ Fc\ - Fq ■ 



+ Gi 



1-1 



■G 



cl- 
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We can now define a classical statistical self-energy as 



^F,ci — —Gc\ -Fci-G^i +Gq -Fq-Gq , (277) 

and find the Schwinger-Dyson equation for the statistical propagator: 

Fc\ = Fq — Gq- Lr^cI ■ Fc\ — Fq ■ ^A,c\ ■ — Gq ■ Lf^ci ■ G^i . (278) 

Acting with (D + m^) on (275) and (278) brings the classical Schwinger-Dyson equa- 
tions in a form which is more suitable for initial-value problems. We make the retarded 
nature of rnanifest by writing 

Z^,ei(^,y) = i:^\x)5'+\x-y) + @{x^ -/)Ep,d(^, , (279) 

and similarly for E^ cil-^jj)- The spectral component of the classical self-energy is 
^pM^^y) = ^RAi^^y) -'^A,ci{^,y) = -Ep,ci(y,-^)- After properly taking into account 
all ©-functions as well as (266), (271) and (267), one finds for the exact time evolution 
equations for and F^: 

[D,+M^,{x)]p,i{x,y) = - dzEp,d(x,z)Pci(z,y), (280) 

[n, + M^^{x)]F,i{x,y) = - dzEp,ei(x,z)Fci(z,);) 

+ / dzEf,ci(^,z)pci(z,y), (281) 
Jo 

where we use the abbreviated notation fl'^ dz = fl^ dz'^ and 

m2(jc) =m2 + E^°)(x). (282) 

We emphasize that the form of these classical time evolution equations is the same as for 
the quantum evolution described by (138). If the initial conditions are chosen to be the 
same then the only difference between classical and quantum theory originates from the 
self-energies entering the evolution equations. As a direct consequence, one concludes 
that LO large-N or any Gaussian/Hartree type approximation exhibits purely classical 
dynamics: since the quantum self-energies E^' and Ep as well as the classical E^^ d and 
Ep ci vanish for these approximations the evolution equations derived from the quantum 
theory are identical to those derived from the classical theory. The tadpole contributions 

E*^*^) and l!"^^ entail no difference if the initial conditions are chosen to be the same. In the 
following we consider approximations with non-vanishing spectral and statistical self- 
energies in order to discuss the differences between classical and quantum dynamics. 

5.2. Classicality condition 

The classical self-energies Ep ^i, T,pci{x,z) and E^^-* could be approximated by pertur- 
bation theory to a given order in the coupling A . However, as for the case of the quantum 
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evolution the perturbative classical self-energies lead to a secular time-evolution and fail 
to provide a reliable description for the classical nonequilibrium dynamics. The ana- 
lytical description can, however, be based on the same 2PI summation techniques as 
introduced for the quantum theory above. We will only state here the result of the re- 
spective calculation within classical statistical field theory, since we will then discuss 
how the same result can be directly obtained from the classical limit of the known quan- 
tum self-energies. 

One finds for the classical statistical 0{N) model employing the 2PI l/N expansion 
to next-to-leading order: 

4T^(-^'3') = -^^ci(^,j)/F,ci(-^,y), (283) 

^fX^\^^y) = -^(^cl(^,j)/p,cl(^,j)+Pcl(^,j)/F,cl(^,j) 



with the classical summation functions 
A 
6 



6 



d^zi J dz''lp,ci{x.z)Fl{z,y)-2j dz"" lFAMF,,{z,y)p,x{z.y) 



Ip,ciM = -^F^iMp^iix.y)-- J d^'zj dz°Ip,ciix,z)F^iiz,y)pciiz,y). (284) 

At this stage we can compare the result from the 2PI l/N expansion for the quantum 
theory, as discussed in Sec. 2.4, with the resummed classical expressions obtained 
here. One observes that the classical self-energies (283) — (284) are obtained from the 
corresponding quantum expressions (178) — (182) by dropping the terms with a product 
of two spectral functions p compared to a product of two statistical propagators F. 
Accordingly, the classical self-energies are obtained from the respective quantum ones 
as 

Ef,ei = Ef(F2>p2) , Ep,ci = Ep(F2>p2). (285) 



The same analysis can be done employing the 2PI loop expansion. For the follow- 
ing analytical discussion we consider the spatially homogeneous case and employ the 
Fourier modes E/7ci(f,f';p) and Ep^cl(^'^';p)• For the classical statistical 0{N) symmet- 
ric theory one obtains to two-loop order the self-energies: 

4'ir^(^^';P) = -A2^|^^Fd(r,r';p-q-k)Fei(r,r';q)Fei(^?';k), 
^fT^i^^t'-^l?) = -^'^ |^^Pci(?,?';p-q-k)Fei(?,?';q)Fei(^?';k).(286) 
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Let us directly confront this with the respective two-loop self-energies for the respective 
quantum 0{N) theory: 



, (21oop) 



t(21oop) 



{t,t';p) 



,2 N + 2 
18A^2 



/ F(f,f';p-q-k) 

Jq.k 



F{t,t';q)F{t,t';k)--p{t,t';q)p{t,t';k) 



/ p(r,?';p-q-k) 

Jq.k 



F{t, t'- q)F , t'- k) - ^ P , t'; q)p (r , k) 



(287) 



From this one infers a sufficient condition for classical evolution: 



\F{t,t';q)F{t,t';k)\ > - |p(r,/;q)p(?,/;k) 



(288) 



This condition has to be fulfilled for all times and all momenta in order to ensure that 
classical and quantum evolution agree. However, we will observe below that it can be 
typically only achieved for a limited range of time and/or momenta. 

One expects that the classical description becomes a reliable approximation for the 
quantum theory if the number of field quanta in each mode is sufficiently high. The 
classicality condition (288) entails the justification of this expectation. In order to il- 
lustrate the condition in terms of a more intuitive picture of occupation numbers, we 
employ the free-field theory type form of the spectral function and statistical propagator 
given in Eq. (211). From this one obtains the following estimates for the time-averaged 
correlators: 



F^it,t';p) 



p2(f,f';p) 



27t Jt-lyT/cOp 



dt'F^{t,t';-p) 



Mt) + \/2f 
2ml{t) 



2(ol{t) 



(289) 



Inserting these estimates in (288) for equal momenta yields 

2 o 



»4 



or 7ip(r) > 0.37. 



(290) 



This limit agrees rather well with what is obtained for the case of thermal equilibrium. 
For a Bose-Einstein distributed particle number hbe = (e^/^ ~ 1)^ with temperature 
T one finds n^^{oi = T) = 0.58, below which deviations from the classical thermal 
distribution become sizeable. We emphasize that in thermal equilibrium the fluctuation- 
dissipation relation (129) ensures \F^^°^\o),i^)\ 3> |p^'^'i^(ft),p)| for high-temperature 
modes T ^ (o: 



F('^'i)(a),p) = -/(nBE(«) + ^) P^'^'i^w, 



p) ^ -i — p^^"^' 

0) 



(w,p). 



(291) 



Introduction to Nonequilibrium Quantum Field Theory 



94 



The nonequilibrium estimate (290) as well as the equilibrium result (291) leads one to 
expect that the quantum evolution described by F and p is well approximated by the 
classical evolution in terms of and pd for large initial occupation numbers np(? = 0). 
However, note that the nonequilibrium classicality condition (288) is typically not ful- 
filled at all times, since the unequal-time correlator F(f,f';p) can oscillate around zero 
with a phase difference to the oscillations of p(f,f';p) (cf. also the free-field expres- 
sions (211)). In particular, we will find that characteristic quantities for the quantum 
late-time behavior such as thermalization times are very badly approximated by the re- 
spective classical estimates. Nevertheless, for time-averaged quantities and not too late 
times t <C ?eq, i-C- before the approach to quantum thermal equilibrium sets in, the clas- 
sical evolution can give an accurate description for sufficiently high initial occupation 
numbers as will be demonstrated below. 

A similar analysis can be performed in the presence of a non-vanishing classical 
macroscopic field <^a,c\{^) defined in (259). For instance, to NLO in the 2PI 1/A'' ex- 
pansion the classical macroscopic field is described by the evolution equation 




= - r d3;E^i„^(x,);; = 0, i^i, Pel) My) ^ (292) 

J 

where ^ciai('^'^ci = O^-^ci) is given by (146) with the replacement F — > F^x, and 
^^(x,3;;0ci = Oj-^chPci) is obtained from the respective quantum self-energy 
Ep ^^(x,};; ^ =0,F,p) given in (185) by employing (285). This can be directly compared 
to the corresponding evolution equation (189) for the case of the quantum field theory. 



5.3. Precision tests and the role of quantum corrections 

The nonequilibrium evolution of classical correlation functions in the 0{N) model 
can be obtained numerically up to controlled statistical errors. Initial conditions for 
the nonequilibrium evolution are determined from a probability functional on classi- 
cal phase-space. The subsequent time evolution is solved numerically using the classical 
equation of motion for the field (258). The results presented below have been obtained 
from sampling 50000-80000 independent initial conditions to approximate the exact 
evolution of correlation functions. The statistical propagator is constructed from these 
individual runs according to (261) and using (263) for the spectral function. Since these 
results include all orders in l/N, they can be used for a precision test of the 2PI 1/N 
expansion implemented in classical statistical field theory. We emphasize that this com- 
pares two very different calculational procedures: the results from the simulation involve 
thousands of individual runs from which the correlators are constructed, while the cor- 
responding results employing the 2PI l/N expansion involve only a single run solving 
directly the evolution equation for the correlators. The accuracy of the simulations man- 
ifests itself also in the fact that the time-reversal invariant dynamics can be explicitly 
reversed in practice for not too late times. The close agreement between full and 2PI 
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FIGURE 16. Left: Unequal-time two-point function F(f,0; p = 0) at zero momentum for = 2, 10,20. 
The full lines show results from the NLO classical evolution and the dashed lines from the exact classical 
evolution (MC). One observes a convergence of classical NLO and exact results already for moderate 
values of A^. Right: Damping rate extracted from F{t,Q;p = 0) as a function of Open symbols 
represent NLO and exact classical evolution. The quantum results are shown with full symbols for 
comparison. The initial conditions are characterized by low occupation numbers so that quantum effects 
become sizeable. One observes that in the quantum theory the damping rate is reduced compared to the 
classical theory. (All in units of niR.) 



NLO results, which will be observed below, is therefore also indicative of the numerical 
precision of the latter. 

We consider a system that is invariant under space translations and work in momen- 
tum space. We choose a Gaussian initial state such that a specification of the initial 
two-point functions is sufficient. As mentioned above, the classical spectral function at 
initial time is completely determined from the equal-time relations (264). For the sta- 
tistical propagator we take F(0, 0; p) = [«p(0) + j]/o}p, with the initial particle number 
representing a peaked "tsunami" (cf. Sect. 4.1.1). 

The initial mode energy is given by (Op = [p^ +M^j)^/^ where Mc\ is the one-loop 
renormalized mass in presence of the nonequilibrium medium, determined from the one- 
loop "gap equation" for in Eq. (283). As a renormalization condition we choose the 
one-loop renormalized mass in vacuum mr = M|„(q)^() = 1 as our dimensionful scale. 

The results shown below are obtained using a fixed coupling constant A/m^ = 30. 

On the left of Fig. 16 the classical statistical propagator Fci(?,0; p = 0) is presented 
for three values of N. All other parameters are kept constant. The figure compares the 
time evolution using the 2PI l/N expansion to NLO and the full Monte Carlo (MC) 
calculation. One observes that the approximate time evolution of the correlation function 
shows a rather good agreement with the exact result even for small values of A^^ (note 
that the effective four-point coupling is strong, X/6N = 2.5m| for N — 2). For N = 20 
the exact and NLO evolution can hardly be distinguished. A very sensitive quantity for 
comparisons is the damping rate 7, which is obtained from an exponential fit to the 
envelope of Fc\{t,0;p = 0). The systematic convergence of the NLO and the Monte 
Carlo result as a function of 1 /N can be observed from the right graph of Fig. 16. The 
quantitatively accurate description of far-from-equilibrium processes within the NLO 
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approximation of the 2PI effective action is manifest. 

The right graph of Fig. 16 also shows the damping rate from the quantum evolution, 
using the same initial conditions and parameters. One observes that the damping in 
the quantum theory differs and, in particular, is reduced compared to the classical 
result. The effective loss of details about the initial conditions takes more time for the 
quantum system than for the corresponding classical one. In the limit N ^ oo damping 
of the unequal-time correlation function goes to zero since the nonlocal part of the self- 
energies vanishes identically at LO large-A'^ and scattering is absent. In this limit there 
is no difference between evolution in a quantum and classical statistical field theory. 
(Cf. Sec. 4.1.1 and the discussion above.) 

For finite A'^ scattering is present and quantum and classical evolution differ in general. 
However, as discussed in Sec. 5.2, the classical field approximation may be expected 
to become a reliable description for the quantum theory if the number of field quanta 
in each field mode is sufficiently high. We observe that increasing the initial particle 
number density leads to a convergence of quantum and classical time evolution at not 
too late times. In Fig. 17 (left) the time evolution of the equal-time correlation function 
F{tj;p) is shown for several momenta p and A'^ = 10. Here the integrated particle 
density / ^np{0)/Mci = 1.2 is six times as high as in Fig. 16. At p = 2pts one finds 
np=2p^^{0) ~ 0.35 and a slightly larger value at this momentum of about ^ 0.5 for late 
times. For these initial conditions the estimate (290) for the classicality condition (288) 
is therefore approximately valid up to momenta p ~ 2pi^, and one indeed observes from 
the left of Fig. 17 a rather good agreement of quantum and classical evolution in this 
range. For an estimate of the NLO truncation error we also give the full (MC) result for 
N — 10 showing a quantitative agreement with the classical NLO evolution both at early 
and later times. 



5.3.7. Classical equilibration and quantum thermal equilibrium 

From the left of Fig. 17 one observes that the initially highly occupied "tsunami" 
modes "decay" as time proceeds and the low momentum modes become more and more 
populated. At late times the classical theory and the quantum theory approach their 
respective equilibrium distributions. Since classical and quantum thermal equilibrium 
are distinct, the classical and quantum time evolutions have to deviate at sufficiently late 
times. Figure 17 shows the time dependent inverse slope parameter 



which has been introduced in Sec. 4.1.2 to study the approach to quantum thermal 
equilibrium.^^ It employs the effective particle number np{t) defined in (201)^^ and 
mode energy ep{t) given by (203). Initially one observes a very different behavior of 



Note that dhog(n-p\t) + l)/d£p{t) ^T-\t,p). 
Here we employ g(f,f';p) = 0, cf. Sec. 4.1.2. 




(293) 
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FIGURE 17. Left: Nonequilibrium evolution of the equal-time two-point function F{t,t; p) foiN = 10 
for various momenta p. One observes a good agreement between the exact MC (dashed) and the NLO 
classical result (full). The quantum evolution is shown with dotted lines. The integrated initial particle 
density is six times as high as in Fig. 16. Right: A very sensitive quantity to study deviations is the 
time dependent inverse slope T{t,p) defined in the text. When quantum thermal equilibrium with a Bose- 
Einstein distributed particle number is approached all modes get equal T{t,p) = Teq, as can be observed 
to high accuracy for the quantum evolution. For classical thermal equilibrium the defined inverse slope 
remains momentum dependent. 



T{t,p) for the low and high momentum modes, indicating that the system is far from 
equilibrium. The quantum evolution approaches quantum thermal equilibrium with a 
momentum independent inverse slope Teq = A.l mr to high accuracy. In contrast, in the 
classical limit the slope parameter remains momentum dependent since the classical 
dynamics does of course not reach a Bose-Einstein distribution. 

To see this in more detail we note that for a Bose-Einstein distribution, nBE(£p) = 
l/[exp(ep/req) — 1], the inverse slope (293) is independent of the mode energies and 
equal to the temperature Teq. During the nonequilibrium evolution effective thermaliza- 
tion can therefore be observed if T{t,p) becomes time and momentum independent, 
T{t,p) Teq. This is indeed seen on the left of Fig. 17 for the quantum system. If the 
system is approaching classical equilibrium at some temperature and is weakly cou- 
pled, the following behavior is expected. From the definition (201) of np{t) in terms in 
two-point functions, we expect to find approximately 

T{t,p)^T,i{\-e^p/T^,), (294) 

i.e. a remaining momentum dependence with T{t,p) < T{t,p') if Cp > e^. Indeed, this 
is what one observes for the classical field theory result in Fig. 17. 

For a classical theory a very simple test for effective equilibration is available. An 
exact criterion can be obtained from the classical counterpart of the "KMS" condition for 
thermal equilibrium discussed in Sec. 3.4.2. In coordinate space the classical equilibrium 
"KMS" condition reads 

^:lo'^r"(--y)=-pir'(--y). (295) 
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FIGURE 18. Nonequilibrium time evolution in 1 + 1 dimensions in the classical limit of the three- 
loop approximation of the 2PI effective action for one scalar field, = 1 . Shown is the effective mode 
temperature Tc\{t,p), defined in the text, with initial T^\{0^p) /m^ ~ 5 for all p. The classical field theory 
approaches classical equihbrium, Tci{t,p) Td, after a very long time. 



and in momentum space 



F}r\k) = -ina{k^)pf>{k), n,,{k^) = ^. (296) 



Differentiating Eq. (295) with respect to j*^ at jc^ = = ? gives 



1 d d 



(eq). 



(297) 



Combining this KMS relation with the equal-time condition (264) for the spectral 
function leads to 

dtd,F^^''\ty-x-y)\r=, = T,i5{x-y). (298) 

In terms of the fluctuating classical fields na{x) = d^)(pa{x) this is of course the well- 
known equilibrium relation (;ra(?,x);r^(?,y))^j'^^ = rci5(x — y)^^^,. Out of equilibrium 
one can define an effective classical mode temperature 

T,^{t,p) = dtdMt,t';p)\r=,> , (299) 

and effective classical equilibration is observed if T^i {t,p) becomes time and momentum 
independent, Tc\{t,p) Tc\. 

Apart from the 2PI 1/A^ expansion the late-time behavior can also be studied from the 
2PI loop expansion. For the quantum theory this has been demonstrated in Sec. 4. 1 .4 for 
the three-loop approximation of the 2PI effective action (cf. Fig. 6, left). The equivalent 
approximation in the classical theory is given by the evolution equations (280) and (281) 
with the two-loop self-energies (286). In order to demonstrate this, we will use the loop 
approximation to verify explicitly the above statements about classical equilibration. 
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In Fig. 18 the nonequilibrium evolution of the effective "mode temperature" Tc\{t,p) 
is shown for various momentum modes in the (1 + 1) -dimensional classical scalar field 
theory for A'^ = 1 . The equations are solved by a lattice discretization with spatial lattice 
spacing niRa = 0.4, time step ao/a = 0.2, and Ns — 24 sites. Without loss of generality 
we use X/m\ = 1. For the initial ensemble we take Fci(0,0;p) = Tq/{p^ + ^^id 
dtdtiPc\{U^'\p)\t=t'=Q = To with To/rtiR = 5. We have observed that at sufficiently late 
times the contributions from early times to the dynamics are effectively suppressed. 
This fact has been employed in Fig. 18 to reach the very late times. One sees that 
at sufficiently late times the system relaxes towards classical equilibrium with a final 
temperature Td/mj^ ^5.5. For the zero-momentum mode we find an exponential late- 
time relaxation of Tc\{t,p = 0) towards T^i with a small rate (~ 2 x lO^^m^). It should 
be emphasized that typical classical equilibration times are substantially larger than the 
times required to approach thermal equilibrium in the respective quantum theory. In 
contrast to the quantum theory statements about equilibration times are, however, in 
general not insensitive to the employed momentum regularization for the classical theory 
because of the classical Ray leigh- Jeans divergences. 



6. A^-PARTICLE IRREDUCIBLE GENERATING 
FUNCTIONALS II: EQUIVALENCE HIERARCHY 

To understand the success and, more importantly, the limitations of expansion schemes 
based on the 2PI effective action we consider in this section n-particle irreducible (nPI) 
effective actions for n> 2. Recall that the description of the 2PI effective action r[0,D] 
employs a self-consistently dressed one-point function and two-point function, which 
for notational purposes we denote here by D:^^ The one-point and two-point functions 
are dressed by solving the equations of motion ST/ 6^=0 and ST/ 5D = for a given 
order in the (e.g. loop) expansion of r[^,D] (cf. Sec. 2). However, the 2PI effective 
action does not treat the higher n-point functions with n > 2 on the same footing as the 
lower ones: The three- and four-point function etc. are not self-consistently dressed in 
general, i.e. the corresponding proper three-vertex V3 and four-vertex V4 are given by 
the classical ones. In contrast, the nFl effective action r[^,D, V3, V4, . . . ,y„] provides 
a dressed description for the proper vertices Vs, V4, . . . , V,, as well, with dr/dV^ = 
0, 5r/5V4 = 0, . . . , 5r/5Vn = 0. 

The use of nPl effective actions with higher n > 2 is not entirely academic. They 
are relevant in the presence of initial-time sources describing a non-Gaussian initial 
density matrix for nonequilibrium evolutions (cf. the discussion in Sees. 3.2 and 3.3). 
They are also known to be relevant in high-temperature gauge theories for a quantitative 
description of transport coefficients in the context of kinetic theory. As an example, the 
calculation of shear viscosity in a theory like QCD can be based on the inclusion of an 
infinite series of 2PI "ladder" diagrams in order to recover the leading order "on-shell" 
results in the gauge coupling g. Further examples where approximation schemes based 



G will denote the ghost propagator in gauge theories below. 
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on higher nPI effective actions are relevant include critical phenomena near second-order 
phase transitions. For instance, the quantitative description of the universal behavior 
near the second-order phase transition of scalar 0^ theory goes beyond a 2PI loop 
expansion.^' It requires taking into account vertex corrections that start with the 4PI 
effective action to four-loop order. The latter agrees with the most general nFl loop 
expansion to that order, which is shown below. 

The evolution equations, which are obtained by variation of the nPI effective action, 
are closely related to known exact identities for correlation functions, i.e. Schwinger- 
Dyson (SD) equations. Without approximations the equations of motion obtained from 
an exact nPl effective action and the exact SD equations have to agree since one can 
always map identities onto each other. However, in general this is no longer the case for 
a given order in the loop or coupling expansion of the nPI effective action. By construc- 
tion, SD equations are expressed in terms of loop diagrams including both classical and 
dressed vertices, which leads to ambiguities of whether classical or dressed ones should 
be employed at a given truncation order. In particular, SD equations are not closed a 
priori in the sense that the equation for a given n-point function always involves infor- 
mation about m-point functions with m > n. These problems are absent using effective 
action techniques. In turn, the nPI results can be used to resolve ambiguities of whether 
classical or dressed vertices should be employed for a given truncation of a SD equation. 
For instance, in QCD the three-loop effective action result leads to evolution equations, 
which are equivalent to the SD equation for the two-point function and the one-loop 
three-point function if all vertices in loop-diagrams for the latter are replaced by the full 
vertices at that order. As mentioned in previous sections, the "conserving" property of 
using an effective action truncation can have important advantages, in particular if ap- 
plied to nonequilibrium time evolution problems, where the presence of basic constants 
of motion such as energy conservation is crucial. 

We will derive below the 4PI effective action for a nonabelian SU (N) gauge theory 
with fermions up to four- loop or ^{g^) corrections, starting from the 2PI effective action 
and doing subsequent Legendre transforms (Sees. 6.1 and 6.2). The class of models 
include gauge theories such as QCD or abelian theories as QED, as well as simple scalar 
field theories with cubic or quartic interactions. In Sec. 6.1.2 we derive an equivalence 
hierarchy for nPl effective actions, which implies that the 4PI results to this order are 
equivalent to those from the nPl effective action up to four-loop or ^{g^) corrections 
for arbitrary n > 4. We derive the non-equilibrium gauge field and fermion evolution 
equations (Sec. 6.2.3), and discuss the connection to kinetic theory in Sec. 6.3. 



^ Critical phenomena can be described using the 1 /N expansion of the 2PI effective action beyond 
leading order. These approximations indeed resum an infinite number of loop diagrams and, for instance, 
the NLO result can be rewritten as a "loop approximation in the presence of an effective four-vertex" 
(cf. Sees. 2.1 and 4.3). 
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6.1. Higher effective actions 



Recall that all information about the quantum theory can be obtained from the ef- 
fective action, which is a generating functional for Green's functions evaluated in the 
absence of external sources, i.e. at the physical or stationary point. All functional rep- 
resentations of the effective action are equivalent in the sense that they are generating 
functional for Green's functions including all quantum/statistical fluctuations and, in 
the absence of sources, have to agree by construction:^^ 

r[0] = r[0,D]=r[0,z),v^3]=r[(^),D,y3,V4] = ... = r[(|),D,y3,V4,...,K] (300) 

for arbitrary n without further approximations. However, e.g. loop expansions of the 
IPI effective action to a given order in the presence of the "background" field differ in 
general from a loop expansion of r[0,D] in the presence of ^ and D. A similar statement 
can be made for expansions of higher functional representations. As mentioned in 
the introduction in Sec. 1.1.2, for applications it is often desirable to obtain a self- 
consistently complete description, which to a given order in the expansion determines 
r[^,D, V3, V4, . . . , V„] for arbitrarily high n. For this it is important to realize that there 
exists an equivalence hierarchy as displayed in the introduction in Eq. (13), which is 
derived in Sec. 6.1.2. For instance at three-loop order one has: 

r(3io°p)[<^] ^r(3i°°p)[,^,D] 7^r(3i°°p)[0,D,y3] = r(3i°°p)[,^,D, 73,^4] (30i) 

= r(3i-p)[(^,D,y3,V4,...,K], 

for arbitrary n > 4 in the absence of sources. As a consequence, there is no differ- 
ence between T^'^^°°^^[^,D,Vi] and r(3'°°P)[0,D,V3,V4] etc. at the stationary point, and 
the 3PI effective action captures already the complete answer for the self-consistent 
description to this order. More explicitly, the equality reads r(-^^°°P)[0,D(0),V3(0)] = 
r(^^°°P) [^,D((^), V3(^), V4(^)] = . . . since at the stationary point of the effective action 
all n-point correlations become functions of the field expectation value 0. At four- loop 
order the 4PI effective action would become relevant etc. For instance, for a theory as 
quantum electrodynamics (QED) or chromodynamics (QCD) the 2PI effective action 
provides a self-consistently complete description to two-loop order or^^ ^{g^)'- For a 
two-loop approximation all nPl descriptions with n > 2 are equivalent in the absence of 
sources. In contrast, a self-consistently complete result to three-loop order or ^{g^) re- 
quires at least the 3PI effective action etc. To go to much higher loop-order can become 
somewhat academic from the point of view of calculational feasibility. 

To present the argument we will first consider the 4PI effective action for a simple 
generic scalar model with cubic and quartic interactions. The formal generalization to 



^ Of course, the nonequilibrium nPI effective action in the presence of initial-time sources representing 
an initial density matrix can differ in general. However, these pose no additional complications since they 
vanish identically for times different than the initial time. Cf. the discussions in Sees. 3.3 and 3.4.3. 

Here, and throughout the paper, g means the strong gauge coupling for QCD, while it should be 
understood as the electric charge e for QED. For the power counting we take ~ ^{1/g) (cf. Sec. 6.1.1). 
The metric is denoted as g^^ = g^v = diag(l, — 1, — 1, — 1). 
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fermionic and gauge fields is straightforward, and in Sec. 6.2 the construction is done 
for SU{N) gauge theories with fermions. We use here a concise notation where Latin 
indices represent all field attributes, numbering real field components and their internal 
as well as space-time labels, and sum/integration over repeated indices is implied. We 
consider the classical action 

1 , P p2 

S[(p] = :^(Pi iD^jj (Pj - —Vo3 J jk(pi(pj(Pk - -^Vo4jjki(Pi(pj(pk(pi , (302) 

where we scaled out a coupling constant g for later convenience. The generating func- 
tional for Green's functions in the presence of quadratic, cubic and quartic source terms 
is given by: 

Z[7,i?,i?3,i?4] = exp(/W[7,i?,i?3,^4]) 

^ J ^(p exp I / (^S[(p] + J, (pi -f- ]^R,j (Pi(pj (303) 

+ J^R3,ijk (Pi(Pj(Pk + ^R4,ijkl (Pi(Pj(Pk(Pl 

The generating functional for connected Green's functions, W, can be used to define the 
connected two-point (D), three-point (Dj) and four-point function (D4) in the presence 
of the sources. 



5W 



<^,, (304) 



l^{D,j + ^4j), (305) 



5Rij 2 



^ {D3,jk + D,j + ^fa- <^>,/ + Djk (pi + (j),(l)j(l)k) , (306) 



SR3,ijk 6 

5W 1 



^ = TTT {D4,ijki + [D3,ijk <^>/ + 3 perm.] + [DjjDki + 2 perm.] 

OKAdjkl ^4 

+ [D,j ^k^i + 5 perm.] + ^i^jMi) ■ (307) 



We denote the proper three-point and four-point vertices by gV3 and g^V4, respectively, 
and define 

D3,ijk = -igDii>DjfDkk>V3/fk' , (308) 
D4,ijki = -ig DiiiDjfDj,;,iDiiiV4j>fk>i> 

+DiiiDjiuiD]^ikDjyiDiii)V3^iijiuiV3^u'v'w' ■ (309) 
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.40 



The effective action is obtained as the Legendre transform of W[7,i?,i?3,i?4]: 

5W 5W 

5W 5W 

R3,ijk — 0-5 R4,ijkl ■ (310) 



,ijkl 

For vanishing sources one observes from (310) the stationarity conditions 

5r 5r dr dr 



0, (311) 



5(1) 8D 5V3 5V4 
which provide the equations of motion for ^, D, V3 and V4. 

6.7.7. 4PI effective action up to four-loop order corrections 

Since the Legendre transforms employed in (310) can be equally performed subse- 
quently, a most convenient computation of r[0,D, V3, V4] starts from the 2PI effective 
action r[0,Z)]. According to (28) the exact 2PI effective action can be written as: 

r[0,D] = 5[0] + ^TrlnD-i + ^TrD^ i((^))D + r2[(^,Z)] +const, (312) 
with the field-dependent inverse classical propagator 

To simplify the presentation, we use in the following a symbolic notation which sup- 
presses indices and summation or integration symbols (suitably regularized). In this no- 



In terms of the standard one-particle irreducible effective action T[<j)] ~W[J]— J(j) the proper vertices 
V3 and V4 are given by 

^ ^ 505050 ' ^ * 50505050 ■ 
Here it is useful to note that in terms of the connected Green's functions D,, one has 

S'W[J] ^ 5^ ^ , 

5757 ' 5050 ' 



575757 505050 ' 

57575757 ^ 50505050 \5(j)5(j)5<p J 
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tation the inverse classical propagator reads 

/Do = iDQ'-g(l)Vo3 - ^/(^)Vo4, (314) 
and to three-loop order one has"^^ (cf. Sec. 2.1) 

+lg^D\gVo3+g^(l>Vo4fVo4 - ^D\gVo3 +g^(l)Vo4)^ 

+ ^{g"{g^mn+.n^6). (315) 

for n, m = 0, . . . , 6. We emphasize that the exact ^-dependence of r2 , D] can be written 
as a function of the combination {gVo^ +^^<^>Vb4). In order to obtain the vertex 2PI 
effective action r[(j),D,V3,V4.] from r[(j),D], one can exploit that the cubic and quartic 
source terms ~ and ~ R4 appearing in (303) can be conveniently combined with the 
vertices gVo3 and g^Vo4 by the replacement: 

gVo3 ^ gVo3-R3 = gV3 , g^Vo4 ^ g\4 ' R4 = g'^ ■ (316) 



The 2PI effective action with the modified interaction is given by 

8W dW 

Since 



ry[^,D]=W[J,R,R3.R4]-^J- -^R . (317) 



5r,7 5W 5r,7 5W 



dR3 dR3 ' dR4 dR4 ' ^^^^^ 
one can express the remaining Legendre transforms, leading to r[^,D, V3,V4], in terms 
of the vertices V3, V4 and V03, V04: 

= r,|*,D]-^I||^(V3-V„3)-^^||^(V4-V„.). (319) 

What remains to be done is expressing V3 and V4 in terms of V3 and V4. On the one hand, 
from (307) and the definitions (308) and (309) one has 

= -U-igD'V3 + 3D<^ + <^'). (320) 
goV3 6 

= -^f-ig2D'V4-3g^D'vi-4igD'V3^ + 3D^ 
g'^oV4 24 V 

+60^ + ^"^). (321) 



Note that for ^ 7^ 0, in the phase with spontaneous symmetry breaking, (j) ^ i^{l/g), and the three-loop 
result (315) takes into account the contributions up to order g^. 
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On the other hand, from the expansion of the 2PI effective action to three-loop order 
with (315) one finds"^^ 

4 

+ <7-(j"(s"c«"|„+,„=5), (322) 

+ Id^{8V3 +g^(l)V4f - ^D^(l){gV3+g^(j)V4f 

+ ^{g"-H8^<l>r\n+m=6). (323) 

Comparing (322) and (320) yields 

8V3 = {gV3+g^^V4)-yg^D^{gV3+g^^V4)V4-D^{gV3+g^^V4)^ 

+ ^{g^{g^mn+m=5) . (324) 
Similarly, for V4 comparing (323) and (321), and using (324) one finds 

g^V4 = g^V4 + G {f-^{g^^TU„r^^) . (325) 
This can be used to invert the above relations as 

gV3+g^(pV4 = gV3 + yg^D^V3V4 + g^D^Vi + G{g^), (326) 

g^V4 = g^V4 + G{g^). (327) 

Plugging this into (319) and expressing the free, the one-loop and the r2 parts in terms 
of V3 and V4 as well as Vo3 and V04, one obtains from a straightforward calculation: 

r[(l),D,V3,V4]^S[(l)] + ^Tv\nD-' + ^TvDQ\(j))D + r2[(l>,D,V3,V4], (328) 

with 

r2[(l),D,V3,V4] = r^_[(l),D,V3,V4]+rf[D,V3,V4], (329) 



Note that since the exact -dependence of r2[0,Z)] can be written as a function of (gVo3 +^^0Vo4), the 
parametrical dependence of the higher order terms in the variation of (315) with respect to (gVoa) is given 
by £^i8"{g^m+>n=5) (cf. (322)). 
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rO[0,D,y3,V4] = -lg^D^Vo4 + ^gD^V3{gVo3+g^(l>Vo4) 

+^/£>V4yo4 + Ig'^D^ViVo^ , (330) 

rf[AV3,\/4] = '{^g'D\'-^g'D'vi-Y/D%' + ^{g'). (331) 

The diagrammatic representation of these results is given in Figs. 19 and 21 of Sec. 6.2.1. 
There the equivalent calculation is done for a SU (N) gauge theory and one has to replace 
the propagator lines and vertices of the figures by the corresponding scalar propagator 
and vertices. Note that for the scalar theory the thick circles represent the dressed 
three-vertex gV^, and four-vertex g^V4, respectively, while the small circles denote the 
corresponding effective classical three-vertex gVbs + g^^VoA and classical four-vertex 
g^Vo4. As a consequence, the diagrams look the same in the absence of spontaneous 
symmetry breaking, indicated by a vanishing field expectation value 0. 

In (328), the action S[^] and Dq depend on the classical vertices as before. The 
expression for Fj, which includes all terms of F2 that depend on the classical vertices, 
is valid to all orders: Tf^ contains no explicit dependence on the field or the classical 
vertices Vqs and Vq4, independent of the approximation for the 4PI effective action. This 
can be straightforwardly observed from (319), where the complete (linear) dependence 
of F on Vo3 and V04 is explicit, together with (320) and (321). 



6.1.2. Equivalence hierarchy for nPI effective actions 

As pointed out in Sec. 1.1.2 of the introduction, for applications it is often desirable 
to obtain a self-consistently complete description, which to a given order of a loop 
or coupling expansion determines the nPI effective action F[^,Z), V3, V4, . . . , V„] for 
arbitrarily high n. Despite the complexity of a general nPl effective action such a 
description can be obtained in practice because of the equivalence hierarchy displayed 
in Eq. (13): Typically the 2PI, 3PI or maybe the 4PI effective action captures already 
the complete answer for the self-consistent description to the desired/computationally 
feasible order of approximation. Higher effective actions, which are relevant beyond 
four-loop order, may not be entirely irrelevant in the presence of sources describing 
complicated initial conditions for nonequilibrium evolutions. However, their discussion 
would be somewhat academic from the point of view of calculational feasibility and we 
will concentrate on up to four-loop corrections or ff{g^) in the following. Below we will 
not explicitly write in addition to the loop-order the corresponding order of the coupling 
g for the considered theory, which is detailed above in Sec. 6.1.1. 

To show (13) we will first observe that to one-loop order all nPI effective actions 
agree in the absence of sources. The one-loop result for the IPI effective action is given 
by (18) for vanishing source R. As has been explicitly shown in Sec. 2 (cf. Eq. (26)), the 
one-loop 2PI effective action agrees with that expression, i.e. 

r(iioop)[,^,z)] =r("°°p)[(^], (332) 



Introduction to Nonequilibrium Quantum Field Theory 



107 



in the absence of sources. The equivalence with the one-loop 3PI and 4PI effective 
actions can be explicitly observed from the results of Sec. 6.1.1. In order to obtain 
the 3PI expressions we could directly set the source R4 = from the beginning in the 
computation of that section such that there is no dependence on V4. Equivalently, we 
can note from Eqs. (328)-(331) that already the 4PI effective action to this order simply 
agrees with (18) for zero sources. As a consequence, it carries no dependence on V3 and 
V4, i.e. 

r(iioop)[0,z),y3,y4] =r("°°p)[(^),AV3] =r(i'°°p)[0,D]. (333) 

For the one-loop case it remains to be shown that in addition 

r(iioop) [0,Z), Vg, V4, . . . , V„] = r(i'°°P) [(j>,D, V3, V4] (334) 

for arbitrary n > 5. For this we note that the number / of internal lines in a given loop 
diagram is given by the number V3 of proper 3-vertices, the number V4 of proper 4- 
vertices, . . . , the number v„ of proper n- vertices in terms of the standard relation: 

2/ = 3v3 + 4v4 + 5v5...-Fnv„, (335) 

where V3 + V5 + V7 + . . . has to be even. Similarly, the number L of loops in such a 
diagram is 

L = / — V3 — V4 — V5 . . . — v„ -t- 1 
1 3 /7 2 

= :^V3+V4+-V5... + —^Vn+l. (336) 

The equivalence (334) follows from the fact that for L = 1 equation (336) implies that 
r(iioop) [<^,D, V3,V4, Vn] cannot depend in particular on V5,... Vn.^^ 

The two-loop equivalence of the 2PI and higher effective actions follows along the 
same lines. According to (328)-(331) the 4PI effective action to two-loop order is given 
by: 

r(2i°°P)[,^,D,y3,V4] = 5[0] + ^TrlnD-i + ^TrDoi(0)D 

-frf°°P)[0,D,y3,V4], (337) 

rf°°p)[(^),AV3,y4] = -lg^D^Vo4 + yD^V3{gVo3+g^^Vo4)-^g^D^vl 

There is no dependence on y4 to this order and, following the discussion above, there is 
no dependence on ys, . . . , y„ according to (336) for L = 2. Consequently, 

r(2i-p)[0,z),y3,y4,...,y„] =r(2i-p)[0,D,y3,y4] =r(2i-p)[0,D,y^ (338) 



Note that we consider here theories where there is no classical 5-vertex or higher, whose presence 
would lead to a trivial dependence for the classical action and propagator 
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for arbitrary n in the absence of sources. The latter yields 



^ 5V3 ^ .?^03 + g 0Vb4 , (339) 

which can be used in (337) to show in addition the equivalence of the 3PI and 2PI 
effective actions (cf. Eq. (315)) to this order: 



rf°°P)[(^),D,y3] = -lg^D^Vo4 + ^D\gVo3 + g^<^Vo4f 

^(21oop) 



= r™[0,D], (340) 

for vanishing sources. The inequivalence of the 2PI with the IPI effective action to this 
order, 

r(2i°°p)[0,D] 7^r(2i°°p)[0], (341) 

follows from using the result of 5rf [(j),D]/5D = for D in (340) in a straightfor- 
ward way."*"^ 

In order to show the three-loop equivalence of the 3PI and higher effective actions, 
we first note from (328)-(331) that the 4PI effective action to this order yields V4 = V04 
in the absence of sources: 

Constructing the 3PI effective action to three-loop would mean to do the same calcu- 
lation as in Sec. 6.1.1 but with V4 V04 from the beginning {R4 = 0). The result of a 
classical four-vertex for the 4PI effective action to this order, therefore, directly implies: 

r(3i-p)[<^,D,y3,V4] =r(3i°°p)[(^,D,y3], (343) 

for vanishing sources. To see the equivalence with a 5PI effective action 
p(3ioop) y^j^ wc uotc that to thrcc-loop order the only possible diagram 

including a five- vertex requires V3 = V5 = 1 for L = 3 in Eq. (336). As a consequence, 
to this order the five-vertex corresponds to the classical one, which is identically zero 
for the theories considered here, i.e. V5 = Vqs = 0. In order to obtain that (to this order 
trivial) result along the lines of Sec. 6.1.1, one can formally include a classical five- 
vertex Vqs and observe that the three-loop 2PI effective action admits a term ~ D^VqsV^. 
After performing the additional Legendre transform the result then follows from setting 
^05 — > in the end. The equivalence with nFl effective actions for n > 6 can again be 
observed from the fact that for L = 3 Eq. (336) implies no dependence on yg, . . .y„. In 
addition to (343), we therefore have for arbitrary n>5: 

r(31oop) [^^ y^^ y^^ _ _ _ ^ y^] ^ r^^'-P) y3, y4] . (344) 



Here r'^'°°P^[(/),D] includes e.g. the summation of an infinite series of so-called "bubble" diagrams, 
which form the basis of mean-field or Hartree-type approximations, and clearly go beyond a perturbative 
two-loop approximation r(^'°°P)[0] (cf. Sec. 2.1). 
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The inequivalence of the three-loop 3PI and 2PI effective actions can be readily observed 
from (328)-(331) and (343): 



^p(31oop) \A Dy^\ . o , 

^ir^^^O ^ gV3 = g{ym+g^yQA)-g^D^Vl (345) 

Written iteratively, the above self-consistent equation for V3 sums an infinite number of 
contributions in terms of the classical vertices. As a consequence, the three-loop 3PI 
result can be written as an infinite series of diagrams for the corresponding 2PI effective 
action, which clearly goes beyond r(^'°°P^ (cf. Eq. (315)): 

r(31oop) ^ ^ y^] ^ p(31oop) ^ ^] _ (34g) 

The importance of such an infinite summation will be discussed for the case of gauge 
theories below. 



6.2. Nonabelian gauge theory with fermions 

We consider a SU (N) gauge theory with Nf flavors of Dirac fermions with classical 
action 

•^eff = S + Sgf + 5fpg 
- / -VA(-/^)V/-f]«<9^(D^T])") , (347) 

where i/a (y), A and rj (fj) denote the (anti-)fermions, gauge and (anti-)ghost fields, 
respectively, with gauge-fixing term ^"(A) = d^A"^ for covariant gauges. The color 

indices in the adjoint representation are a,b,... = l,...,N^ — 1, while those for the 
fundamental representation will be denoted by j, . • • and run from 1 to N. Here 

F;, = ^^A^-^vA^^-^r^^A'^A^,, (348) 
(D^T])" = d^'n"-gf"'''A^'''q'\ (349) 
ip = f(^d^+igAy^, (350) 

where [?^?*] = tr{t"t'') = 5"^ /2. For QCD, f = jl with the GeU-Mann 

matrices {a= 1, . . . , 8). We will suppress Dirac and flavor indices in the following. It 
is convenient to write (347) as 

^eff = UA^"{x)iD-l^f{x,y)A^'{y)+[ r{x)iG,'^'{x,y)v'iy) 

Z Jxy Jx}' 

+ I %{^)i^o!M.y)Wj{y)-\gl VS^';,,^{x,y,z)A^^{x)A'\y)Ay^{z) 

Jxy D Jxyz 

- I <'?vr5(^.3'.^.>^M^'^W^''(3'M'^'«A^^(w) 

Z^ Jxyzw 
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Jxyz 

- s[ V^'^; {x,y;zmx)Vj{y)A^%z). (351) 

Jxyz 

with the free inverse fermion, ghost and gluon propagator in covariant gauges given by 

i\!ji^^y) = i^^ij^ix-y), (352) 
iGo'"'{x,y) = -a5«*S^(^-j), (353) 
iDo'^fM = [g^vD-{\-^-')d^dy]J"'^d<^{x-y), (354) 

where we have taken the fermions to be massless. The tree-level vertices read in coordi- 
nate space: 

Vo%%ri^.y,z) = 

giiv[d^Ay-z)dy5^{x-y)-5^{x-z)d^y&^{y-x)] 
+ g^iy[^{x-y)dydc^{z-x)-5<(^{y-z)d^5^{x-z)] 
+ gvy[^{x-z)dl5^{y-x)-&g{x-y)di,5^{z-x)]^, 

(355) 

+ r''f''''[g^vgy8-g^^8gvr] +r'''f%^rS8v-g^vgr8]) 
d^{x-y)5<^{x-z)&^{x-w), (356) 

VolT^'i^^y^z) = -r''d-;d^ix-z)d^iy-z). (357) 
yiZij^^^y^^^ = r^itfM^-z)&rf{z-y). (358) 

Note that VQ^'Ibci^^y^z) is symmetric under exchange of {ji,a,x) ^ {v,b,y) ^ (7,c,z). 

Likewise, ^04^afocd(^'^'^'^) symmetric in its space-time arguments and under ex- 
change of ill, a) ^ (v,Z?) ^ (7,c) ^ (5, J). 

In addition to the linear and bilinear source terms, which are required for a construc- 
tion of the 2PI effective action, following Sec. 6.1 we add cubic and quartic source terms 
to (351): 

^'source = U Rf^yy{x.y,z)Af'''{x)A''iy)Ay^{z) 
D Jxyz 

+ i / <;?r5(^'3''^'^M^"W^''(3'M'''(^M^^(w) 

^'T' Jxyzw 

+ [ R[^l;^"'^'\x,y-,z)ri^)r]'iy)A^'\z) 

Jxyz 
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+ / Rf;,jix,y;z)n^)¥j{y)A^''{z). (359) 

Jxyz 

where the sources R2 4 obey the same symmetry properties as the corresponding classical 
vertices Vbs and V04 discussed above. The definition of the corresponding three- and four- 
vertices follows Sec. 6.1. In particular, we have for the vertices involving Grassmann 
fields: 

5W 



-ig [ D^'^'-'(z,zOA,,(j,y)y3(|],'^ .,(y,y;z')A/,(y,..) ,(360) 

Jx'y'z' 



for the case of vanishing "background" fields (A) = (i//) = (i//) = (r]) = (r]) = 0, which 
we will consider in the following. 



6.2.7. Effective action up to four-loop or &{g ) corrections 

Consider first the 2PI effective action with vanishing "background" fields, which 
according to Sec. 2 can be written as 

r[D, A, G] = ^TrlnD^i -t- ^TrD^ - zTrlnA^^ - zTr Aq 

-zTrlnG"^ -zTrGo^G-Fr2[AA,G]. (361) 

Here the trace Tr includes an integration over the time path ^, as well as integration 
over spatial coordinates and summation over flavor, color and Dirac indices. The exact 
expression for Y2 contains all 2PI diagrams with vertices described by (355)-(358) and 
propagator lines associated to the full connected two-point functions D, G and A. In 
order to clear up the presentation, we will give all diagrams including gauge and ghost 
propagators only. The fermion diagrams can simply be obtained from the corresponding 
ghost ones, since they have the same signs and prefactors."^^ For the 2PI effective action 
of the gluon-ghost system, r[D,G], to three-loop order the 2PI effective action is given 
by (using the same compact notation as introduced in Sec. 6.1.1): 

+\g'D'vlV,, - ^/dX, + y'D^G\^f\3 
-K^/D2GV(fV^(/). (362) 



Note that to three-loop order there are no graphs with more than one closed ghost/fermion loop, such 
that ghosts and fermions cannot appear in the same diagram simultaneously. 
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FIGURE 19. The figure shows together with Fig. 20 the diagrammatic representation of 

r2[Z),G,V3,y3^''\V4] as given in Eq. (365). Here the wiggled lines denote the gauge field propagator 
D and the unwiggled lines the ghost propagator G. The thick circles denote the dressed and the small ones 
the classical vertices. This functional contains all terms of that depend on the classical vertices gVoT,, 

gV^f^^ and g^VQ4 for an SU{N) gauge theory. There are no further contributions to Fj appearing at higher 
order in the expansion. For the gauge theory with fermions there is in addition the same contribution as in 
Fig. 20 with the unwiggled propagator lines representing the fermion propagator A and the ghost vertices 
replaced by the corresponding fermion vertices Vqj and Vj (cf. Eq. (358)). 



- 1 



FIGURE 20. Ghost/fermion part of F^ . 



The result can be compared with (315) and taking into account an additional factor 
of (—1) for each closed loop involving Grassmann fields (cf. Sec. 2.3). Here we have 
suppressed in the notation the dependence of r2[£>, G] on the higher sources (359). The 
desired effective action is obtained by performing the remaining Legendre transforms: 

riz,.c,V3.vr.v.l^riz,.G]-||.3-^<'-|^... (363, 

The calculation follows the same steps as detailed in Sec. 6.1.1. For the effective action 
to ^{g^) we obtain: 

r[D,G,V3,V^^'^\v4\ = ^Tr\nD-^ + ^TrDQ^D-iTr\nG-^-iTrGQ^G 

+r2[AG,y3,V3^''\V4], (364) 



with 



r2[AG,y3,V3^''\v4] = r«[D,G,y3,\/3(«*'\y4]+rf[AG,y3,vi''V4], 

li[D,G,V,,vi^\v,] = -l8'D^Vo, + y'D\,Vo3-ig'DGWi'%^f 

+:^g^D^V4Vo4 + ^g^D^ViVo4 , (365) 

rf[D,G,y3,vf\v4] = -{^s'D'vi + y^DG^v^^''^'-^g'DX 



+Lg^D^G'vj^''^' + ff{g'). (366) 
4 
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FIGURE 21. The figure shows together with Fig. 22 the diagrammatic representation of 

rf^[D,G,V3,V^^^\V4\ to three-loop order as given in Eq. (366). For the gauge theory with fermions, 
to this order there is in addition the same contribution as in Fig. 22 with the unwiggled propagator lines 
representing the fermion propagator A and the ghost vertex replaced by the corresponding fermion vertex 

(f) 

V3 . This functional contains no explicit dependence on the classical vertices independent of the order of 
approximation. 

FIGURE 22. Ghost/fermion part of to three-loop order 



The contributions are displayed diagrammatically in Figs. 19 and 20 for F^, and in 
Figs. 21 and 22forrf . 

The equivalence of the 4PI effective action to three-loop order with the 3PI and nPI 
effective actions for « > 5 in the absence of sources follows along the lines of Sec. 6. 1 .2. 
As a consequence, to three-loop order the nPl effective action does not depend on higher 
vertices Vs, Ve, . . . V„. In particular with vanishing sources the four-vertex is given by the 
classical one: 

gr(3ioop)[D,G,y3,vf\v4] grf°°p)[D,G,V3,vi^^\v4] ^ ,,,,, 

dV4 ^ W, = 0^y4 = Vo4. (367) 

If one plugs this into (365) and (366) one obtains the three-loop 3PI effective action, 
p(3loop) [0 , £), V3 , V^^^^ ] . Similarly, to two-loop order one has 



5r(2i°op)[D,G, 73,^3^''^] _ 5rf°°P)[D,G,V3,V3^ 



(gh)i 



5V3 5V3 



^ = V03 : 



(f) 

and equivalently for the fermion vertex V3 . To this order, therefore, the combinatorial 
factors of the two-loop diagrams of Fig. 19 and 21 for the gauge part, as well as of 
Fig. 20 and 22 for the ghost/fermion part, combine to give the result (362) to two-loop 
order for the 2PI effective action. 

We have seen above that to two-loop order the proper vertices of the nPl effective 
action correspond to the classical ones. Accordingly, at this order the only non-trivial 
equations of motion in the absence of "background" fields are those for the two-point 
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FIGURE 23. The self-energy for the gauge field (H) and the ghost/fermion (E) propagators as obtained 
from the self-consistently complete two-loop approximation of the effective action. Note that at this order 
all vertices correspond to the classical ones. 



functions: 

dr dr dr 

55-°' SG=°' SA-°' <^**' 
for vanishing sources. Applied to an nPl effective action (n> 1), as e.g. (364), one finds 
for the gauge field propagator: 

D i=Doi-n , n = 2ij^. (369) 
The ghost propagator and self-energy are 

G-' = Gq'-L , ^=-i^, (370) 

and equivalently for the fermion propagator A. (Cf. also Sec. 2 for the same relations 
in the context of 2PI effective actions.) The self-energies to this order are shown in 
diagrammatic form in Fig. 23. In contrast, for the three-loop effective action the three- 
vertices get dressed and the stationarity conditions, 

5r 5r 5r 

0' ^::^ = ^^ ^?v = 0' (371) 



8V3 ' 5y^(gh) ' gyif) 

applied to (364)-(366) lead to the equations shown in the left graph of Fig. 25. Here 
the diagrammatic form of the contributions is always the same for the ghost and for the 
fermion propagators or vertices. We therefore only give the expressions for the gauge- 
ghost system. If fermions are present, the respective diagrams have to be added in a 
straightforward way. 

The self-energies to this order are displayed in Fig. 24. It should be emphasized that 
their relatively simple form is a consequence of the equations for the proper vertices. 
Fig. 25. To see this we consider first the many terms generated by the functional 
derivative of (365) and (366) with respect to the gauge field propagator: 

^p(31oop) . ^ ^ 
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tt(3) ^ _i ^.J^^w^ -1 
^^2 2 



+ 







FIGURE 24. The self-energy for the gauge field (11) and the ghost/fermion (E) propagators as obtained 
from the self-consistently complete three-loop approximation of the effective action. (Cf. Fig. 25 for the 
vertices.) 



The short form for the self-energy of Fig. 24 is obtained through cancellations by 
replacing in the above expression 

(373) 

as well as 

The latter equations follow from inserting the expressions for the dressed vertices of 
Fig. 25. Noting in addition that the proper four- vertex to this order corresponds to 
the classical one (cf. (367)) leads to the result. Along the very same lines a similar 
cancellation yields the compact form of the ghost/fermion self-energy displayed in 
Fig. 24. 



6.2.2. Comparison with Schwinger-Dyson equations 

The equations of motions of the last section are self-consistently complete to two- 
loop/three-loop order of the nPI effective action for arbitrarily large n. We now compare 
them with Schwinger-Dyson (SD) equations, which represent identities between n-point 
functions. Clearly, without approximations the equations of motion obtained from an 
exact nPI effective action and the exact (SD) equations have to agree since one can 
always map identities onto each other. However, in general this is no longer the case for 
a given order in the loop expansion of the nPI effective action. 

By construction each diagram in a SD equation contains at least one classical vertex. 
In general, this is not the case for equations obtained from the nVl effective action: The 
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Gauge three-vertex as well as ghost/fermion 
vertex from r(3i°°P): 



i 





Compare: exact SD equation (ghost/fermion 
diagrams not displayed): 




FIGURE 25. Left: The gauge field three-vertex as well as the ghost (fermion) vertex as obtained from 
the self-consistently complete three-loop approximation of the effective action. Apart from the isolated 
classical three-vertex, all vertices in the equations correspond to dressed ones since at this order the four- 
vertex equals the classical vertex. The result reflects the proper symmetry of the three-vertex. Right: 
Schwinger-Dy son equation for the proper three-vertex V3 . Additional diagrams involving ghost or fermion 
vertices are not displayed for brevity. We show it for comparison with the three-loop effective action result 
displayed on the left. A naive truncation of the Schwinger-Dyson equation at the one-loop level does not 
agree with the latter and symmetries (the second and third diagram contain a classical three-vertex instead 
of a dressed one). 



loop contributions of F™' in Eq. (366) or Figs. 21-22 are solely expressed in terms of 
full vertices. However, to a given loop-order cancellations can occur for those diagrams 
in the equations of motion which do not contain a classical vertex. For the three-loop 
effective action result this has been demonstrated above for the two-point functions. 
Indeed, the equations for the two-point functions shown in Fig. 24 correspond to the SD 
equations, if one takes into account that to the considered order the four- vertex is trivial 
and given by the classical one (cf. Eq. (367)). However, such a correspondence is not 
true for the proper three- vertex to that order. 

As an example, we show on the right of Fig. 25 the SD equation for the proper 
three-vertex, where for brevity we do not display the additional diagrams coming from 
ghost/fermion degrees of freedom. One observes that a naive neglection of the two-loop 
contributions of that equation would not lead to the effective action result for the three- 
vertex shown in Fig. 25. Of course, the straightforward one-loop truncation of the SD 
equation would not even respect the property of V3 being completely symmetric in its 
space-time and group labels. This is the well-known problem of loop-expansions of SD 
equations, where one encounters the ambiguity of whether classical or dressed vertices 
should be employed at a given truncation order. 

We emphasize that these problems are absent using effective action techniques. The 
fact that all equations of motion are obtained from the same approximation of the ef- 
fective action puts stringent conditions on their form. More precisely, a self-consistently 
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complete approximation has the property that the order of differentiation of, say, r[D, V] 
with respect to the propagator D or the vertex V does not affect the equations of motion. 
Consider for instance: 



dr[D,v = v{D)] _ dr 

5D ~ 5D 



5r 

dV 



5V 

(375) 

D oD 



If y = V{D) is the result of the stationary condition ST/ 5V = then the above corre- 
sponds to the correct stationarity condition for the propagator for fixed V. ST/ 5D = 0. 
In contrast, with some ansatz V = f{D) that does not fulfill the stationarity condition 
of the effective action, the equation of motion for the propagator would receive addi- 
tional corrections ~ 5V / SD. In particular, it would be inconsistent to use the equation 
of motion for the propagator 5r/ 5D = (cf. e.g. Fig. 24 which corresponds to the SD 
equation result) but not the equation dT/ dV = for the vertex (cf. Fig. 25). 

In turn, one can conclude that a wide class of employed truncations of exact SD equa- 
tions cannot be obtained from the nFl effective action: this concerns those approxima- 
tions which use the exact SD equation for the propagator but make an ansatz for the 
vertices that differs from the one displayed in Fig. 25. The differences are, however, typ- 
ically higher order in the perturbative coupling expansion and there may be many cases, 
in particular in vacuum or thermal equilibrium, where some ansatz for the vertices is 
a very efficient way to proceed. Out of equilibrium however, as mentioned above, the 
"conserving" property of the effective action approximations can have important conse- 
quences, since the effective loss of initial conditions and the presence of basic constants 
of motion such as energy conservation is crucial. 



6.2.3. Nonequilibrium evolution equations 

Up to G'{g^) corrections in the self-consistently complete expansion of the effective 
action, the four-vertex parametrizing the diagrams of Figs. 24-25 corresponds to the 
classical vertex. At this order of approximation there is, therefore, no distinction between 
the coupling expansion of the 3PI and 4PI effective action. To discuss the relevant 
differences between the 2PI and 3PI expansions for time evolution problems, we will 
use the language of QED for simplicity, where no four-vertex appears. However, the 
evolution equations of this section can be straightforwardly transcribed to the nonabelian 
case by taking into account in addition to the equation for the gauge-fermion three- 
vertex those for the gauge-ghost and gauge three-vertex (cf. Fig. 25). In the following 
the effective action is a functional of the gauge field propagator Dyiy{x,y), the fermion 

propagator A(x,);) and the gauge-fermion vertex V-^^{x,y\z), where we suppress Dirac 

(f) 

indices and we will write V3 = V. According to Eqs. (364) — (366) one has in this case 

r2 [D, A, V] = [D, A, V] 4- rf [D, A, y] , (376) 

with 

ro = -ig" [ tr [r^Aix,y)Vy{y,z;u)A{z,x)D^''ix,u)] , (377) 

J xyzu 
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where the trace acts in Dirac space. For the given order of approximation there are two 
distinct contributions to Fj"': 

Ff = F^"^+Ff + ^(/) , (378) 



= tr[V^{x,y-z)A{y.u)Vy{u,v-w)AMD^\z.w 

^ J xyzuvw 

F^') = V/ ix[V^{x,y;z)A{y,u)Vy{u,v;w)A{v,^) 



xyzuvwxl y'zlu'v'w' 

Vp{^ ,y'-z!)A{y\uya{u\v'-w')A{v\x)D^P{z,^)D^''{w,W^ . 

The equations of motions for the propagators and vertex are obtained from the station- 
arity conditions (368) and (371) for the effective action. To convert (369) for the photon 
propagator into an equation which is more suitable for initial value problems, we convo- 
lute with D from the right and obtain for the considered case of vanishing "background" 
fields, e.g. for covariant gauges (cf. also the discussion in Sec. 3.4): 

[g%U-{\-^-')d^dy]Dy-{x,y)-ijTl%{x,z)D'^{z.y) 

= ig^''&^{x-y). (379) 

Similarly, the corresponding equation of (370) yields the evolution equation for the 
fermion propagator given in Eq. (141). Using the above results the self-energies are 

EW(x,y) = -g" f D^v(z',j)V^^,z";z')A(/,j)/, (380) 
■ 'z!z" 

Tl^\x,y) = f tr/A(x,zOV^(z',z";y)A(/,x). (381) 

J z'z" 

Note that the form of the self-energies is exact for known three-vertex. To see this 
within the current framework, we note that the self-energies can be expressed in terms 
of Fj only. The latter receives no further corrections at higher order in the expansion 
(cf. Sec. 6.2.1), and thus the expression is exactly known: With 

/.</.(...)A(.,)^-,-^(^.^)a(.,). (382) 

and since Fj"' is only a functional of VAD^/^ (cf. Sec. 6.2.1) one can use the identity 

^rf , _ /■ , . 5Ff 



f 01 f 

/ gA/^ x Mz,y) = / V^{x,z;z) , 
JzOA{z,x) Jzz! V^{y,z;z') 



51^ 
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to express everything in terms of the known"^^ T^. The last equality in (383) uses that 
5(r2 + Tf^) /5A = 0. A similar discussion can be done for the photon self-energy. As 
a consequence, all approximations are encoded in the equation for the vertex, which is 
obtained from (378) as 

Jvwx'y'u'w' 

A{y',u')V,y{u',y;w')D''''{w',w) + ff{g^) , (384) 

where 

VQ^{x,y;z)=f5{x-z)5{z-y). (385) 

For the self-consistently complete two-loop approximation the self-energies are given 
by 

L^f\x,y) = -g^D^y{x,y)fA{x,y)f + ^{g'), (386) 

n^\x,y) = ghrfA{x,y)fAiy.x) + ff{g^). (387) 

Following the discussion of Sec. 3.4.1, we decompose the two-point functions into spec- 
tral and statistical components using the identities (143) for gauge fields and (121) for 
fermions. Then po corresponds to the gauge field spectral function and Fjj is the statis- 
tical two-point function, while p^-^^ and F^-^^ are the corresponding fermion two-point 
functions. The same decomposition can be done for the corresponding self-energies: 

n^'ix^y) = nf;)(x,y)-^nf;)(x,j)sign(/-/), (388) 

e(/)(x,j) = E^^)(x,3;)-^Ej,-^)(x,3;)sign(/-/), (389) 

and similarly for the fermions as described by (126). Since the above decomposition for 
the propagators and self-energies makes the time-ordering explicit, we can evaluate the 
r.h.s. of (379) along the time contour following the discussion of Sec. 3.4.3. One finds 
the evolution equations: 

[g%D-{\-^-')d^dy\^pj;'{x,y) = j^^'dznjj;(x,z)pz,/(z,j) , (390) 
[g\D-{\-^-')d^dy\/^{x,y) = r\zn'^J{x,z)Fo/{z,y) 

Jto 

- / dznf/(x,z)pz),/(z,);) , (391) 

where we use the abbreviated notation Jl^ dz = Jl^ dz'^ d'^z- Here the initial time is 
denoted by to, which was taken without loss of generality to be zero in the respective 



° This can also be directly verified from (377) to the given order of approximation. 
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equations (138) for scalars. The equations of motion for the fermion spectral and statis- 
tical correlators are obtained in a similar way from (141) as described in Sec. 3.4.3 and 
are given in Eq. (142). 

A similar discussion as for the two-point functions can also be done for the higher 
correlation functions. For the three- vertex we write 

V^{x,y;z)=V^{x,y-z)+V^{x,y-z). (392) 

and the corresponding decomposition into spectral and statistical components reads 

U^P^{x,y;z) sign(/ -/) sign(z« -x«) - \u^pP,y;z) sign(/ -z°) 

+ Vl^^{x,y-z) sign(xO -zO) sign(/ -z^) _ '-Vl^^^ix^y-z) sign(/ -/) 

+ j;z) sign(zO -/) sign(xO -/) - j;z) sign(zO (393) 

To discuss this in more detail we use the short-hand notation 

0(xO,/,zO) =0(xO-/)0(/-zO). (394) 

With the separation of Eq. (392), the time-ordered three- vertex can be written as (cf. also 
the corresponding discussion for two-point functions in Sec. 3.4.1) 

V^{x,y-z) = y(';)(x,j;z)0(xO,/,z«)-t-\/(^)(x,3;;z)0(/,zO,/) 
+ y(^)(x,3;;z)0(z«,xO,/)-fy(^)(x,y;z)0(z«,/,xO) 

y(^)(x,);;z)0(/,z«,/)-Fy(^)(x,3;;z)0(/,/,zO), (395) 

with 'coefficients' Vl^.^{x,y;z), i= 1, . . . , 6. These coefficients can be expressed in terms 

of three spectral vertex functions U^^^ {x,y;z), Vj^^^ {x,y;z) and ^^p^ ju). as well as the 

corresponding statistical components U^p~^{x,y;z), y^^^(x,y;z) and W^p-^{x,y;z) that have 
been employed in Eq. (393). One finds, suppressing the space-time arguments: 

^1) - ^(F) + V) V) 2V (p) (p) (p)y' 
^2) - ^(F) V) + V) 2V (p) V)^^(p)y' 
^(3) - -^(')+^f)+^(n-K-^(p)+^('p)+^(p))' 

^(5) - ^r.)+^(')-^(')+K^(p)+c-^(p))' ^^^^^ 

^5) - ^(F) V)^V)^2V (p) ^p)^'^(p)y' 
y^ = -u^ +y^ +-( -u^ ) 
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In terms of the coefficients V^r^ these are given by: 

= Ifv^ +v^) =-(v^ -v^ -v^) 

^(F) 4l,^i) + ^2) + V)^^5); ' ^(p) 2V (1) (2) (4) ^5);' 

V) 4i,^i)+^3)+V)^>);' ^p) 2V (i)^ (3) V) ^6);' 

_ l/'yM , yM , yM , yM A _ 1 ("yM , yM _ yM _ yM A 

Insertion shows the equivalence of (395) and (393). 

6.3. Kinetic theory 

To make contact with frequent discussions in the literature, we will consider for the 
above gauge field and fermion nonequilibrium equations a standard "on-shell" approx- 
imation which is typically employed to derive kinetic equations for effective particle 
number densities. This part can be viewed as a continuation of Sec. 4.1.3, where scalar 
fields and the limitations of "on-shell" approximations have been discussed. 

6.3.1. "On-shell" approximations 

The evolution equations (390)-(142) to order and higher contain "off-shell" and 
"memory" effects due to their time integrals on the r.h.s. (cf. also Sec. 4.1). To simplify 
the description one conventionally considers a number of additional assumptions which 
finally lead to effective kinetic or Boltzmann-type descriptions for "on-shell" particle 
number distributions. The derivation of kinetic equations for the two-point functions 
F^'^{x,y) and p^'^{x,y) of Sec. 6.2.3 can be based on (i) the restriction that the initial 
condition for the time evolution problem is specified in the remote past, i.e. Iq 
(ii) a derivative expansion in the center variable X = {x+y) /2, and (iii) a "quasiparticle" 
picture. In contrast to the discussion for scalar fields in Sec. 4.1.3, within this approach 
one first sends the initial time Iq to the remote past in the equations (390)-(142). This 
allows one to use standard derivative expansion techniques in a straightforward way. The 
procedure of Sec. 4.1.3 has the advantage that one can discuss which contributions to 
the evolution are lost in this limit. On the other hand, the advantage of the derivative ex- 
pansion is that it may in principle be used to include higher order corrections. However, 
the complexity of a derivative expansion grows rapidly beyond the lowest order. 

For the sake of simplicity (not required), we consider the Feynman gauge ^ = I 
in the following. We will also consider a chirally symmetric theory, i.e. no vacuum 
fermion mass, along with parity and CP invariance. Therefore, the system is charge 
neutral and, in particular, the most general fermion two-point functions can be written in 
terms of vector components only: F^f\x,y) = 'Y^F^f^^{x,y), p'^f\x,y) — Y^p^-f^^^{x,y), 
with hermiticity properties F(-^^^(;c,y) = [F'^-f^^{y,x)]*, p^f^^{x,y) = -[p^f^^{y,x)]* .For 
the gauge fields the respective properties of the statistical and spectral correlators read 
F^'{x,y) = [F;^{y,x)r, p^\x,y) = -[p'/iy^x)]*. 
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In order to Fourier transform with respect to the relative coordinate = — , we 
write 

F^iX^) = |d4..*Tr(x + ^,X-0 , (397) 

p^'iX,k) = -/|d4.e*>r(^ + |.^-0 . (398) 

and equivalently for the fermion statistical and spectral function, F^f\X,k) and 
p^f\X,k). Here we have introduced a factor —i in the definition of the spec- 
tral function transform for convenience. For the Fourier transformed quantities 
we note the following hermiticity properties, for the gauge fields: [F^^(X,fc)]* = 
F^^{X,k) , [p^^{X,k)]* = p^^{X,k), and for the vector components of the fermion 
fields: [F^f'^^'{X,k)]* = F^f^^{X,k) , [p^f^^'iX.k)]* = p(-^)^(X,;t). After sending 
to — °° the derivative expansion can be efficiently applied to the exact Eqs. (390) — 
(142). Here one considers the difference of (390) and the one with interchanged 
coordinates x and y, and equivalently for the other equations. We use 

Jd^se'"' Jd\nx,z)g{z,y) = nx,k)giX,k) + ... (399) 

jd'se^'^' jd\jd'^f{x,z)g{z.^W,y) = f{X,k)g{X,k)h{X,k) + .. . 

where the dots indicate derivative terms, which will be neglected. E.g. the first derivative 
corrections to (399) can be written as a Poisson bracket, which is in particular important 
if "finite-width" effects of the spectral function are taken into account. However, a 
typical quasiparticle picture which employs a free-field or "zero-width" form of the 
spectral function is consistent with neglecting derivative terms in the scattering part. 
We also note that the quasiparticle/free-field form of the two-point functions implies 

F^\X,k)^-g^^FD{X,k) , p^J{X,k)^-g^'^pD{X,k). (400) 

At this point the only use of the above replacement is that all Lorentz contractions can be 
done. This doesn't affect the derivative expansion but keeps the notation simple. Similar 
to Eq. (398), we define the Lorentz contracted self-energies: 

-4%)(X,fc) ^ |d4.e'^-^nJ,)^(x-K^,X-0, (401) 
-4n(p)(X,fc) ^ -i j ASe'ks^l^^(^X+'-,X-':^i . (402) 

Without further assumptions, i.e. using the above notation and applying the approxima- 
tion (399) and (400) to the exact evolution equations one has"^^ 

2k^^FD{X,k) = n^p)iX,k)FoiX,k)-n(^p)iX,k)pDiX,k), (403) 



The relation to a more conventional form of the equations can be seen by writing: 

(n(p)/b-n(f)pD)(x,^) = 
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(404) 



One observes that the equations (403) and (404) have a structure reminiscent of that 
for the exact equations for vanishing "background" fields, (390) and (391), evaluated 
at equal times = y*^- However, one should keep in mind that (403) and (404) are, in 
particular, only valid for initial conditions specified in the remote past and neglecting 
gradients in the collision part. 

From (404) one observes that in this approximation the spectral function receives 
no contribution from scattering described by the r.h.s. of the exact equation (390). As a 
consequence, the spectral function obeys the free-field equations of motion. In particular. 



Pd i^^y) have to fulfill the equal-time commutation relations [p^ {x,y)]^Q^^{) = and 

-g^'^5{x — y) in Feynman gauge. The Wigner transformed free- 



o_,,o 



field solution solving (404) then reads po{X,k) = poik) = 2;rsign(fc^) d{k^). A very 
similar discussion can be done as well for the evolution equations (142) for fermions, 
which is massless due to chiral symmetry as stated above. Again, in lowest order in 
the derivative expansion the fermion spectral function obeys the free-field equations of 



motion and one has p^J>{X,k) = p^^>{k) = Inljt sign(fc") d{k^). 

Assuming a "generalized fluctuation-dissipation relation" or so-called "Kadanoff- 
Baym ansatz": 



FoiX.k) 
F^f\x,k) 



1 



+ nD{X,k) 



--n^f\x,k) 



pD{X,k), 

P^f\x.k). 



(405) 



one may extract the kinetic equations for the effective photon and fermion particle num- 
bers no and n^f\ respectively. Considering spatially homogeneous, isotropic systems for 
simplicity, we define the on-shell quasiparticle numbers {t = X'^) 



(406) 



and look for the evolution equation for n/)(?,k). Here it is useful to note the symmetry 
properties 

FD{t,-k) = Foit^k) , pD{t,-k) = -pD{t,k), 
F^f\t,-kY = ~F^^\t,kf rp^^\t,-kY = P^^\t,k)^. (407) 



Applied to the quasiparticle ansatz (405) these imply 

nD{t,-k) = -[nD{t,k) + l\ , n^^\t,-k) = - 



,(/), 



t,k)-\ 



(408) 









1 








1 




( 


) + -Ti(p) 




Fd-^Pd 




) - -n(p) 




Fd + ^Pd 


) 



{X,k). 



The difference of the two terms on the rh.s. can be directly interpreted as the difference of a so-called 
"loss" and a "gain" term in a Boltzmann-type description. 
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FIGURE 26. Infinite series of self-energy contributions with dressed propagator lines and classical 
vertices. 



This is employed to rewrite terms with negative values of k^. To order the self- 
energies read (cf. Eq. (387)): 



dV 



F^f\x,k + p)^F^-f\x,p)^ 



1 



p'^f\x,k + p)^'p'^^){X,p 



n(p)(x,fc) = 2g 



dv 



F(/)(Z,fc + p)^p(/)(X,p)^ 



-p(/)(X,fc + p)^F(/)(X,p)^ 
From the equations (403) and (405) one finds at this order: (q = k — p) 

dV 1 



^'"-^^''^)=^''V(2;r)3 2k2p2q 

(n(^)(?,p)/^)(f,q)[no(?,k) + l] 



n(^)(r,p)-l n(-^^(?,q)-l nz)(^k) )2;r5(k-p-q) 



,(/) 



+ 2 



+ 



n(^)(r,p)-l n^^^(r,q) [nz,(r, k) + l] 



,(/) 



n(^)(r,p) n^-^^(f,q)-l nz)(r,k) )2;r5(k + p - q) 



,(/), 



?,q)-l 



[/iz)(r,k) + l] 



-r(^) (r, p) q) nD(/', k)) 2;r5(k + p + q 



(409) 



(410) 



The r.h.s. shows the standard "gain term" minus "loss term" structure. E.g. for the case 
k^ > 0, k*^ > the interpretation is given by the elementary processes ee — > 7, e ^ ey, 
e ^ ey and "0" — > eey from which only the first one is not kinematically forbidden. 
From (410) one also recovers the fact that the "on-shell" evolution with k^ = vanishes 
identically at this order. A nonvanishing result is obtained if one takes into account 
"off-shell" corrections for a fermion line in the loop of the self-energy (409). As a 
consequence the first nonzero contribution to the self-energy starts at &'{g'^). 

Since the lowest order contribution to the kinetic equation is of ff{g^), the 3PI 
effective action provides a self-consistently complete starting point for its description. 
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At this order the self-energies and vertex are given by Eqs. (380), (381) and (384). 
Starting from the three- vertex (384) consider for a moment the vertex resummation for 
the photon leg only, i.e. approximate the fermion-photon vertex by the classical vertex. 
As a consequence, one obtains: 

V^'{x,y;z) - Y''S{x-z)5{z-y) (411) 
-g2 / fH^.z)V^'{x\y';z)A{y\y)fD^,{y,x). 

Using this expression for the photon self-energy (381), by iteration one observes that 
this resums all the ladder diagrams shown in Fig. 26. Here propagator lines correspond 
to self-energy resummed propagators whereas all vertices are given by the classical 
ones. In the context of kinetic equations, relevant for sufficiently homogeneous systems, 
the dominance of this sub-class of diagrams has been discussed in detail in the weak 
coupling limit in the literature (cf. the bibliography at the end of this section). One 
may decompose the contributions to the kinetic equation into 2^2 particle processes, 
such as ee ^ 77 annihilation in the context of QED, and inelastic "1 ^ 2" processes, 
such as the nearly coUinear bremsstrahlung process. For the description of "1 ^ 2" 
processes, once Fourier transformed with respect to the relative coordinates, the gauge 
field propagator in (41 1) is required for space-like momenta. Furthermore, as seen from 
(410), the proper inclusion of nonzero contributions from 2^2 processes requires to 
go beyond the naive on-shell limit. In the context of the evolution equations (390) and 
(391) this can be achieved by employing the following identities: 



F^\x,y) = lim r dz r dz'[pD{x,z)n^F)iz,z')pD{z,y 



dzdz'[D^(x,z)n(^)(z,z')£>A(z,y)]^ 
Po'i^^y) = lini / / dz [pD{x,z)n^p){z,z)pD{z,y 



/oo 
dzdz! [D^(x,z)n(p)(z,zO^A(z , (412) 

-00 

written in terms of the retarded and advanced propagators, D^(jc,y)^^ = ®{x^ — 
y^)pD{x,y)^^ and DA{x,y)^^ = —©{y*^ — x^)pi){x,y)^^ , in order to have an unbounded 
time integration. The above identity follows from a straightforward application of the 
exact evolution equations and using the anti-symmetry property of the photon spectral 
function, p^^(x,y) \jcO=yO = 0. We emphasize that the identity does not hold for an initial 
value problem where the initial time to is finite. Similarly, one finds from (142) for the 
fermion two-point functions using 'fp'^-^\x,y) \xQ=yQ = i5{x — y): 

/oo , . 

Jzdz'AR{x,z)L^^p\{z,z')AA{z,y), 

/oo 
Jzdz'ARix,z)T-\l]{z,z')AAiz,y) , (413) 
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with AR{x,y) = 0(x° -/)p(-^)(jc,);) and Aa(x,j) = -<d{y^ - x^)p^f\x,y). Neglecting 
all derivative terms, i.e. using (399), and the above notation these give:"^^ 



PoiX.k) ~ DR{X,k)n(^p){X,k)DA{X,k), 



(414) 



and equivalently for the fermion two-point functions. Applied to one fermion line in the 
one-loop contribution of Fig. 26, it is straightforward to recover the standard Boltzmann 
equation for 2^2 processes, using the ff{g^) fermion self-energies: 



dV 



FD{X,p)F^f\x,k-pY 



1 



+Lp^{x,pyp^f\x,k-pf 



dV 



+PD{X,p)F^f\x,k-p)^ 



FD{X,p)p^f\x,k-pY 



(415) 



For the Boltzmann equation A^; and A^ are taken to enter the scattering matrix element, 
which is evaluated in (e.g. "hard thermal loops" resummed) equilibrium, whereas all 
other lines are taken to be "on-shell". The contributions from the 1^2 processes may 
be obtained with the help of (414) with the ff{g^) photon self-energies (409). Of course, 
simply adding the contributions from 2^2 processes and 1^2 processes entails the 
problem of double counting since a diagram enters twice. This occurs whenever the 
internal line in a 2 ^ 2 process is kinematically allowed to go on-shell and has to be 
suppressed. 



6.3.2. Discussion 

In view of the generalized fluctuation-dissipation relation (405) employed in the above 
"derivation", one could be tempted to say that for consistency an equivalent relation 
should be valid for the self-energies as well: 



tl^F){X,k) 



^- + nD{X,k) 



n(p)(x,fc). 



(416) 



Such a relation is indeed valid in thermal equilibrium, where all dependence on the 
center coordinate X is lost. Furthermore, the above relation can be shown to be a 
consequence of (405) using the identities (412) in a lowest-order derivative expansion: 
Together with Eq. (414) the above relation for the self-energies is a direct consequence 



As for the spectral function p(X,k) in Eq. (398), the Fourier transform of the retarded and advanced 
propagators includes a factor of — /. 
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of the ansatz (405). However, clearly this is too strong a constraint since the evolution 
equation (403) would become trivial in this case: Eq. (405) and (416) lead to a vanishing 
r.h.s. of the evolution equation for Fd{X, k) and there would be no evolution. 

The above argument is just a manifestation of the well-known fact that the kinetic 
equation is not a self-consistent approximation to the dynamics. The discussion of 
Sec. 6.3.1 takes into account the effect of scattering for the dynamics of effective 
occupation numbers, while keeping the spectrum free-field theory like. In contrast, the 
same scattering does induce a finite width for the spectral function in the self-consistent 
nPI approximation discussed in Sec. 6.2.3 because of a nonvanishing imaginary part of 
the self-energy (cf. also the discussion and explicit solution of a similar Yukawa model 
in Sec. 4.2 and the discussion in Sec. 4.1). 

Though particle number is not well-defined in an interacting relativistic quantum 
field theory in the absence of conserved charges, the concept of time-evolving effective 
particle numbers in an interacting theory is useful in the presence of a clear separation 
of scales. Much progress has been achieved in the quantitative understanding of kinetic 
descriptions in the vicinity of thermal equilibrium for gauge theories at high temperature, 
which is well documented in the literature and for further reading we refer to the 
bibliography below. 

For gauge theories the employed on-shell limit circumvents problems of gauge in- 
variance or subtle aspects of renormalization. We emphasize that renormalization for 
linear symmetries as realized in QED can be treated along similar lines as discussed in 
Sec. 2.2. The generalization to nonabelian gauge theories is, however, technically more 
involved and needs to be further investigated. 

A derivative expansion is typically not valid at early times where the time evolution 
can exhibit a strong dependence on X, and the homogeneity requirement underlying 
kinetic descriptions may only be fulfilled at sufficiently late times. This has been exten- 
sively discussed in the context of scalar and fermionic theories in Sec. 4. Homogeneity 
is certainly realized at late times sufficiently close to the thermal limit, since for thermal 
equilibrium the correlators do strictly not depend on X. Of course, by construction ki- 
netic equations are not meant to discuss the detailed early-time behavior since the initial 
time tQ is sent to the remote past. For practical purposes, in this context one typically 
specifies the initial condition for the effective particle number distribution at some fi- 
nite time and approximates the evolution by the equations with tQ The role of 
finite-time effects has been controversially discussed in the recent literature in the con- 
text of photon production in relativistic plasmas at high temperature. Here a solution of 
the proper initial-time equations as discussed in Sec. 6.2.3 seems mandatory. 
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8. BIBLIOGRAPHICAL NOTES 
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annotated list below, which concentrates on relativistic quantum field theory applications 
related to the content of the text presented above. 
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• The nonperturbative 2PI l/N expansion is derived in: Controlled nonperturbative 
dynamics of quantum fields out of equilibrium, J. Berges, Nucl. Phys. A699 (2002) 
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D63 (2001)096003 
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J. Berges, Sz. Borsanyi and C. Wetterich, Phys. Rev. Lett, in print, http://arXiv:hep- 
ph/0403234 

• Far-from-equilibrium dynamics of macroscopic fields with large fluctuations using 
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large-A^ description has been given in: Analytic and numerical study of preheating 
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